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Chapter 1 

Introduction 



Bose-Einstcin condensation is often considered to be a macroscopic quantum phenomenon. This is 
because bosons occupy the same single-particle ground state below a critical temperature through 
Bose-Einstein condensation to form a macroscopic wave function (order parameter) extending over 
the entire system. As a direct result of the formation of a macroscopic wave function, quantized 
vortices appear in the Bose-condensed system. A quantized vortex is a vortex of inviscid super- 
flow, and any rotational motion of a superfluid is sustained by quantized vortices. A quantized 
vortex is a stable and well-defined topological defect, very different from ordinary vortices in a 
conventional fluid. Hydrodynamics dominated by quantized vortices is called quantum hydrody- 
namics, and turbulence comprised of quantized vortices is known as quantum turbulence (QT). The 
studies of quantized vortices originally began in 1950s using superfluid *He, and much theoretical, 
numerical, and experimental effort has been devoted to the field. Superfluid "^He, discovered in 
1972, presented a system with a variety of quantized vortices characteristic of p-wave superfluids. 
However, quantum hydrodynamics has become very important again in recent years, for two rea- 
sons. The first reason was the appearance of new research activity on QT [1] [^. Before 1990, 
most studies of QT were limited to thermal counterflow in ^Hc in which the normal fluid and the 
superfluid flow oppositely through the injection of heat current. Since counterflow turbulence has 
no classical analog, understanding the relationship between QT and traditional classical turbu- 
lence (CT) has been deserted. However, some experiments without counterflow appeared in the 
1990s, demonstrating the Kolmogorov law of energy spectra which is one of the most important 
statistical laws in turbulence, thus opening a new stage of QT research. The second reason is that 
Bose-Einstein condensation of cold atoms was realized in 1995. This system enables the control 
and direct visualization of condensates, which was impossible in other quantum condensed systems, 
such as superfluid helium and superconductors. The present article will review these recent studies 
on quantized vortices in both superfluid helium and cold atoms. 

The contents of this article are as follows. Section 2 deals with quantized vortices in superfluid 
helium. A brief review of the research history and the dynamics of quantized vortices is followed 
by an important modern topic, QT, which includes the energy spectra issue, dissipation at very 
low temperatures, QT created by vibrating structures, and the visualization of quantized vortices 
and turbulence. In Section 3, which begins with the basics on atomic Bose-Einstein condensates 
(BECs), we review the topics of vortices in single-component, two-component, spinor and dipolar 
BECs, and then discuss quantum turbulence. Section 4 is devoted to a summary and conclusions. 
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Chapter 2 

Quantized vortices in superfiuid 
helium 

Quantum hydrodynamics and quantized vortices have been long studied in superfiuid "^He. This 
section begins by describing the research history and some basic dynamics of quantized vortices, 
and then discusses the modern topics, QT. 

2.1 Research history 

Liquid *He enters a superfiuid state below the A point (T\ —2.17 K) with Bose-Einstein con- 
densation of the ''He atoms. The characteristic phenomena of superfluidity were experimentally 
discovered in the 1930s by Kapitza [3] and Allen et The hydrodynamics of superfiuid helium 

are well described by the two-fluid model proposed by Landau [S] and Tisza[51. According to the 
two-fluid model, the system consists of an inviscid superfiuid (density ps) and a viscous normal 
fluid (density pn) with two independent velocity fields Vg and Vn. The mixing ratio of the two 
fluids depends on the temperature. As the temperature is reduced below the A point, the ratio 
of the superfiuid component increases, and the fluid becomes entirely superfiuid below about 1 K. 
The two-fluid model successfully explained the phenomena of superfluidity, while it was known in 
1940s that superfluidity breaks down when it flows fast [7; and this phenomenon was not explained 
through the two-fkuid model. This was later found to be caused by turbulence of the superfiuid 
component due to random motion of quantized vortices. 

The A transition is closely related to the Bose-Einstein condensation of '*He atoms, as first 
proposed by London [8]. The Bose-condensed system exhibits the macroscopic wave function 
5'(x,t) = \'^{x,t)\e^^^^'*^ as an order parameter. The superfiuid velocity field is given by Vg ~ 
{h/m)\I9, with boson mass m, representing the potential flow. Since the macroscopic wave function 
should be single-valued for the space coordinate x, the circulation F = ^ v • c?£ for an arbitrary 
closed loop in the fluid is quantized by the quantum k = h/m. A vortex with quantized circulation 
is called a quantized vortex. Any rotational motion of a superfiuid is sustained only by quantized 
vortices. 

A quantized vortex is a topological defect characteristic of a Bose-Einstein condensate, and 
is different from a vortex in a classical viscous fluid. First, the circulation is quantized, which is 
contrary to a classical vortex that can have any circulation value. Second, a quantized vortex is a 
vortex of inviscid superflow. Thus, it cannot decay by the viscous diffusion of vorticity that occurs 
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in a classical fluid. Third, the core of a quantized vortex is very thin, of the order of the coherence 
length, which is only a few angstroms in superfluid ^He. Because the vortex core is very thin and 
does not decay by diffusion, it is always possible to identify the position of a quantized vortex 
in the fluid. These properties make a quantized vortex more stable and definite than a classical 
vortex. 

The idea of quantized circulation was first proposed by Onsager, for a series of annular rings 
in a rotating superfluid [9]. Feynman considered that a vortex in a superfluid can take the form 
of a vortex filament, with the quantized circulation k and a core of atomic dimension |10j . Early 
experimental studies on superfluid hydrodynamics focused primarily on thermal counterflow. The 
flow is driven by an injected heat current, and the normal fluid and superfluid flow in opposite 
directions. The superflow was found to become dissipative when the relative velocity between the 
two fluids exceeds a critical value [7]. Gorter and Mellink attributed the dissipation to mutual 
friction between two fluids, and considered the possibility of superfluid turbulence. Feynman 
proposed a turbulent superfluid state consisting of a tangle of quantized vortices [TO]. Hall and 
Vinen performed the experiments of second sound attenuation in rotating ^He, and found that the 
mutual friction arises from interaction between the normal fluid and quantized vortices 112) ; 
second sound refers to entropy wave in which superfluid and normal fluid oscillate oppositely, and 
its propagation and attenuation give the information of the vortex density in the fluid. Vinen 
confirmed Feynman's findings experimentally, by showing that the dissipation arises from mutual 
friction between vortices and the normal fiow [TOl [TH [121 [IS] • Vinen also succeeded in observing 
quantized circulation using vibrating wires in rotating superfluid ^He [171 . Subsequently, many 
experimental studies have examined superfiuid turbulence (ST) in thermal counterflow systems, 
and have revealed a variety of physical phenomena [TO]. Since the dynamics of quantized vortices 
are nonlinear and non-local, it has not been easy to quantitatively understand these observations 
on the basis of vortex dynamics. Schwarz clarified the picture of ST based on tangled vortices 
by numerical simulation of the quantized vortex filament model in the thermal counterflow |19l 
[20] • However, since the thermal counterflow has no analogy in conventional fluid dynamics, this 
study was not helpful in clarifying the relationship between ST and classical turbulence (CT). 
Superfluid turbulence is often called quantum turbulence (QT), which emphasizes the belief that 
it is comprised of quantized vortices. 

Turbulence has long been one of the great mysteries in nature, with discussion dating back to the 
era of Leonardo da Vinci. Turbulence has been intensely studied in a number of fields, but it is still 
far from completely understood. This is primarily because turbulence is a complicated dynamical 
phenomenon with strong nonlinearity. Comparing QT and CT demonstrates the importance of 
studying QT. Turbulence in a classical viscous fiuid appears to be comprised of vortices, as pointed 
out by Da Vinci. However, these vortices are unstable, and appear and disappear repeatedly. 
Moreover, circulation is not conserved, nor is it identical for each vortex. QT consists of a tangle 
of quantized vortices that have the same conserved circulation. Thus, QT may present an easier 
system for study than CT because the elements are more definite and clear, which is confimred at 
relatively large scales by many works described later. However, the situation is not so simple. The 
severe constraint of quantum mechanics will be effective at small scales and lead to some definite 
difference from CT, and it is not so trivial what happens there. 

Based on these considerations, QT research has headed in a new direction since the mid 90s. 
One primary interest has been to understand the relationship between QT and CT [T|. The 
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energy spectrum of fully developed CT is known to obey the Kolmogorov law in an inertial range. 
The energy transfer in an inertial range is believed to be sustained by the Richardson cascade 
process, in which large eddies are broken up self-similarly into smaller eddies. Recent experimental 
and numerical studies have supported the Kolmogorov spectrum, even in QT. Another important 
problem is the dissipative mechanism in QT. At a finite temperature, mutual friction works as a 
dissipative mechanism. However, it is not so easy to understand how the energy cascade is for the 
full range of the scales and what mechanism causes dissipation at very low temperatures, at which 
the normal fluid component is negligible. 



2.2 Dynamics of quantized vortices 

Quantum hydrodynamics, including QT, is reduced to the motion of quantized vortices. Hence, 
understanding the dynamics of quantized vortices is a key issue in quantum hydrodynamics. Two 
formulations are generally available for studying the dynamics of quantized vortices. One is the 
vortex filament model, and the other is the Gross-Pitaevskii (GP) model. We will briefly describe 
these two formulations. 



2.2.1 Vortex filament model 

As described in Section 2. 1, a quantized vortex has quantized circulation. The vortex core is 
extremely thin, usually much smaller than other characteristic scales of vortex motion. These 
properties allow a quantized vortex to be represented as a vortex filament. In classical fluid 
dynamics |21j . the vortex filament model is just a convenient idealization; the vorticity in a realistic 
classical fluid flow rarely takes the form of clearly discrete vorticity filaments. However, the vortex 
filament model is accurate and realistic for a quantized vortex in superfluid helium. 

The vortex filament formulation represents a quantized vortex as a filament passing through 
the fluid, having a definite direction corresponding to its vorticity. Except for the thin core region, 
the superflow velocity field has a classically well-defined meaning, and can be described by ideal 
fluid dynamics. The velocity at a point r due to a filament is given by the Biot-Savart expression: 



K /■ (si — r) X dsi 
Att Jc |si - r|3 



where k is the quantum of circulation. The filament is represented by the parametric form s — s(<;, i) 
with the one-dimensional coordinate along the filament. The vector Si refers to a point on the 
filament, and the integration is taken along the filament. Helmholtz's theorem for a perfect fluid 
states that the vortex moves at the superfluid velocity. Calculating the velocity Vg at a point r = s 
on the filament causes the integral to diverge as Si — >■ s. To avoid this divergence, we separate the 
velocity s of the filament at the point s into two components [T^ : 

The first term is the localized induction field arising from a curved line element acting on itself, 
and and are the lengths of the two adjacent line elements after discretization, separated 
by the point s. The prime denotes differentiation with respect to the arc length The mutually 
perpendicular vectors s', s", and s' x s" are directed along the tangent, the principal normal, and 
the binormal, respectively, at the point s, and their respective magnitudes are 1, R^^, and i?^^, 
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where R is the local radius of curvature. The parameter ao is the cutoff, corresponding to the core 
radius. Thus, the first term represents the tendency to move the local point s in the binormal 
direction with a velocity inversely proportional to R. The second term represents the non-local 
field obtained by integrating the integral of Eq. (|2.f I) along the rest of the filament, except in the 
neighborhood of s. 

Neglecting the non-local terms and replacing Eq. (|2.2I) by s = /3s' x s" is referred to as the local- 
ized induction approximation (LIA). Here, the coefficient /3 is defined hy P — {n / An) \n {c{R) / ao) , 
where c is a constant of order f and {i+i-Y^'^ is replaced by the mean radius of curvature (i?) along 
the length of the filament. This approximation is believed to be effective for analyzing isotropic 
dense tangles due to cancellations between non-local contributions. However, the ILA lacks the 
interaction between vortices, not suitable for the description of a realistic vortex tangle. 

A better understanding of vortices in a real system is obtained when boundaries are included 
in the analysis. For this purpose, a boundary-induced velocity field Vs,b is added to Vg, so that the 
superflow can satisfy the boundary condition of an inviscid flow, that is, the normal component 
of the velocity should disappear at the boundaries. To allow for another, presently unspecified, 
applied field, we include Vg^a- Hence, the total velocity Sq of the vortex filament without dissipation 
is 

At finite temperatures, it is necessary to take into account the mutual friction between the vortex 
core and the normal flow Vn. Including this term, the velocity of s is given by 

s = So + as' X (vn - So) - a's' x [s' x (vn - So)], (2.4) 

where a and a' are temperature-dependent friction coefficients |19| . and So is calculated from Eq. 



The numerical simulation method based on this model has been described in detail elsewhere 
[T9l [20l [22l [23] . A vortex filament is represented by a single string of points separated by a distance 
A<^. The vortex configuration at a given time determines the velocity field in the fiuid, thus moving 
the vortex filaments according to Eqs. (j2.3p and (|2.4p . Vortex reconnection should be properly 
included when simulating vortex dynamics. A numerical study of a classical fluid shows that the 
close interaction of two vortices leads to their reconnection, primarily because of viscous diffusion 
of the vorticity |24| . Schwarz assumed that two vortex filaments reconnect when they come within 
a critical distance of one another, and showed that statistical quantities such as the vortex line 
density were not sensitive to how these reconnections occur [19l[20]. Even after Schwarz 's study, 
it remained unclear as to whether quantized vortices can actually reconnect. However, Koplik 
and Levine directly solved the GP equation to show that two closely quantized vortices reconnect, 
even in an inviscid superfluid [2F. More recent simulations have shown that reconnections are 
accompanied by emissions of sound waves having wavelengths on the order of the healing length 
UllIT]. 

Starting with several remnant vortices under thermal counterfiow, Schwarz numerically studied 
how these vortices developed into a vortex tangle ^U\. The tangle was self-sustained by the com- 
petition between excitation due to the applied fiow and dissipation through mutual friction. The 
numerical results were quantitatively consistent with typical experimental results. Since Schwarz 
used the LIA, he could not obtain the statistical steady state without an artificial procedure. 
Performing the full nonlocal calculation, Adachi et al. succeeded in obtaining a steady state in 
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counterflow turbulence that was consistent with typical observations. This was a significant ac- 
complishment in numerical research |23j . 

Thus, the vortex filament model is very useful for QT, although it cannot describe phenom- 
ena directly related to vortex cores, such as reconnection, nucleation, and annihilation. These 
phenomena can be analyzed only by the GP model. 

2.2.2 The Gross-Pitaevskii model 

In a weakly interacting Bose system, the macroscopic wave function ^{r,t) appears as the order 
parameter of Bose-Einstein condensation, obeying the Gross-Pitaevskii (GP) equation [28l 129) : 

^f^^^^={-^'^'+9m^^t)\'-f^)nr,t). (2.5) 

Here g = Anh^ra/a represents the strength of the interaction characterized by the s-wave scattering 
length a, m is the mass of each particle, and /i is the chemical potential. In = l^fj exp(z6'), 
the squared amplitude is the condensate density, and the gradient of the phase gives the 
superfluid velocity Vg = {h/m)V6 is a frictionless flow of the condensate. This relation causes 
quantized vortices to appear with quantized circulation. The only characteristic scale of the GP 
model is the coherence length, defined by ^ = h/ {y/2mg\'^\), which gives the vortex core size. 

The GP model can explain not only the vortex dynamics but also vortex core phenomena such 
as reconnection and nucleation. However, strictly speaking, the GP equation is not applicable to 
superfluid '^He, which is not a weakly interacting Bose system. The GP model does not feature, 
for example, short-wavelength excitations such as rotons, which are present in superfluid "^He. The 
GP equation is applicable to Bose-Einstein condensation of a dilute atomic Bose gas [30] . 

2.3 Quantum turbulence 

In this section, we review recent developments in the study of QT. "Quantum turbulence" is defined 
as the turbulent state of a quantum fluid, and this term is often used to emphasize the state is 
dominated by the behavior of quantized vortices at low temperatures, with little thermal effect P] . 
From a theoretical point of view, the properties of QT are often discussed in reference to statistical 
quantities such as the energy spectrum, which is important in turbulence. Since quantized vortices 
in QT are definite and stable topological defects, and can be positively identified, the energy 
spectrum of QT is usually related to the dynamics of vortices, such as vortex cascades [31) . 

The energy spectrum and the vortex cascade process are usually divided into two regions in 
wavenumber space. The first region is called the classical region, and exists at wavenumbers below 
the inverse of the mean intervortex spacing. The dynamics of vortices in the classical region are 
dominated by the Richardson cascade, in which large vortices are broken up self-similarly into 
smaller ones, or collective dynamics of aggregated quantized vortices at scales larger than the 
intervortex spacing. Such behavior of vortices supports the analogy of QT to CT, namely, the 
Kolmogorov energy spectrum, which has been confirmed by several theoretical and experimental 
efforts [32) . The second region is called the quantum region, in which vortex dynamics are dom- 
inated by the effects of the quantized circulation, specifically the Kelvin wave cascade of vortices 
[331 134j . which does no appear in CT. The Kelvin wave cascade is also a very important concept 
in understanding the dissipation mechanism of QT at very low temperatures. 
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At finite temperatures, turbulent flow is carried by not only quantized vortices, but also by 
normal fluid. On larger scales, superfluid and normal fluids are coupled together by the mutual 
friction between them, and behave as a classical fluid to show the Kolmogorov energy spectrum. 
On small scales, turbulent flow is dissipated by the viscosity of the normal fluid, and Kelvin waves 
do not exist. 

Another topic in QT is the decay mechanism of quantized vortices P]. The total vortex line 
length, or the vortex line density, can be experimentally observed, and the decay of vortices provides 
important information regarding QT. Since the decay of vortices strongly depends on the initial 
QT state, we can indirectly determine the statistical properties, such as the energy spectrum, of 
vortices in the initial state. 

This section is organized as follows. First, we review the energy spectrum of QT, from the clas- 
sical region to the quantum region, and the connection between these two regions, which remains 
an open question. Next, we discuss the decay of quantized vortices in QT, recent experimental 
developments, QT created by vibrating structures and visualization of quantized vortices. 

2.3.1 Energy spectra in classical turbulence 

To discuss energy spectra of QT, we begin with a brief review of the energy spectrum in CT. The 
dynamics of classical fluids are usually described by the Navier-Stokes equation |55] : 

-^+v-Vv = — VP + j/V^v (2.6) 
ot p 

Here, v = v(x, t) is the fluid velocity, P = P(x, t) is the pressure, p is the fluid density, and v is 
the kinematic viscosity. The flow of this fluid can be characterized by the ratio of the second term 
of the left-hand side of Eq. (|2.6p , hereafter the inertial term, to the second term of the right-hand 
side, hereafter the viscous term. This ratio is the Reynolds number: R — vD /u. Here v and D are 
the characteristic velocity of flow and the characteristic scale, respectively. When v increases to 
the point that the Reynolds number exceeds a critical value, the system enters a turbulent state 
in which the flow is highly complicated, with many eddies. 

The energy spectrum can be obtained from the spatial Fourier transformation of the equal-time 
two-point velocity correlation: 

P(k, i) = ^ y rfx e-^''-^ j dy v(y, t) ■ v(y + x, i). (2.7) 

Here k is the wave number from the Fourier transformation. Kolmogorov proposed the concept 
of a globally steady state of fully developed turbulence [SHI EZ] in which energy is injected into 
the fluid at scales comparable to the system size in the energy-containing range. In the inertial 
range, this energy is transferred to smaller scales without being dissipated. In the inertial range, 
the system is assumed to be locally homogeneous and isotropic, and the angular dependence of 
_E(k, t) in wavenumber space becomes less important. The integration over the angle is defined as 
the energy spectrum: 

E{k, t) - j dcj)^de^ k^F{k, t). (2.8) 

The energy spectrum holds the following relation with the spatial integration of the kinetic energy: 

j dk Eik, t) = j dk |v(k, t)\^^\j d^ v(x, tf = E{t) (2.9) 
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Here, E{t) is the total kinetic energy per unit mass, and v(k, t) is the Fourier transformation of the 
fluid velocity. In the steady state, the averaged value of E{k, t) over t is regarded as the statistical 
value, together with a statistical law known as the Kolmogorov law: 

E{k) = Ce^/S/fc-s/a^ (2.io) 

The energy transferred to smaller scales in the energy-dissipative range is dissipated at the Kol- 
mogorov wave number /ck = (s/'^)^''* through the viscosity of the fluid, at the dissipation rate 
e in Eq. (I2.10p . This spectrum of Eq. (|2.10p is easily derived by assuming that E{k) is locally 
determined by only e and k. The Kolmogorov constant C is a dimensionless parameter of order 
unity. 

The inertial range is thought to be sustained by the self-similar Richardson cascade, in which 
large eddies are broken up into smaller ones with many reconnections 38 . In CT, however, the 
Richardson cascade is not completely understood because it is impossible to definitively identify 
each individual eddy. 

The Kolmogorov law is based on the assumption that the turbulence is locally homogeneous 
and isotropic. However, actual turbulence is neither homogeneous nor isotropic and hence the 
energy spectrum deviates from eq. (|2.10p . This phenomenon is called "intermittency" , which is 
closely related to the coherent dynamics of eddies |35) . 

2.3.2 Energy spectra in quantum turbulence: overall picture 

In QT, it is important to consider the energy spectrum for the following two reasons: First, like CT, 
QT is a complicated dynamical phenomenon with many degrees of freedom, and it is very useful to 
consider the energy spectrum as the statistical law which extracts the important imformations for 
the dynamical properties of the turbulence. Secondly, in QT, unlike in CT, vortices are definite and 
stable topological defects, which enables the study of actual vortex dynamics. In particular, when 
there is no normal fiuid component at near-zero temperatures, any rotational motion is carried 
by quantized vortices and their dynamics are not affected by thermal components, so the inherent 
dynamics of vortices in turbulence are revealed. If QT has an analogy to CT, vortex dynamics 
in QT should show the real vortex dynamics in the Richardson cascade. Therefore, QT is an 
ideal prototype for the study of the statistics of turbulence, such as the relationship between the 
Kolmogorov law in wavenumber space and the Richardson cascade in real space [2] . Whether this 
scenario is true can be confirmed by investigating the energy spectrum and energy transfer in QT. 

As in CT, the energy spectra in QT can also be obtained from the spatial Fourier transformation 
of the equal-time two-point superfluid velocity correlation: Eq. (|2.8p . where v(x, t) is replaced with 
the superfluid velocity Vs(x, t). 

Here, we summarize the overall picture of the energy spectrum of QT at zero temperature 
on the basis of theoretical and numerical studies (Fig. 12. ip [1]. In QT, vortices form a highly 
complicated tangled structure. If the tangle is assumed to be homogeneous and isotropic, there are 
two characteristic length scales: one is the mean intervortex spacing I = L^^/^ with a vortex line 
length density L (the vortex line length per unit volume), and the other is the coherence length ^ 
corresponding to the size of the vortex core. Generally, I is much larger than ^; / 3> ^. Using the 
length scales I and ^, we can define the characteristic wave numbers ki = 2tt/1 and k^ = 27r/^. 

At length scales larger than I, the dynamics of QT are dominated by a tangled structure of many 
vortices. Because vortex dynamics becomes collective at large scales, quantization of the circulation 
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cascade 
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27r/C log k 



Figure 2.1: Overall picture of the energy spectrum of QT at zero temperature. The energy spectrum 
depends on the scale, and its properties change at about the scale of the mean intervortex spacing 
I. When k < ki = 27t/1, a Richardson cascade of quantized vortices transfers energy from large to 
small scales, maintaining the Kolmogorov spectrum E{k) — Ce^^^k^^^^ . When k > ki, energy is 
transferred by the Kelvin-wave cascade, which is a nonlinear interaction between Kelvin waves of 
different wavenumbers. Eventually, energy is dissipated at scales of ^ by the radiation of phonons. 



is not relevant, and the dynamics are similar to those of eddies in CT. Therefore, this region can be 
referred to the classical region. As a result, the energy spectrum E{k) in the range oi k < ki obeys 
the Kolmogorov law (j2.10l) . Vortices in QT sustain a Richardson cascade that transfers energy 
from smaller wavenumbers to larger ones without dissipation. The Richardson cascade can be 
understood as large vortices breaking up into smaller ones in real space. In the Richardson cascade 
process, the key vortex dynamic is reconnection; when two vortex lines approach each other, they 
become locally antiparallel and reconnect [15J [2S1 HZ] • Reconnection causes topological changes of 
the vortex lines in the tangle, formations of distortion waves on the vortex lines (Kelvin waves), 
or fissions of vortex loops through self-reconnections, all of which are responsible for the cascade 
of energy towards smaller scales [351 13S1 131] • The Richardson cascade and the Kolmogorov energy 
spectrum have been confirmed by several numerical studies [40l |4T1 [421 113 Ull US] ■ 

At length scales comparable to I, reconnection is the dominant dynamic. At these scales, two 
vortex lines separated by a distance ~ I reconnect, forming small cusps on the vortex lines, which 
are regarded as the primary source of Kelvin waves in QT. The wavelength of the created Kelvin 
waves is on the order of I. 

At length scales smaller than which is referred to as the quantum region, the Richardson 
cascade is no longer dominant, and the quantized circulation and the motion of each vortex line 
becomes significant [331 US] HT] US] (311 [SD] [SIl [HS] ■ In this range, vortex dynamics are characterized 
by the cascade process of the Kelvin-waves which are formed by reconnection. The nonlinear 
interaction of the Kelvin waves is the origin of the cascade from smaller to larger the wavenumbers. 
The cascade dynamics of the Kelvin waves along a single vortex line has been confirmed, both 
numerically and theoretically. The Kelvin-wave cascade in QT has also been studied. The energy 
spectrum in the quantum region ki < k < k^ is theoretically predicted to obey a Kolmogorov-like 
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power law: E{k) cx fc''. Several values have been proposed for the exponent j], which originates 
from both the configuration of the Kelvin waves and their cascade process. The power-law behavior 
of the energy spectrum in this region has been confirmed numerically [351 H?! SHI Ell ISHl ISI] ■ 

In the region k ^ k^, the Kelvin waves with wavelength of ^ change to elementary excitations, 
such as phonons and rotons, in the primary decay process of QT near zero temperature |34| . 

There is one open question regarding the energy spectrum in the region k ^ ki, which is the 
transitional region between the Richardson cascade and the Kelvin-wave cascade. A theoretical 
study proposed a bottleneck effect connecting the spectrum (|2.10p of the classical region to that 
of the quantum region |52j . Another theoretical prediction was based on vortex reconnection 
dynamics on the scale ~ ^, in which the transitional cascade process was predicted to occur by 
reconnection of the vortex bundle |53] . 

2.3.3 Energy spectra in quantum turbulence: classical region 

The energy spectra in the classical region k < ki has been a major problem in terms of the analogy 
of QT to CT. In this region, the energy spectrum is determined by the collective behavior of many 
vortices, such as the vortex tangle and the aggregated bundle structure at scales larger than I. 
Several numerical studies have calculated the energy spectrum in this region by simulating QT at 
zero temperature. 

There are two formulations for studying the vortex dynamics in QT, as described in Sec. 
12.21 The first is the vortex filament model, in which the time evolution of each element s of a 
vortex filament follows Eq. (j2.4p . At zero temperature there is no mutual friction, and we obtain 
a = a' = in (|2.4p and s = Sq. In this model, the energy spectrum can be calculated directly 
from the configuration of the vortex filaments, but not through the superfluid velocity Vs(f ) of Eq. 
(|2.ip . The vorticity <j^){r) = V x Vs(r) is concentrated on the vortex filament: 

a;(r) = k [ s'(?)5(s(0 - r). (2.11) 



Here, <; is the arc length along the vortex filament. The Fourier transformation of the vorticity 
a;(k) = ik X Vs(k) becomes 



(k) = K y exp[-ik • s(i;)]s'(0. (2.12) 



Assuming fluid incompressibility V • Vs(r) = and Vs(k) = «k x a;(k)/|kp, we obtain the energy 
spectrum as the integrated form over <;: 

E{k) = J dc^i^dOk J d^id^2 exp[-ik • (s(?i) - s(<;2))]s'(^i)s'(?2). (2.13) 

This formula enables us to calculate the energy spectrum directly from the configuration of the 
vortex filament [32] ■ 

Using the vortex filament model, Araki et al. simulated QT with no mutual friction, starting 
from a Taylor Green vortex in a cube of size D = l.Q cm [32]. A cutoff was introduced for the 
smallest vortices, the size of which was comparable to the space resolution = 1.83 x 10~^, which 
was the only dissipation mechanism in the simulation. The energy spectrum E{k) below ki had a 
peak at the smallest fc, and strongly depended on the initial configuration. After about 70 sec, the 
system exhibited a nearly homogeneous and isotropic turbulent state, and lost its dependence on 
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Figure 2.2: Comparison of the energy spectrum (solid line) at < 70 sec after the Taylor Green 
vortex with the Kolmogorov law (j2.10p (dotted line) with C = 1 and e — 1.287 x 10^^ cm^/sec'^ 
[Araki, Tsubota, and Nemirovskii: Phys. Rev. Lett 89 (2002) 145301, reproduced with permission. 
Copyright 2002 by the American Physical Society]. 



the initial configuration. The energy dissipation rate e — ~dE/dt became less dependent on time, 
and the system behaved as a quasi-steady turbulent state. The energy spectrum obtained from Eq. 
(j2.13p agreed with the Kolmogorov law (|2.10p aX k < ki (Fig. 12. 2[) . They also calculated the vortex 
length distribution n(x), where n{x)lS.x is the number of vortex loops, the length of which ranged 
from X to X -\- Ax. In the range from I to Z?, n{x) obeyed the scaling property n{x) oc 2;-i-34±o.i8^ 
which supports the self-similar Richardson cascade from larger to smaller scales. The energy 
transfer rate through the Richardson cascade is equal to the energy dissipation rate e = ~dE/dt 
in Eq. (|2.10p . Furthermore, the energy spectrum E{k) in the region A: < fc; is determined by 
the energy dissipation rate, although the dissipation mechanism is effective only at wavenumbers 
corresponding to the cutoff 27r/A^. In the region of fc > fc;, the energy spectrum deviates from 
the Kolmogorov law, showing a different power-law behavior. Although this behavior may come 
from the velocity field near each vortex line, the simulation was not appropriate to examine the 
spectrum in this region because it did not have sufficient spatial resolution. 

The other formulation is the GP model. The GP equation (j2.5p can be obtained from the 
energy functional: 



H = dyi 



|^|Vvl/(x,i)|^ + ||*(x,t)r-HvI/(x,i)|^ 



(2.14) 



Because the vorticity a;(x, t) = V x Vs(x, t) vanishes everywhere in a singly-connected region of 
the fluid, all rotational flow is carried by quantized vortices. In the core of each vortex, \l/(x, i) 
vanishes so that the circulation ^Vs(x, <) • ds around the core is quantized by h/m. The vortex 
core size is given by the healing length ^ = h/ ^2mgps^ where the superfluid density ps is defined 
as the spatially averaged condensate density |5'(x, t)p. When considering the energy spectrum, it 
should be noted that the hydrodynamics are compressible in the GP model. The total energy of 
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the GP equation per unit mass: 



E{t) I dx 



|l|VM'(x,i)P + f|*(x,i)r (2.15) 
Zm 2 



= / (V|*(x,t)|)2 (2.17) 



can be separated into the interaction energy Eintit), the quantum energy Eq{t), and the kinetic 
energy Ei,in{t); 

E;^,{t) = ^ I dx\^M\' (2.16) 
"27V 

Eunit) - j dx(|vI/(x,t)|V0(x,i))^ (2.18) 

Here, N = J dx |^'(x, t)p is the total number of particles. The kinetic energy can be further 
divided into a compressible part ^'^in(0 '^^'^ ^'^ compressible excitations, and an incompressible 
part -Bkin(t) due to vortices [101 liT] : 

i^kiW = ^^^^ j dx[(|vI/(x,t)|V0(x,t))-f . (2.19) 

Here, [■ ■ - Y denotes the compressible part and V x [• • • J'^ = 0, and [•••]' denotes the incompressible 
part V •[•••]' = 0. Corresponding to each energy, there are several kinds of energy spectra. The 
most important of these is the energy spectrum of the incompressible kinetic energy 

4in(fc, t) = J dcj^kde^ [p'(k, ^)]^ (2.20) 

because it should obey the Kolmogorov law with the Richardson cascade of quantized vortices. 
Here p'(k, i) is the Fourier transformation of the momentum density {|^(x, i)| V6'(x, t)}' 

Starting from a Taylor Green vortex, Nore et al. simulated QT by numerically solving the GP 
model [301 mi- After some time, the initial vortices became tangled, and the calculated spectrum 
obeyed the power- law behavior El^^^{k,t) cx fc^''^*^ When vortices formed a tangle, the exponent 
ri(t) was about 5/3, but this value did not hold for long because the turbulence was decaying with 
the conservation of total energy. Because of energy conservation, the energy of the vortices El^^^{t) 
was transferred to compressible excitations £'^i„(t) through repeated reconnections and disappear- 
ances of small vortex loops with sizes on the order of ^. This thermalization process hindered the 
vortex cascade process through interactions between vortices and compressible excitations, which 
caused El.i^{k,t) to deviate from the Kolmogorov law. 

To avoid this thermalization, Kobayashi and Tsubota proposed a modified GP model, in which a 
dissipation term was introduced to remove the compressible excitations |43[ 144] . The characteristic 
wavelength of the compressible excitations is on the order of ^, so the introduced dissipation was 
set to act only at scales smaller than ^ and to not affect vortices. The resulting GP equation in 
wavenumber space becomes: 

- 7(fc)] = K:- - *(k, t) 



dt V 2m 



(27r) 



. ' (2.21) 

^ / dkidk2*(ki,t)**(k2,t)*(k-ki+k2,0. 
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Here, ^(k, t) is the Fourier transformation of 4'(x, t). The dissipation term has the form j{k) — 
7o^(A:— 27r/^) with the step function 9, being effective only at fc > k^. Introduction of 7(fc) conserves 
neither the energy E nor the number of particles N. However, when studying the hydrodynamics 
of turbulence, it is realistic to assume that the number of particles is conserved. Hence, the time 
dependence of the chemical potential was introduced to conserve the total number of particles 
TV. By numerically analyzing Eq. (j2.2ip . they confirmed the Kolmogorov spectrum in two kinds 
of QT. 

The first analysis was of decaying turbulence [43]. To obtain a turbulent state, an initial 
configuration was set to have a uniform density |'I'(x, t = 0)p = 1 and the phase had a random 
spatial distribution. Because the initial superfluid velocity Vs(x, t = 0) = {h/m)'S/9{'x,t ~ 0) was 
also random, the initial wave function was dynamically unstable and soon produced fully developed 
turbulence with many quantized vortex loops. They confirmed that only the compressible kinetic 
energy was decreased by the dissipation term, and that the incompressible kinetic energy 
dominated the total kinetic energy i^kinl^), demonstrating that the thermalization was effectively 
suppressed. The obtained energy spectrum of the incompressible kinetic energy was consistent with 
the Kolmogorov law in the period during which the vortex tangle was formed, and this consistency 
lasted longer than it did without the dissipation term. In this period, the energy dissipation 
rate e — —dEl^^{t)/dt became nearly constant, and the system can be considered to occupy a 
quasi-steady turbulent state with a stationary Kolmogorov spectrum of e^^^^k^^^^. 




Figure 2.3: Incompressible energy spectrum El^^^{k) obtained by Kobayashi and Tsubota [H]. The 
plotted points are from an ensemble average of 50 randomly selected states. The solid line is the 
Kolmogorov law [Kobayashi and Tsubota: J. Phys. Soc. Jpn. 74 (2005) 3248, reproduced with 
permission] . 

The second analysis was of steady turbulence, with not only energy dissipation but also large- 
scale energy injection [44^ . The energy injection was accomplished by introducing the stochastic 
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external potential ^(x, t) into the GP equation (j2.21D . V{x, t) has a Gaussian two-point correlation: 
(F(x,t)t/(x',t')> ^V^o'exp 

where Xq and Tq are the characteristic spatial and time scales of the potential. By introducing 
T^(x, t), energy was injected at the spatial scale of Xq, and quantized vortices of radius Xq were 
effectively nucleated. Starting from the uniform state \E'(x, t = 0) = 1, they developed the GP 
equation (I2.2ip with the potential (j2.22p . and obtained a steady QT in which E{t), E\^in{t), £'kin(*)' 
and El.-^^{t) were nearly constant. In the simulation, the incompressible kinetic energy El.-^^{t) was 
always dominant in the total kinetic energy -EkinC^); the introduced potential contributed to the 
nucleation of vortices rather than that of compressible waves. The obtained energy spectrum was 
consistent with the Kolmogorov law in the range 2tt/Xo < which is regarded as the inertial 
range of Richardson cascade of vortices. (Fig. 12. 3p . They also calculated the flux Il{k,t) of the 
incompressible kinetic energy from smaller to larger wavenumbers: 

n(fc,i) = ^^^ J dxLfc[{p(x,t).Vv,(x,t)r]L,[{p(x,i)r], (2.23) 

where Lj~ is the operator for the low-pass filter: 

ifc[^W] = 7^/ dk [ d^' e^<---'\si^), (2.24) 

for an arbitrary function s(x). Il{k,t) can be obtained from the scale-by-scale energy budget 
equation derived from the GP equation (|2.2ip . Il{k,t) was nearly constant in the inertial range, 
and almost the same as the energy dissipation rate e = —dEl^^{t)/dt after switching off the moving 
random potential. This confirmed the picture of the inertial range and the Richardson cascade of 
quantized vortices in QT; the incompressible kinetic energy steadily flows in wavenumber space 
through the Richardson cascade at the constant ratell, and finally dissipates at the rate e ~ H. 

In the above three numerical studies, the system size was not so large that the inertial range was 
less than one order in wavenumber space. Furthermore, the mean intervortex spacing I was close 
to the healing length ^, and was too short to study the Kelvin wave cascade. To obtain the energy 
spectrum of a wider range of wavenumber space, Yepez et al. performed a large-scale simulation 
of the GP model, the size of which was (11.25)'^ times larger than the above three simulations |45) . 
They also used a highly accurate numerical code for the time development; the unitary quantum 
lattice gas algorithm, in which the time development of the order parameter ^'(x, t) in Eq. (12. 5p 
could be performed as the unitary time evolution of qubits on a cubic lattice. Starting from the 
initial state with 12 straight vortices, which consisted of three groups of four vortices aligned along 
X, y, and z axes, they obtained QT with highly tangled vortices. In this simulation, two length 
scales were defined; the inner radius of a vortex core f , and its outer radius tt^. On an grid, 
the wavenumber corresponding to the core's inner radius was dinner — (a/3/2)L/^, while the outer 
radius wavenumber was fcoutcr — dinner /tt. They found that the incompressible kinetic energy 
spectrum (|2.20p had three distinct power-law k^" regions that ranged from the classical turbulent 
regime of Kolmogorov a = —5/3 at large scales k < /coutcr to the quantum Kelvin- wave cascades 
a = 3 at small scales k > fcinncr- There was a semiclassical region 6.34 < a < 7.11 connecting 
the Kolmogorov and Kelvin- wave spectra /coutcr < k < fcjnnor- Compared to previous simulations, 
this simulation supplied the Kolmogorov spectrum over a much wider inertial range, with about 2 



^'\2 



\' {t-t'Y 



2Xi 



2Ti 



(2.22) 
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orders in wavenumber space, and the Kelvin- wave spectrum in the quantum region. Although they 



in this region is smaller than the vortex core size. The one-dimensional picture of Kelvin wave is 
relevant in the scales much larger than the vortex core, and a nature of short waves propagating in 
three-dimensional space is considered to be dominant in the length scale smaller than the vortex 
core size. They also did not discuss the mean intervortex spacing I, which is the key length scale to 
distinguish the classical and quantum resions. Therefore, it is not so clear that the spectrum 
reflects the Kelvin-wave cascade. 

It should be noted that the Richardson cascade of quantized vortices is genuine in QT, and 
does not completely imitate that of CT. The Kolmogorov constant C in Eq. (|2.10p in QT is not 
necessarily equal to that of CT. The obtained Kolmogorov constant was C ~ 0.7 for the vortex 
filament model [H], and C ~ 0.32 ^ and 0.55 for the decaying and steady QT in the GP 
model, which was smaller than that of CT, in which C ~ 1.5. These smaller Kolmogorov constants 
may be characteristic of QT. 

2.3.4 Energy spectra in quantum turbulence: quantum region 

In the region of A; > fc;, the picture of aggregated vortices is no longer effective, and the motion 
of each vortex line becomes essential. The most probable dominant dynamic of vortices is Kelvin 
waves, which originate from distortion waves on the vortex lines after their reconnection. A Kelvin 
wave is a transverse, circularly-polarized wave motion, with the approximate dispersion relation 
for a rectilinear vortex: 



with a dimensionless constant c ~ 1. A: is the wavenumber of the Kelvin wave, and was different from 
that used for the energy spectrum. Kelvin waves were theoretically proposed to exist in inviscid 
fluids [54], and were first observed by inducing torsional oscillations in a rotating superfluid '^He 



At finite temperatures where there is a significant fraction of normal fluid, Kelvin waves are 
damped by mutual friction. On the other hand, at very low temperatures they can be damped only 
by the radiation of phonons. Vinen estimated the rate of radiation, and found that it is extremely 
low unless the frequency is very high, typically on the order of 4 GHz for superfluid '*He |34| . In 
these circumstances, low-frequency Kelvin waves can lose energy only by nonlinear coupling to 
waves of a different frequency, which is the Kelvin wave cascade from small to large wavenumbers. 

An earlier numerical work showed the nonlinear interaction of Kelvin waves in the framework of 
the vortex filament model. Samuels and Donnelly performed a simulation of a full Biot-Savart law 
(j2.2p and the local induction approximation (LIA), in which the second Biot-Savart term in (|2.2p 
was neglected, and found sideband instability in both cases above a critical amplitude of the Kelvin 
wave [13]: Aq > X/ (27rn), where Aq is the amplitude of the main helical wave, A is the wavelength, 
and n is the number of half waves on the vortex. The occurrence of sideband instability indicates 
the mixing of Kelvin waves of the main wavelength with those of smaller wavelengths via nonlinear 
interactions, and the transfer of energy from small to large wavenumbers. 

Although the Kelvin-wave cascade seems to be a very important mechanism in QT at scales 
smaller than I, it is non-trivial problem about the actual cascade process. Vinen et al. performed 
a numerical simulation of a Kelvin wave excited along a single vortex line using the vortex filament 



related the k spectrum at small scales k > ki 



^mncr 



to the Kelvin-wave cascade, the length scale 




(2.25) 



[35] [55]. 
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model |47j . They considered a model system in which the helium was contained in a space between 
two parallel sheets separated by a distance £b = 1 cm. A single, initially rectilinear (along the z 
axis), vortex was stretched between the two sheets. The allowed wavenumbers of the Kelvin waves 
were fc = 27rn/^B, where n is a positive integer. A Kelvin wave with a small wavenumber of no 
was continuously driven, and all modes with n exceeding a large critical value ric were strongly 
damped. The simulation was based on the full Biot-Savart law, and the force driving one mode 
was of the form V pnsvwikfiZ — wo)j where /cq — 27rno/^B, P is the density of helium, and wq is 
related to the frequency by the dispersion relation (j2.25l) . Damping was applied by a periodic 
smoothing process. Starting from a straight vortex line, the total length of the vortex line evolved 
to reach a steady average value after application of the driving force, which suggests the existence 
of a steady state. In the typical simulation, with V = 2.5 x 10^^ cm s^^, feg = IOtt cm~^, and 
fee = ^-Kncl^-Q — 1207r cm~^, they calculated the root mean square amplitudes nk{t) = (CfcCfe)^^^ of 
the Fourier components of the displacement of the vortex. Initially, only the mode that resonated 
with the drive was excited. However, as time passed, nonlinear interactions led to the excitation of 
all modes. The spectrum nfc(t) also reached a steady state, and in this stage energy was injected 
at a certain rate at wavenumber fco and dissipated at the same rate at wavenumber h^. For large 
fc, where the modes nearly formed a continuum, the steady state was observed to have, to a good 
approximation, a spectrum of the form 



with the dimensionless parameter A of order unity. They continued their simulations, changing 
the driving amplitude V and the drive wavenumber fco to show that there was no effect on the 
steady state, which suggests that the scaling property of the spectrum (|2.26p is universal. The 
mean energy per unit length of a vortex in mode k is related to by the equation [34) 



Prom Eqs. (|2.26p and (|2.27p . the energy spectrum of the single Kelvin wave becomes EK{k) = 



When two vortices reconnect, they twist to become locally antiparallel at the reconnection 
point and create small cusps or kinks after the reconnection, which were initially confirmed nu- 
merically by Schwarz [Tni HO] ^ and later by Tsubota et al. in detail [32] using the vortex- filament 
model. Svistunov suggested that the relaxation process of these cusps or kinks causes the emission 
of Kelvin waves, and plays an important role in the decay of QT at low temperatures |33j . Fol- 
lowing Svistunov, Vinen analyzed the energy spectrum of the Kelvin-wave cascade by introducing 
a "smoothed" length of vortex line per unit volume after all the Kelvin waves were removed, and 
considered EK{k)dk, the energy per unit length of the smoothed vortex lines associated with Kelvin 
waves in the range k to k + dk ^3^. By dimensional analysis, E^ik) was estimated as 



Ilk — ^ 



(2.26) 



EK{k) = CKk^Uk, 



(2.27) 



where ex is an effective energy per unit length of vortex, given by 




(2.28) 



Exik) = ApK^k 



(2.29) 



with a constant A. 
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Kivotides et al. numerically confirmed the generation of Kelvin waves through reconnections 
using the vortex filament model |39j . In their simulation, four vortex rings were placed symmet- 
rically on opposite sides and oriented so that they all moved toward the center. The four rings 
approached each other and underwent reconnections which induced cusp relaxation and the gener- 
ation of large-amplitude Kelvin waves. They calculated the energy spectrum E{k) defined by Eq. 
(|2.8p for the superfluid velocity, and found that E{k) developed approximately a form after 
the reconnections occurred and the reconnected vortex lines began to move farther apart. This 
supports the nonlinear transport of energy between different wavenumbers. Considering that the 
fluctuations of the velocity field were induced by the Kelvin waves on the filament, this result is 
consistent with Vinen's analysis of Eq. (j2.29p with EK{k) ^ E{k). 

Kozik and Svistunov analyzed the Kelvin- wave cascade using weak-turbulence theory [48] . They 
employed the Hamiltonian representation of the vortex line motion; a Kelvin wave was excited along 
the z direction and the position of the vortex line was specified in the parametric form x — a;(z), 
y — y{z). In terms of the complex canonical variable w{z,t) = x{z,t) + iy{z,t), the Biot-Savart 
dynamics equation (j2.2l) . acquires the Hamiltonian form iw = SH[w,w*]/6w* with 

= ^ [dz.dz. l-HReK-(.rX(..)] 
47ri V(zi-Z2)2 + |ii;(zi)-«;(z2)P 

where w'{z) = dw{z)/dz. Under the assumption that the amplitude of the Kelvin wave is small 
compared to its wavelength, the Hamiltonian (|2.30p can be expanded in powers of — 
w{z2)\/\zi — Z2I: H — Eq + Ho + Hi + ■ ■ ■ [Eo is just a number). Hq can be diagonalized by 
the Fourier transformation wt — J dz w{z)e^^^^- to 

^^0 = ^ X! ^kw'lwk, (2.31) 
fe 

where a;^ is the dispersion given in Eq. (|2.25p . By introducing kelvons, which are the quanta 
Wk for the amplitude of Kelvin waves Wk, the Hamiltonian can be canonically quantized with the 
kelvon annihilation operator d^ — \/ K,p/2hwk- Because of the conservation of momentum and 
energy, the two-kelvon scattering channel is suppressed and three-kelvon scattering becomes the 
dominant process. The effective vertex ¥^2'^ for the three-kelvon scattering process consists of 
two parts; a two-kelvon vertex in the second order of perturbation associated with Hi, and the 
bare three-kelvon vertex associated with H2. In the classical- field limit, the kinetic equation of the 
averaged kelvon occupation number Uk — {a\ak) over the statistical ensemble is 

hi = 2167r ^ |<2'3T<5(Ac.)5(Afc)(/i;5%3 - ftf^). (2.32) 

with Afc — fci + fc2 + fca — fc4 — fcs — fcg. Aw = wi + a;2 + W3 — a;4 — tJs — ajg, and fi'2'3 — 
"■i"-2"-3("-4"-5 + n^nQ + ngn^). Here, (1, • • • ,6) denotes the set of the wave number (fci, • • • , kg). 
The kinetic equation (|2.32p supports the energy cascade, if two conditions are satisfied. The first 
is that kinetic time grows progressively smaller in the limit of large wavenumbers, which can be 
checked by a dimensional estimate. The second is that the collision term is local in wavenumber 
space, which has been verified numerically. By the dimensional analysis, Eq. (j2.32p yields 

hk (XLu-^^nlk^^, (2.33) 

with fci ~ • • • ^ fcg ^ fc and \ V\ ^ k^. The energy flux 9k per unit vortex line length at momentum k 
is defined as 9k — L^^ J2k'<k'^k'hk' , where L is the system size. This implies 9k ~ khkOJk ~ n^.k^'^ 
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with Eq. (j2.33p . The cascade requirement that 6^ is k independent, without energy dissipation 
in the inertial range, leads to the spectrum: rik oc k^^"^/^ and cx ^^.w/cn/j ^ X^fc 
To investigate the above spectrum in more detail, Kozik and Svistunov performed a numerical 
simulation of the full Biot-Savart law of Equation (|2.2[) over a wide range of scale 49 . Starting 
from a single vortex line with the initial distribution Uk oc fc~^, they found that the waves converted 
into a new power-law distribution <x k~^'' . 

For the spectrum of Kelvin waves, Nazarenko presented a nonlinear differential equation model, 
pointing out that turbulence displays a dual cascade behavior of both the direct energy cascade 
and an inverse cascade of wave action, which correspond to two constants of motion: 

and 



s X ds (2.35) 

for the full Biot-Savart equation (|2.2p |50) . The nonlinear differential equation was constructed to 
preserve the main scaling of the original closure cx fc"^''/^ as 



where C is a dimensionless constant, lo = a;(fc) = nk"^ /An is the approximate dispersion of the 
Kelvin wave. This equation preserves the energy E = J duj nuj^^'^, and the wave action N = 
J duj mjj~^l'^. Equation (j2.36p has both the direct cascade solution cx fc~^^/^ and the inverse 
cascade solution cx k~^ . It also has the same thermodynamic Ray leigh- Jeans solutions n = 
T I (w -I- [s) with constants representing temperature T and chemical potential \x. This model was 
developed to include the effect of generation of Kelvin waves by reconnection, dissipation caused 
by mutual friction, and phonon radiation. 

Boffetta et al. suggested a simpler model for the Kelvin-wave cascade than that of Svistunov el 
al. [35j, using a truncated expansion of the LIA [ST]. As is the case in Eq. (|2.30p . a Hamiltonian 
form for the LIA can be obtained as 

U{w\ = 2^ In (^0 y dz v/rTRW - 2/3L[H, (2.37) 

where £ is a length on the order of the curvature radius, and /? = (K/47r) ln(£/^). The Hamiltonian 
is proportional to the vortex length = J dz ^\ + |w'(z)p. The equation of motion is 

(2.38) 




with /3 — 1/2. As a consequence of the invariance under phase transformation, the total wave action 
N\w\ = J dz jwp is also conserved. In the framework of the LIA, the system becomes integrable 
with infinite numbers of conserved quantities, so that the energy and the wave action cannot 
cascade. The integrability is broken by considering a truncated expansion of the Hamiltonian in 
power of wave amplitude w'{z) such as 

ilexpN ^H^ + Hi+H2= / 1 + - ^ . (2.39) 
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Taking into account only the six-wave resonant condition, and taking the Taylor expansion of w{k, t) 
around w{k,t) = c(fc,0) {w{k,t) = J dze~'^'^^w{z,t) is the Fourier transformation of w{z,t)), one 
can obtain the resulting Hamiltonian expressed in and the equation of motion for = (jcfcp), 



rffcl23456 C'l23456'5456CiC2'^3C4C5C6, 



•hk — IStt / dh 



<''23456 |C'fe23456 



|2 rfe23 c/, ,fc^34\ r- 
'^456 ('^456 )/' 



k23i\ 



fe23456, 



where C123456 is the interaction coefficient, and 
/fc23456 = nkn2n3n4n5ne — 



1 

"-2 



1 

"3 



1 



1 

n5 



1 

ne 



5(^456) = S{ujk + UI2 + UJ3 - uji ~ UJ5 - We)- 



(2.40) 
(2.41) 

(2.42) 
(2.43) 



The kinetic equation (|2.41l) is the almost same as that of the Biot-Savart formulation (|2.32p by Svis- 
tunov, and a simple dimensional analysis of Eq. p. 411) gives hk oc Besides the equilibrium 

Rayleigh- Jeans solution Uk —T / {cok + fJ-), there are non-equilibrium steady state solutions of the ki- 
netic equation (I2.4ip in wave turbulence theory, which rely on a constant flux in some inertial range 
and are known as Kolmogorov-Zakharov solutions |57j . Using a dimensional analysis, one can ob- 



tain the energy flux Ilj^. 

existe] 

to the spectra: 



/ dk'hk'iOk' ~ k^'^n\ and the wave action flux n^,^-* = j dk'hk' '■^ 
Requiring the existences of two ranges of scales in which III^^ and n[,^^ are /c-independent leads 



uk-k-^y^ 

Hk ~ k~^, 



(2.44a) 
(2.44b) 



respectively. As in the analysis by Nazarenko [5D], the first spectrum is expected to be the direct 
cascade of energy fiowing to large k, and the second is expected to be the inverse cascade of wave 
action fiowing to small k. Finally, we mention the energy spectra estimated from the spectra in 
(12.44p . By using Eq. (|2.27p and the relation E]<^{k) ^ E{k) proposed by Kivotides et ai, we can 
obtain E{k) oc for the direct cascade and E{k) oc k~^ for the inverse cascade. 

Yepez et al. numerically investigated the energy spectrum in the quantum region |45] using the 
Gross-Pitaevskii model. The obtained spectrum E{k) oc fc^'^ is inconsistent with those of the above 
analytical works, and the energy spectra in the quantum region remains as an open question. As 
described in the previous subsection, however, the relationship between the Kelvin-wave spectrum 
and the mean intervortex spacing I was not mentioned, and it is not clear whether the obtained 
spectrum can appropriately be regarded as the spectrum of the Kelvin-wave cascade. 



2.3.5 Energy spectra in classical-quantum crossover 

As discussed in the previous sections, there are two different types of energy spectra in the classical 
[k < ki) and quantum {ki < k < k^) regions. An important question arises: How do these two 
energy spectra connect to each other at the length scale 17 Although there are several theoretical 
and numerical reports on this region, consistency among this work has not yet been obtained, and 
the problem remains controversial. In analysis of classical-quantum crossover, A = ln(^/^) appears 
to be an important parameter. In typical ^He experiments, A is about 15. 

L'vov et al. suggested a bottleneck crossover between the two regions ^2]. For a given mean 
intervortex spacing I, the mean vorticity in the system (|w|) is given by (|w|) ~ kI^^. The energy 
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spectrum in the classical region is of the Kolmogorov form (Eq. (I2.10p '): Ec\{k) ~ e^/'^k^^/'^ ^ 
where the Kolmogorov constant C is approximately unity. Assuming that is dominated by 
the classical-quantum crossover scale, we obtain 



f ~ PE.iik) ~ e2/3^-4/3_ (2.45) 



In the quantum region, the reformulated kelvon occupation number from the analysis by Kozik 
and Svistunov is rik ~ (Ze)^/^^^/^^"^^/^ from dimensional analysis, and the corresponding energy 
spectrum becomes 
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E^^{k)^A k-y\ (2.46) 



which is consistent with Eq. (j2.44a|) . This equation also gives the mean vorticity as 

(M>^ - - A f^V' . (2.47) 



If the energy flux e is the same in both the classical and the quantum regions, the ratio between 
Eci{l/l) and Eq^{l/l) at the crossover of A; ~ l/l is Equ{l/l)/E^i{l/l) ~ A^o/^ > 1. This large 
mismatch indicates that the energy flux carried by classical hydrodynamic turbulence cannot fully 
propagate through the crossover region, and that larger scale hydrodynamic motion will increase in 
energy up to the level £^qu(l/0; this is called the bottleneck effect (Fig. 12. 4p . To get a qualitative 
picture of the bottleneck, L'vov et al. used the warm cascade solutions following from the model: 

for the energy flux of hydrodynamic turbulence [SS] . Here Fk is the three-dimensional spectrum 
of turbulence. For a constant energy flux, the solution of Eq. (|2.48p is 



Fk = 




2/3 

(2.49) 



where p is the fluid density. While (|2.49p coincides with the Kolmogorov energy spectrum at low 
fc, the large k range comprises a thermalized portion of the spectrum with equipartition of energy, 
characterized by an effective temperature T. 

The analysis by L'vov et al. was based on the assumption that the coarse-grained macroscopic 
description of quantized vorticity is effective down to the scale of I. Based on the LIA, Kozik 
and Svistunov suggested a different picture for the crossover region, in which the locally-induced 
motion of the vortex lines emerges at the scale of rp ^ A^/'^Z, and the crossover range is divided 
into three subranges, r^^ < k < \^^, A^^ < k < A^T"^, and A~^ < k < A~^, where Ab A^/**/, 
Ac l/A^/^, and A* //A^/^ (Fig. [13]) 53 . In the first region, r^^ < k < Ab, polarized 
vortex lines are organized in bundles and reconnect with other bundles to form Kelvin waves with 
amplitude bk ~ r^^fc"^, where bk is defined as bj, ^ ^q^ki^q'^i) ^ ^^^X^q^fc"? ~ ^''^k- 
Here, is the kelvon annihilation operator and nq is the kelvon occupation number. In the second 
region, X^^ < k < A~^, the cascade is supported by nearest-neighbor reconnections in a bundle, 
and bk ~ l{Xhk)^^^^. In the third range, A^T^ < k < A^^, the cascade is driven by self-reconnection 
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Figure 2.4: The energy spectra E(k) in the classical, k < l/l, and quantum k > l/l ranges, 
proposed by L'vov et al. |52] . The two straight solid lines in the classical range indicate pure 
Kolmogorov scaling E(k) cx k~^^^ and pure thermodynamic scaling E{k) oc fc^. In the quantum 
range, the dashed line indicates the Kelvin wave cascade spectrum (slope —7/5, whereas the dash- 
dotted line marks the spectrum corresponding to the noncascading part of the vortex tangle energy 
(slope -1) [L'vov, Nazarenko, and Rudenko: Phys. Rev. B 76 (2007) 024520, reproduced with 
permission. Copyright 2007 by the American Physical Society]. 

of vortex lines, giving bk ^ k~^. The quantum region for the Kelvin- wave cascade of a single vortex 
line begins from A"-'^, giving bk c>c k~^^^. The Kelvin- wave spectrum bk smoothly connects these 
ranges, and there is no bottleneck effect in the model. Although they emphasized that the energy 
spectrum E{k) is practically meaningful only in the classical region, we can estimate E{k) from 
their model: E{k) ^ fc~^ in ^ < fc < E{k) fc° in A^^ < fc < A^:^, and E{k) ^ k"^ in 
A-i < k < A-i. 

The theoretical works of both L'vov et al. and Kozik et al. were based on the idea that 
quantized vortices are locally polarized in a tangle, forming bundles. However, the nature of 
their polarization is not clear. Furthermore, the energy spectra estimated from their model were 
inconsistent with those given from the numerical work of the GP model by Yepez et al [45] which 
also remained a problem in this wavenumber region as we commented in the subsection l2.3.3l More 
quantitative analysis and detailed numerical work is needed to clarify the nature of the crossover 
region. 

2.3.6 Energy spectra in quantum turbulence at finite temperatures 

Experimental studies measured the energy spectrum of QT at finite temperatures, and supported 
the Kolmogorov spectrum directly or indirectly [SHI HOI ISB ISH IMl [Ml ISi] . These experiments were 
also consistent, in the sense that they showed similarities between QT and CT. Vinen theoretically 
considered this similarity, and proposed that the superfluid and the normal fluid are likely to be 
coupled together by mutual friction at scales larger than the intervortex spacing I, and would thus 
behave like a classical fluid and show a Kolmogorov energy spectrum where the mutual friction 
does not cause dissipation |66l . This idea was confirmed by Kivotides et al. through numerical 
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Figure 2.5: Spectrum of Kelvin waves in the quantized regime proposed by Kozik et al. |53) . The 
inertial range consists of a chain of cascades driven by different mechanisms: (1) reconnections 
of vortex-hne bundles, (2) reconnections between nearest-neighbor vortex lines in a bundle, (3) 
self-reconnections on single vortex lines, and (4) nonlinear dynamics of single vortex lines with- 
out reconnections [Kozik and Svistunov: Phys. Rev. B 77 (2008) 060502 (R), reproduced with 
permission. Copyright 2008 by the American Physical Society]. 

simulation of coupled dynamics of a vortex filament and a normal fluid [67], and by L'vov et al. 
through theoretical analysis of the two fluid equations [68 . Although they also showed a decoupling 
of the two fluids at small scales, this remains a controversial topic. 

2.3.7 Experimental study of quantum turbulence 

Early experimental studies of QT focused on thermal counterflow, in which the normal fluid and 
superfluid flow in opposite directions |69| . However, as thermal counterflow has no analogy with 
conventional fluid dynamics, it has not enabled an understanding of the relationship between QT 
and CT. In the mid-90s, QT experiments were performed that did not involve thermal counterflow. 
Maurer and Tabeling studied QT of superfluid ^He that was produced in a cylinder 8 cm in 
diameter and 20 cm high and driven by two counter-rotating disks [60J. They observed local 
pressure fluctuations, and converted these into an energy spectrum. The experiments were done at 
three temperatures: 2.3 K {ps/p = 0), 2.08 K {ps/p ~ 0.05), and 1.4 K {ps/p = 0.9), which are both 
above and below the A point. In all Kolmogorov energy spectrum was observed (Fig. 12. 6p . 

Above the A point, the obtained spectrum is reasonable because the system is a classical viscous 
fluid. Below the A point, however, it is debatable whether the obtained spectra are consistent 
with the Kolmogorov law and independent of the ratio between ps and pn, namely, the dissipative 
mechanism. Their observations were understood on the basis of the idea that the two fluids were 
probably coupled together by mutual friction, and behaved as a one-component fluid |66j . 

Afterwards, a series of experiments with superfluid "^He above 1 K were performed by the 
Oregon group [59l |61] [62l |63l [64l |65]. The helium was contained in a channel with a square (1 cm 
X 1 cm) cross section, along which a grid was pulled at a constant velocity. This method of creating 
QT is similar to a method of generating homogeneous, isotropic turbulence in a classical fluid [35) . 
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Figure 2.6: Energy spectra obtained using counter-rotating disks at (a) 2.3 K, (b) 2.08 K, and (c) 

1.4 K m. 



A pair of second-sound transducers placed in the walls of the channel detected the appearance of 
a vortex tangle by second-sound attenuation [32] . The decay of the vorticity of the tangle behind 
the grid was observed. To obtain the energy spectrum from this result, the authors made some 
assumptions. In a fully developed CT, the energy dissipation rate becomes e = t^(w^), where (w^) 
is the mean square vorticity |70) . A similar formulation could be satisfied in fully developed QT 
above 1 K, in which k^L^ is a measure of (uj^). They assumed that the dissipation rate in the 
experiment became 

e = v'K'^L^, (2.50) 

with an effective kinematic viscosity v' . They also assumed that the fiows of the superfluid and 
normal fluid were coupled together to behave as a one-component fluid by mutual friction at scales 
larger than and Eq. (|2.50p could be directly connected with the observation of second-sound 
attenuation by choosing a suitable value of v' as a function of temperature [64] . The result suggests 
that kL decays as kL oc which is indirectly consistent with the Kolmogorov energy spectrum 

in QT. 

The relationship between kL oc t^'^/^ and the Kolmogorov energy spectrum was reviewed by 
Vinen and Niemela [32]. When we assume that the flow motion in the inertial range < k < 
with <§; kd dominates the overall fluid dynamics and gives a Kolmogorov energy spectrum 
(|2.10p , the total energy can be approximated by 

E C£2/3fc-5/3 ^ ^Ce2/3D2/3. (2.51) 

If the turbulence decays slowly, the dissipation rate e can be written as 

s = -^ = -Ce-'/'D''^'^. (2.52) 
dt dt ^ ' 
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The solution of this differential equation is £ = 27C^D^ /{t + t^Y ^ with a constant t^. This result 
and Eq. ((^3(1)) give 



The behavior L ^ i^^/^ for large t suggests that the Kolmogorov energy spectrum in QT is 
independent of the value of v' , namely the temperature. 

This type of decay L ~ was also observed in turbulence by an impulsive spin down for 

superfluid *He |71j . in turbulence by injecting negative ions into superfluid ^He |72| . and for grid 
turbulence in superfluid ^He-B 73 , with little normal fluid at low temperatures, which supports 
the classical analogue of QT with no normal fluid component. 

The Helsinki group experimentally studied vortex dynamics of propagation into a region of 
vortex- free flow in a rotating superfluid ^He-B |74) . They measured the velocity of the vortex front 
toward the metastable region, and determined the rate of dissipation as a function of temperature. 
The results showed a transition from laminar through quasiclassical turbulent to quantum turbulent 
flow with decreasing temperature. Below 0.25Tc with a superfluid critical temperature Tc, there 
was a peculiar decrease of dissipation. As one possibility, the authors suggested that the energy 
flux toward small scales propagates to the scale Z, and vortex discreteness and quantization effects 
become important. The dominant part of the energy loss was the Kelvin- wave cascade below the 
scale I. As proposed by L'vov et al., Kelvin waves are much less efficient in downscale energy 
transfer than classical turbulence, which leads to a bottleneck accumulation of kinetic energy |52) . 

2.3.8 Decay of vortices in quantum turbulence at low temperatures 

In this section, we consider the decay process of QT at low temperatures, at which the normal ffuid 
component is negligible and mutual friction does not occur. In this case, there is no dissipation in 
vortex lines at large scales and dissipation occurs only at scales comparable to the core radius ^. 
Energy at large scales cannot dissipate, but ffows to smaller scales via the Richardson cascade and 
the Kelvin- wave cascade |32) . 

Feynman first proposed a dissipation mechanism of QT at zero temperature |10J. In this 
mechanism, large vortex loops are broken up into smaller loops via the Richardson cascade, and 
the smallest vortex rings, with radii comparable to the atomic scale, decay into excitations, such as 
rotons. This mechanism is not currently accepted. Later, Vinen considered the decay of superfluid 
turbulence at finite temperatures in an examination of his experimental results from thermal 
counterflow, and proposed Vinen's equation |15) . At finite temperatures, vortices in turbulence are 
randomly spaced with no polarization, and there is only one length scale /. The energy and the 
decay of vortices spreads over a wide range of scales. The decay of the total energy is expressed in 
terms of a characteristic velocity Vs — n/inl and a characteristic time constant t = l/vg as, 



where x is a temperature dependent dimensionless parameter p3j . Incorporating the vortex line 
density L = this equation becomes 
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This is a type of Vinen's equation without a generation term proportional to {vs ~ v^jL'^^'^ due to 
thermal counterflow. The solution of this equation is 

T = ^+X^t, (2.56) 

where Lq = L{t = 0). At large t, L behaves as L oc t^^. The experimental observations of thermal 
counterflow can be explained by this solution, and x can be obtained as a function of temperature. 
Although the decay of QT at finite temperatures is due to mutual friction, this solution can describe 
QT at zero temperature without mutual friction. Using a vortex filament model, Tsubota et al. 
performed numerical simulations of QT without mutual friction, and estimated the value of x [2 2) . 
Equation (I2.56P can also be formulated by assuming that the motion of the quantized vortices 
dominates the fluid dynamics and their energy A/(47r)K^L is their primary contribution to the 
total energy [72]. By using Eq. (I2.50p . we obtain 

dE A o dL , 9 ■> , ^ 

— = —K^-r = -v'k^L^. 2.57 
dt An dt ^ ^ 

The solution of this differential equation is 

1 1 A 



L Lq Anv' 



r\ (2.58) 



which is the same result as Eq. (j2.56p . The effective kinematic viscosity v' in Eq. (|2.53p and that 
in Eq. (|2.58p are naturally different. We denote them as v'y^ for Eq. (|2.53p and v'y for Eq. (j2.58p 

m- 

Two important questions arise: What causes the difference in decay between Eqs. (j2.58l) {L cx 
t-^) and (1^35)1 (L cx t-^'"^), and what is the origin of the decay mechanism at zero temperature? 
The answer to the first question comes from the structural differences of vortex tangles in QT. 
When a vortex tangle supports the Kolmogorov energy spectrum, vortices form an inertial range 
in which the tangle is self-similar, creating polarized vortex bundle structures that are different 
from a completely random distribution of vortex lines in QT at finite temperatures. In this case, 
QT decays as i cx t^^/"^ . However, when a vortex tangle is dilute and random, with no correlation, 
there is only one length scale I in the tangle and QT decays as L cx even at zero temperature. 
This different behavior of L enables us to understand the structure of a vortex tangle in QT. The 
behavior of i cx has been observed in turbulence by injecting negative ions into ^He |72| and by 
a vibrating grid in "^He-B [73]. As discussed in the previous subsection, the behavior of L cx t~'^/^ 
has been observed in turbulence by an impulsive spin down of superfluid ^He, by injecting negative 
ions into superfluid "^He [71) . and by a vibrating grid in superfluid ^He-B [73j . Bradley et al. also 
observed the crossover from L cx t^^/'^ to L cx behavior of grid turbulence in '^He [73]. The 
structure of a generated vortex tangle depends on the velocity of the grid, and the decay of the 
turbulence becomes L cx t^^/'^ for a dense vortex tangle with a rapidly oscillating grid, and L cc t^^ 
for a dilute vortex tangle with a slowly oscillating grid. By using Eqs. (|2.53l) and (|2.58p . Walmsley 
and Golov estimated the effective kinematic viscosity and z/y for turbulence in ^He and showed 
that they are almost same at high temperatures T > 1.0 K for a normal fluid, and that becomes 
much smaller than at low temperatures with no normal fluid [72) . 

For the second question, there are several possible answers. The first is acoustic emissions at 
vortex reconnections, similar to that of eddies in a classical fluid. Numerical simulations of the 
GP model support acoustic emissions at every reconnection [^HJ [37] . The energy dissipation for 
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each reconnection is about 3^ times the vortex hne energy per unit length: ~ AK^/47r. In the 
case of superfluid ^He, however, ^ is quite small and the dissipation due to reconnections can be 
negligible. This dissipation mechanism may emerge for the turbulent state of atomic Bose-Einstein 
condensates [30] ■ By using the GP model, Kobayashi and Tsubota investigated the decay from 
fully developed turbulence and reported both L cx t^^/"^ to L (x behaviors [751. In their 
simulation, ^ is close to the system size, and acoustic emission upon reconnection is the main 
dissipation mechanism. Another possible answer is the radiation of phonons from high-frequency 
Kelvin waves 0^ . Vortex reconnections excite Kelvin waves whose wavelength is on the order of I 
[ini [331 [221 132 • Although Kelvin waves with wavelength ^ I cannot cause effective radiation, the 
Kelvin- wave cascade creates waves with much shorter wavelengths [331 SZl SHI SS [SOI [S] ■ 

In actual experiments, one must also consider vortex diffusion as an origin of the decrease of 
L. When a vortex tangle is inhomogeneous or local and the probe observing the turbulence is also 
local, a vortex tangle may escape from the observable region. Using the vortex filament model, 
Tsubota et al. numerically studied the diffusion of an inhomogeneous vortex tangle [7B|. The 
effects of diffusion can be quantitatively evaluated by the equation: 

= -X^i(x, tf + DV'L{^, t). (2.59) 

Here, L{yL,t) is the space-dependent line length density, and D is the diffusion constant. The 
numerical simulation indicated that D 0.1k 

2.3.9 Quantum turbulence created by vibrating structures 

Recently, vibrating structures, such as discs, spheres, grids, wires, and tuning forks, have been 
widely used for research into QT [77]. Despite detailed differences between the structures, the 
experiments have shown some surprisingly common phenomena. This trend started with the 
pioneering observation of QT on an oscillating microsphere by Jager et al. [78^ . Subsequently, 
many groups have experimentally investigated the transition to turbulence in superfluid ^He and 
3He-B by using grids [721 [501 IH El (HI [73] , wires [SI [53 [Ml ISZl [HS] , and tuning forks [HH [HD]- 
The details of these observations were described in a review article [77] . 

Here, we will describe briefly the essence of the observations by referring to a typical result 
from Yano et a^.|91|. A thin superconducting (typically NbTi) wire with a micron-size radius is 
formed into a semicircle, and its two edges are attached to the wall of a vessel. The wire vibrates 
in resonance with a Lorentz force due to an alternating electric current under a static magnetic 
field. Figure [2771 shows the response velocity of the wire as a function of the driving force. When 
the driving force is relatively low, the wire moves smoothly. However, if the driving force exceeds 
some value, the velocity suddenly drops, indicating that the wire does not vibrate so much despite 
an increase of the driving force. The energy injected from the drive must go somewhere, and the 
only possible escape is the excitation of quantum turbulence. The response of the wire clearly 
shows a hysteresis between the upward and downward sweeps of the driving force. 

Experimental studies conducted by many groups reported common behaviors, independent of 
the structural details, such as type, shape, and surface roughness. However, there were also some 
differences. The key points are the critical velocity for the onset of turbulence, the shift of the 
resonance frequency, and the hysteresis and the drag coefficient. 

The observed critical velocities were in the range from 1 mm/s to approximately 200 mm/s, 
much lower than the velocity necessary for intrinsic nucleation of quantized vortices, which is on 
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Figure 2.7: Typical response of a vibrating wire in superfluid ^Hc at 30 mK. Courtesy of H. Yano. 

the order of 10 m/s. It is known that there are initially some remnant vortices in superfluid ^He[92|. 
and the transition to turbulence should come from their extension or amplification. This scenario 
was confirmed by numerical simulations using the vortex filament model [94) and an experiment 
using a wire free from remnant vortices [88], which will be described later. However, it should be 
noted that this scenario is applicable only to ^He. The transition to turbulence observed in the 
B phase of superfluid '^He [53J [73] can be related with the intrinsic nucleation of vortices, which 
is known to occur near a solid surface at very small velocities, on the order 4 mm/s, depending 
on the surface roughness This issue was discussed in detail in a review article [77]. The 

dependence of the critical velocity Vc on the oscillation uj is very important. Most experiments 
have provided evidence for a universal scaling property Vc ~ y/noj, independent of the geometry 
of the vibrating structure. This relation is trivial from dimensional analysis, but it is necessary 
to investigate its physical origin. Hanninen and Schoepe discussed the scaling property [95 . This 
relation was obtained from the "superfluid Reynolds number" Re = v£/k with some characteristic 
length scale £ [53]. It is also possible to understand the scaling by extending the scenario of the 
dynamical behavior of a vortex tangle in counterflow turbulence at constant velocity |97j to the 
case of oscillating flows. However, more detailed investigation of the issue is required. 

The transition to turbulence is accompanied with a characteristic shift of the resonance fre- 
quency. In the wire experiments Yano et al. observed that the resonance frequency decreased with 
increasing drive force in the laminar regime, but increased in the turbulent regime |86] . Bradley et 
al. also observed an upward shift of the resonance frequency at the onset of turbulence [82j . Such 
an increase in the resonance frequency is surprising, because it indicates either an increased stiff- 
ness or a decreased effective mass of the wire. If vortices were amplifled and a tangle surrounds the 
wire, the effective mass should increase. The reasons for the frequency shift are not yet resolved, 
and more study is required. 

Many groups have reported some hysteresis in the transition to turbulence, which arises from 
the stability of the laminar and turbulent states. By using an oscillating microsphere, Jager et al. 
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observed that the transition from laminar to turbulent response was accompanied by significant 
hysteresis at low temperatures [78 . Below 0.5 K, turbulence was observed to be unstable and the 
flow switched intermittently between turbulent and laminar flow. Schoepe analyzed the lifetime 
of turbulence in connection with the statistical fluctuations of the vortex line length density [S5] . 
Yano et al. [M] and Bradley et al. |:82 also observed hysteresis when using vibrating wires. 

The net drag force is another important quantity when we consider the transition to QT. The 
drag force F is often expressed in terms of a drag coefficient Cd, defined by the equation [99] 

F=^Cj,pAU^ (2.60) 

where p is the density of the fluid, A is the projected area of the object normal to the flow, and U 
is the flow velocity. The behavior of Cd is well known in a steady classical flow past an obstacle. In 
the case of laminar potential fiow of an ideal fiuid past a cylinder, the flow is symmetric about the 
plane through the center of the cylinder and normal to the flow. As a result, the net force vanishes 
with Cd = (the d'Alembert paradox). For laminar viscous fiow, the drag force is approximately 
proportional to U so that Cd ~ (the Stokes' formula). In the case of fiow at a high Reynolds 
number, the fiow behind the obstacle accompanies a wake, which destroys the symmetry and leads 
to the d'Alembert paradox, and Cd becomes of order unity. The value of Cd depends on the 
geometry of the obstacle [95]; Cd is approximately unity for a disc or cylinder, and about 0.3 for 
a sphere. Also, in quantum turbulence created by vibrating structures, it is possible to obtain 
the drag coefficient from the dependence of velocity on the driving force, as in Fig. 12.71 Much of 
the data on microspheres [75] , mm-scale spheres [100] , grids [TS] IHO] , and quartz forks [90 show 
the classical analogue behavior in which the drag coefficient is about U^^ in the laminar regime 
and of order unity in the turbulent regime. This is another classical analogue of QT. Although 
the experimental results are trivial, we need a clearer understanding of the phenomena. The drag 
coefficient characteristic of the turbulent regime was numerically confirmed by simulation of the 
vortex filament model jlOl) . 

The remnant vortices play several important roles in these phenomena, which are discussed in 
the remainder of this subsection. 

Hanninen et al. performed a numerical simulation using the vortex filament model, and de- 
scribed how remnant vortices develop into a tangle under an AC superfiow, as shown in Fig. 12.81 
[94] . A smooth solid sphere of radius 100 fim was placed in a cylindrical vessel filled with super- 
fiuid ^He. Generally, we do not know how remnant vortices remain in a given geometry, but here 
they were assumed to initially extend between the sphere and the vessel wall. When an oscillating 
superfiow of 150 mm^^ at 200 Hz was applied, Kelvin- waves resonant with the fiow were gradually 
excited along the remnant vortices. The amplitude of the Kelvin- waves grew large enough to lead 
to self-reconnection and the emission of small vortex rings. These vortices gathered around the 
stagnation points, repeatedly reconnected, and were amplified by the fiow, eventually developing 
into a vortex tangle. This simulation likely captured the essence of what occurred in the experi- 
ments. However, the simulation was not necessarily satisfactory, and some problems remain. The 
first is the critical velocity. The situation of Fig. 12.81 is quite similar to the experiment using 
a microsphere [78] and gives a critical velocity of about 120 mm~^, which is much larger than 
the observed value of about 40 mm^^. The second problem is hysteresis. The simulation does 
not produce any hysteresis between upward and downward sweeps. If the velocity is taken above 
the critical value for the formation of a tangle, and then reduced below this value, the vortices 
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immediately decay, eliminating the turbulent state. These two difficulties may arise from the arti- 
ficially smooth surface of the sphere. A solid surface is generally rough, and even tiny bumps can 
act as pinning sites for quantized vortices with such a thin core [191 1102j . If such pinning effects 
on the surface are taken into account, they should induce additional disturbances of vortices and 
help the transition to turbulence, thereby reducing the critical velocity. Regarding the hysteresis, 
some vortices may be trapped by pinning cites through the downward sweep, which can make the 
physics different from the upward sweep. However, we do not know how to consider the effect of 
surface roughness, namely multi-pinning cites for numerical vortex dynamics, and the development 
of applicable methods is urgently needed. 



(a) (b) (c) 




Figure 2.8: Evolution of the vortex line near a sphere of radius 100 /im in an oscillating superflow of 
150 mm^i at 200 Hz [Hanninen, Tsubota and Vinen: Phys. Rev. B 75 (2007) 064502, reproduced 
with permission. Copyright 2007 by the American Physical Society]. 

The roles of remnant vortices were revealed by several interesting experiments by the Osaka 
City University group. First, Hashimoto et al. succeeded in setting a vibrating wire free from 
remnant vortices using a unique idea|87|. They added a small chamber with a pinhole in a sample 
cell, and filled the chamber extremely slowly with superfluid ''He. The vortices were likely filtered 
out by the pinhole, providing a vibrating wire free from remnant vortices. The wire never caused 
a transition to turbulence, even when the velocity exceeded 1 m/s. Next, Goto et al. set two wires 
in the small chamber. By using the same procedure, they obtained one wire (labeled A) free from 
remnant vortices and another (labeled B) with remnant vortices [55] • Wire B can create turbulence 
by itself from remnant vortices and the resulting tangle emits vortex rings. Although wire A alone 
never caused turbulence, it could create turbulence if it received seed vortices from wire B. This 
observation clearly demonstrated that the turbulence arises from remnant vortices. This behavior 
was confirmed by numerical simulation of the vortex filament model [851 1101) . 

2.3.10 Visualization of quantized vortices and turbulence 

There has been little direct experimental information regarding the flow in superfluid ^He. This 
is mainly because usual flow visualization techniques are not applicable to cryogenic superfluids. 
However, this situation is rapidly changing |103j . For QT, one can seed the fluid with tracer 
particles in order to visualize the flow field and possibly quantized vortices, which are observable 
by appropriate optical techniques. 

A significant contribution was made by Zhang and Van Sciver |104j . Using a particle image 
velocimetry (PIV) technique with 1.7-/im-diameter polymer particles, they visualized a large-scale 
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turbulent flow both in front of and behind a cylinder in a counterflow in superfluid "^He at finite 
temperatures. In classical fiuids, such large-scale turbulent structures are seen downstream of 
objects such as cylinders, with the structures periodically detaching to form a vortex street. In the 
present case of "^He counterflow, the locations of the large-scale turbulent structures were relatively 
stable, and they did not detach and move downstream, although local fluctuations in the turbulence 
were evident. 

Another significant contribution was the visualization of quantized vortices by Bewley et al. 
|105j . In their experiments, the liquid helium was seeded with solid hydrogen particles smaller 
than 2.7 /itm at a temperature slightly above T\, after which the fiuid was cooled to below T\. 
When the temperature was above Tx, the particles were seen to form a homogeneous cloud that 
dispersed throughout the fiuid. However, on passing through T\, the particles coalesced into web- 
like structures, as shown in Fig. 12.91 Bewley et al. suggested that these structures represent 
decorated quantized vortex lines. They reported that the vortex lines appear to form connections 
rather than remaining separated, and were homogeneously distributed throughout the fluid. The 
fork-like structures may indicate that several vortices were attached to the same particle, as in- 
dicated by numerical simulations of vortex pinning [1(16 . Using the same technique, Paoletti et 
al. obtained the trajectories of tracer particles visualized thermal counterflow |107j . The observed 
trajectories showed two distinct types of behavior. One group consisted of trajectories that moved 
in the direction of the heat flux, while the other consisted of those that opposed this motion. The 
former trajectories were smooth and uniform, but the latter could be quite erratic. Particles of 
the former trajectories were probably dragged by the normal fiuid, while those of the latter were 
trapped in vortex tangles. 

Here, it is necessary to understand whether such tracer particles follow the normal fiow or 
the superflow, or an even more complex flow. Poole et al. studied this problem theoretically 
and numerically, and showed that the situation changes depending on the size and mass of the 
tracer particles [TUS]. However, the situation is so complicated that many issues remain unanswered 

[inniiiiniiiiiiiiia. 
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Figure 2.9: Visualization of quantized vortices by intensity-inverted images. Each image shows the 
full field of view of 6.78 x 6.78 mm^. Multiple particles are often trapped in each visible vortex, 
and trapped particles tend to be uniformly spaced along the vortex core [1071 . Courtesy of E. 
Fonda, K. Gaff, M. S. Paoletti, and D. P. Lathrop. 
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Chapter 3 



Quantized vortices in atomic 
Bose-Einstein condensates 

3.1 Introduction 

The achievement of Bose-Einstein condensation in trapped atomic gases at uhra-low temperatures 
|113l 11141 1115) has stimulated intense experimental and theoretical activity in modern physics. In 
this section, we review the physics of quantized vortices in atomic-gas BECs. While quantized 
vortices have been thoroughly studied in superfluid helium |116) . there has been a resurgence of 
interest in vortices in atomic BECs because of the peculiar features of this system: 

1. Small gas parameter 

Because the gas is dilute, the GP equation (|2.5p gives a quantitatively accurate description 
of the static and dynamic properties of the atomic condensate [501 1117] . Therefore, the vor- 
tex structure and dynamics can be discussed by a more fundamental approach than with 
superfluid helium, in which the interaction cannot be described in such a simple local form. 
Furthermore, the diluteness of a gas leads to a relatively large healing length that character- 
izes the vortex core size, thus enabling the direct experimental visualization of vortex cores 

2. Trapping potential 

The finite size effect of the trapping potential and the associated density inhomogeneity 
yield new characteristics of vortex physics. This also provides a new spectroscopy to study 
characteristic vortex dynamics through collective motion of the inhomogeneous condensate 

[US [ml [na [126]. 

3. Laser manipulation and rotation 

The manipulation of a condensate wave function via external (magnetic or optical) fields pro- 
vides a versatile scheme to control the vortex states. Laser beams can yield effective repulsive 
or attractive forces on atoms, creating confining potentials, periodic lattice potentials, and 
impurity /obstacle potentials for atomic condensates. Vortex formation can be achieved by 
the mechanical rotation of a laser-created optical spoon [1191 11201 1121] or the transfer of 
orbital angular momentum from a Laguerre-Gaussian laser beam [127] . The tunability of 
the rotation over a wide range is useful to study various vortex phases and their dynamics 
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in trapped BECs. This provides an opportunity to study a rapidly rotating condensate with 
a dense vortex lattice [1201 11281 1129] , which is not easy to accomplish in superfluid helium 
systems. Such rapidly rotating BECs present a so-called mean-field quantum-Hall regime 
described by the orbitals of the lowest Landau levels (LLL), and also a regime dominated by 
strong many-body correlation |130[ 11311 1132| . Moreover, laser- induced coherent transition 
of atoms between different hyperfine levels can control the phase profile in such a way that 
the condensate has a circulation [118) . Recently, the laser field has played the role of an 
artificial vector potential |133l I134j , which could be an alternative means of manipulating 
the rotational properties of condensates. 

4. Detection 

The vortex cores can be directly visualized through the observed density profile by a time- 
of-flight (TOE) technique [TTl [111 [HOI [M] ■ The TOE technique involves switching off the 
trapping field (magnetic or optical) at time t — and taking an image of the BEC several 
(typically 5 to 25) milliseconds later. Switching off the trap allows the atomic gas to expand 
until a laser beam probe becomes available to observe the density profile. Images are taken 
by shining a resonant laser beam into the atomic gas and using a CCD camera to observe 
the shadow cast by the absorption of photons, from which the integrated atomic density can 
be determined. 

The matter- wave interference technique is also applicable to vortex detection [1351 1136] . The 
presence of vortices can be revealed by finding the dislocations in the interference fringes 
between expanding condensates |137j . In addition, the Bragg scattering of matter waves from 
laser beams can be used to detect the vortex state, because the behavior of the scattered 
condensate is sensitive to the spatial phase distribution of the initial state [1381 1139) . 

5. Multicomponent condensates 

Since it is possible to load and cool atoms in more than one hyperfine spin state or more 
than one atomic element in the same trap, multicomponent BECs with internal degrees of 
freedom can be created experimentally HHl IMl 1113 [HI [113 [Hi [Wl [HI [US] ■ Mul- 
ticomponent condensates allow the formation of various unconventional topological defects 
with complex properties that arise from interactions between different components, provid- 
ing a new platform for the study of unconventional vortices. This content relates closely to 
other fields of physics, such as superfluids ^He |150) . unconventional superconductors |151| . 
quantum Hall systems |152j . as well as high-energy physics and cosmology jl53) . Atomic 
BEC is an extremely flexible system for the study of such topological defects, since optical 
techniques allow the control of the condensate wave functions. An external field can couple 
the internal sublevels of the condensate and cause a coherent transition of the population. 
This coherent transition can be used to control the spatial variation of the condensate wave 
functions, resulting in an "imprinting" of a phase pattern onto the condensate. In most 
schemes, the spatial configuration of the field, the intensity and detuning of the laser fields, 
and the phase relationship between the different fields must be carefully controlled to create 
the correct phase pattern of topological excitations through the complex internal dynamics. 

6. Feshbach resonance 

The most salient feature of the cold atom system is that a field-induced Feshbach resonance 
can tune the s-wave scattering length between atoms [154] . which determines the strength 
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of atom-atom interaction. A Feshbach resonance occurs when a quasi-bound molecular state 
in a closed channel has an energy equal to that of two colliding atoms in an open channel. 
This technique can change the scattering length over a wide range, from negative values to 
positive ones, thus creating condensates with strongly attractive interactions [155) . strongly 
repulsive interactions |156| . or dominant long-range dipole-dipole interactions |157| . It also 
provides new physics associated with the formation of quasi-bound molecules, provided by 
the coherent transition between an atomic condensate and a molecular one |158| . In the 
case of fermions, this has been used to induce pairing of two fermions, which results in the 
formation of fermion condensates by controlling the character of the two-body interactions 
from strongly-bound bosonic molecules to weakly- interacting fermion (Cooper) pairs jl59| . 
Vortices in these exotic condensates are expected to show quite intriguing properties. A 
review of Feshbach resonance was recently given by Chin et al., fl60] . 

There are several excellent reviews of quantized vortices in ultracold atomic BECs. The basics 
of the theoretical formulation and the early research on the physics of quantized vortices were re- 
viewed by Fetter and Svidzinsky |161) . Recent progress of the extensive study of quantized vortices 
is described in Ref. |162l 1163] . These reviews were mainly concerned with the properties of a con- 
ventional single-component BEC, whose order parameter is scalar and interaction is characterized 
only by an s-wave scattering length. In this chapter, therefore, we will also discuss vortex physics 
in unconventional BECs, which was not covered by these reviews. We begin with the basic theory 
and experiments on quantized vortices in conventional BECs. Then, we describe vortices in some 
unconventional condensates realized in ultra-cold atomic systems, such as spinor condensates and 
dipolar condensates. 



3.2 Vortices in single- component Bose-Einstein condensates 

First, we review vortices in a single-component BEC in a trapping potential from both theoretical 
and experimental points of view. A more detailed review and references can be found in Ref. 
|1621 1163j . The contents of this section will give the basic information necessary for subsequent 
sections. After reviewing the theoretical formulation for trapped BECs, we describe the properties 
of a single vortex and a vortex lattice in trapped BECs, together with various interesting current 
topics. 



3.2.1 Theoretical description of ultra-cold atomic-gas BECs 
Gross-Pitaevskii equation 

The many-body Hamiltonian for bosons in a trapping potential T4x(r) is given by 



H = / drV/(r) 



2m 



'ip{Y) + - dv dv'i^\r)'4}\Y')V{v~T')i,{v')i,{r) (3.1) 



with the bosonic field operator -0(r) and the two-body interaction V{y — r'). The dynamics of ip 
are governed by the Heisenberg equation ihdip/dt — [ij}, H]. In the low temperature limit, the field 
operator can be approximated by a classical field known as "condensate wave function" ^ = (?/;), 
and the trapped BECs can thus be described by the CP equations: 



ih 



a*(r,i) 
'di 



2m 



+ l/cx + .9|4'(r,t)|' 



*(r,t). 
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The trapping potential typically has the form of a harmonic oscillator Vex = {''n/2){u!'^x'^ +LUyy'^ + 
ijj'j.z^). Comparison of the kinetic and trap energies introduces a characteristic length scale Oho = 
yjhjrnu) with oj = (w^Wj/Wz)^/'^. The time-dependent GP equation p.2p can be used to explore the 
dynamic behavior of the condensate, which is characterized by variations of the order parameter 
over distances larger than the mean distance between atoms. This equation is valid when the s- 
wave scattering length is much smaller than the average distance between atoms, and the number 
of atoms in the condensate is much larger than unity. The success of the quantitative description 
of the trapped BECs is described in the books by Pethick and Smith [30], and Pitaevskii and 
Stringari pTT] . 

To determine the ground state of a trapped BEG, we can write the condensate wave function 
as 5'(r,i) — $(r)e"*''*/'', where $(r) obeys the time-independent GP equation 

$(r) = ^^(r). (3.3) 

Here, $ is normalized to the number of condensed particles /c?r|<i>(r)p — Nq, which determines 
the chemical potential fj,. In the dilute limit (the gas parameter n|ap is typically less than 10""^, 
where n is the average density of the gas) the condensed particles A'o can be approximated as the 
total number N, because depletion of the condensate is small as N' — N — Nq oc y/ n\a\^N <C N. 
The time-independent GP equation p.3|) is also derived by minimizing the GP energy functional 

£;[$, $*] = JdVi>* (^-^^ + + $ = Skin + Et, + Sint, (3.4) 

subject to the constraint of a fixed particle number N. This constraint is taken into account by 
the Lagrange multiplier /i as S{E — fiN)/6^* = 0, where fj, is the chemical potential that ensures 
a fixed N. 

Although the exact solutions of the ground state can be generally obtained only by solving Eq. 
(j3.3l) numerically, an approximate analytic solution can be found when the interaction energy Ei,^t 
is much larger than E'kin |164j . To see this, let us neglect the anisotropy of the harmonic potential 
and assume that the cloud occupies a region of radius ^ R, so that n ^ N/R^. Thus, the scale 
of the harmonic oscillator energy per particle is ^ mui^R^ /2, while each particle experiences an 
interaction with the other particles of energy ~ gN/R^. By equating these energies, the radius 
is found to be i? ^ aho(87rAfa/aho)^^^- Since the kinetic energy is of order fP/2mR^, the ratio of 
the kinetic to interaction (or trap) energies is ^ (iVa/aho) "^^^^ In the limit iVa/oho 3> 1, which is 
relevant to most of the experiments, the repulsive interactions significantly expand the condensate, 
so that the kinetic energy associated with the density variation becomes negligible compared to 
the trap and interaction energies. As a result, the kinetic-energy operator can be omitted in Eq. 
(|3.3p . which gives the Thomas-Fermi (TP) parabolic profile for the ground-state density 

n(r) ^ |$TF(r)P = ^"^''"^''^ 9 [m - V;,(r)] , (3.5) 

where Q{x) denotes the step function. The resulting ellipsoidal density in three-dimensional (3D) 
space is characterized by two types of parameters: the central density no — i-t/g and the three 
condensate radii i?| = 2^/mu]j [j = x,y,z). The chemical potential fi is determined by the 
normalization / drn(r) = TV as /i = {nuj/2){15Na/ai,oy/^. In the TF regime, the time-dependent 
GP equation (|3.2p reduces to the incompressible hydrodynamic equation for superfluids |117j . 



-^ + K. + ff|<i>(r)P 
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The Bogoliubov-de Genne equation 



The spectrum of elementary excitations of the condensate is an essential ingredient in calculations 
of thermodynamic properties. To study the low-lying collective-excitation spectrum of trapped 
BECs, the Bogoliubov-de Gennes (BdG) equation coupled with the GP equation is a general 
formalism. Let us consider the equation of motion for a small perturbation around the stationary 
state $, which is a solution of Eq. p. 31) . The wave function may take the form ^(r, t) ~ ['l'(r) -I- 
itj(r)e^'"5* — u* (r)e*'^J*]e^*'^*. Inserting this ansatz into Eq. p.2p . and retaining terms up to the 
first order in u and u, we obtain the BdG equation: 



' = f^A ? :( ), (3.6) 



where L{r) = ~h'^V^/2m + V^^ - fi + 2g|$(r)p, and LUj are the eigenfrequencies related to the 
quasiparticle normal- mode functions Uj(r) and Vj{r). The mode functions are subject to the 
orthogonality and symmetry relations / dr[uiU* — ViV*] = 6ij and / dr[uiV* — ViU*] = 0. 

Since the energy fkoj of this quasiparticle is defined with respect to the condensate energy 
(Eq. p.4p with the stationary solution $), the presence of quasiparticles with negative frequencies 
implies an "thermodynamic" instability for the solution $; if there is energy dissipation in the 
system, the excitation of negative-energy modes lowers the total energy, and the modes grow 
spontaneously to relax $ into a more stable state. In other words, there is a path in configuration 
space along which the energy decreases. This argument is closely related to the Landau criterion 
for superfluidity (Landau instability). On the other hand, one can note that, since the matrix 
element of Eq. (13. 6p is non-Hermitian, the eigenfrequencies can assume complex values. Therefore, 
the small-amplitude fluctuations of the corresponding eigenmodes grow exponentially, even for the 
energy-conserving time development. This is known as "dynamical instability". These instabilities 
play a key role in vortex dynamics. 

3.2.2 Vortex experiments in rotating BECs 

The first experimental detection of a vortex in an atomic BEG was made by Matthews et al. in 
1999 |118) . Their study involved condensates of *^Rb atoms residing in two hyperfine states. We 
describe this experiment in the next section, and summarize here the experiments on vortices in 
single-component BEGs under rotation. 



Figure 3.1: Typical density profiles of a rotating condensate in the plane perpendicular to the 
rotation axis, taken by TOF measurement. From left to right, the rotation frequency incrases. 
Gourtesy of J. Dalibard. 

An intuitive method of creating vortices is to stir a condensate mechanically with a rotating 
"bucket". Using this scheme, Madison et al. observed the formation of vortices in a single- 
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component elongated BEC |119) [see Fig. 13. 1| . The condensate was trapped in a static axisym- 
metric magnetic trap which was deformed by a nonaxisymmetric attractive dipole potential created 
with stirring laser beams. This combined potential produced a cigar-shaped harmonic trap with 
a slightly anisotropic transverse profile. By rotating the orientation of the transverse anisotropy 
at a frequency fl, they observed the dynamic formation of a vortex above a certain critical value 
of Q. When D, was increased further, multiple vortices appeared, forming a triangular lattice. 
The quantized vortices could be directly visualized as "dips" in the transverse density profile of 
the TOF image. Following the experiments of the ENS group, other groups have also observed 
quantized vortices using slightly different methods based on the concept of the rotating bucket. 
Abo-Shaeer et al. observed a vortex lattice consisting of up to 100 vortices in a ^"^Na BEC [120) . 
Hodby et al. created a vortex lattice by rotating the anisotropic magnetic trap directly without 
using a laser beam, which is similar to the rotating bucket experiment |122| . Rotating an optical 
spoon made by multiple-spot laser beams or narrow focusing beams has been also used to nucleate 
vortices [165) . 

Haljan et al. created a vortex state by cooling an initially rotating thermal gas in a static 
confining potential [128) . Thermal gas above the transition temperature was rotated by a slightly 
anisotropic trap. After recovering the anisotropy of the potential, the rotating thermal gas was 
evaporatively cooled until most of the atoms were condensed. Although the number of atoms 
decreased through the evaporative cooling, the condensate continued to rotate because the angular 
momentum per atom did not change, and thus the vortex lattice was created. Since atoms can be 
selectively removed during evaporation, spinning up of the condensate can be efficiently achieved 
by removing atoms extending in the axial direction, as opposed to the radial direction, allowing a 
BEC with a high rotation rate to be obtained. Various properties of rapidly rotating BECs were 
studied, e.g., equilibrium properties [166) . collective dynamics of a vortex lattice |167l 1168] . vortex 
aggregation |169) . and the lowest Landau level regime [170] . 

3.2.3 Structure of a single vortex 

As a simple example, let us consider the structure of a single vortex in a condensate trapped by 
an axisymmetric harmonic potential Vcv.{t,z) = muj'^(r'^ + X'^z'^)/2 with the transverse trapping 
frequency uj± and the aspect ratio A — uJz/ujj_. The condensate wave function with a straight 
vortex line located along the z-axis takes the form $(r) = (j){r, z)e^'^^ with a winding number q 
and cylindrical coordinate {r,9,z). The velocity field around the vortex line is Vg = {qh/mr)9. 
Equation p.3p for a real function becomes 

-_f^ fd^ + + ^\ + f]^ + (r z) + .gJ ci,~^l(j, (3 7) 

2m \dr-^ r dr dz"^ J 2mr'^ 

Here, the centrifugal term q^h^ /2mr'^ arises from the azimuthal motion of the condensate. Equation 
p. 71) . which must be solved numerically, gives the structure of the vortex state. Asymptotically, 
the solution takes n{r, z) ~ no{r/£_y'^^ for r <^ ^, and n{r, z) = no(l — r"^ / R\ — z^ jB?^ — q^S^ jr^') 
for ^ ^ r < i?_L. The latter can be obtained in the TF limit Na/a^^o > 1 for Eq. (j3.7p . where 
no = /i/g is the density at the center of the vortex- free TF profile. We have defined the TF radii 
R\ = 2^/mw\, i?2 = '^n/mojl and the healing length ^ — {fP /2mgno)^/^ . The condensate density 
vanishes at the center, out to a distance of order ^, whereas the density in the outer region has the 
form of an upward-oriented parabola. Hence, the healing length ^ characterizes the vortex core 
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size; for typical BEC parameters, ^ ~ 0.2 /im. In the TF limit, the core size is very small because 
£,/Rj_ = hujj_/2fi ~ (15iVa/aho)~^^^ ^ 1- Increasing the winding number q widens the core radius 
due to centrifugal effects. 

This axisymmetric solutions holds for a BEC in a pancake trap with A 3> 1, but for a cigar- 
shaped trap with A <C 1 the axisymmetry is spontaneously broken as a result of vortex bending. 
Numerical simulation of the 3D GP equation revealed that the vortex bending is held stationary 
for the ground state of a rotating cigar-shaped condensate [1711 11721 I173j . This bending is a 
symmetry breaking effect, which happens even in a completely axisymmetric setup. Evidence of 
vortex bending in the ground state was observed by the ENS group [174j . 

3.2.4 Vortex stability 

Stability of the vortex state should be ensured to demonstrate that a BEC exhibits the superfluidity 
characterized by a "persistent current". At zero temperature, the BdG analysis shows that the 
excitation spectrum of a condensate with a centered vortex in a cylindrical system has at least one 
negative energy mode with a positive norm J dr(|up — \v\'^) > 0, localized at the core (so-called 
"anomalous modes" ) [175j . This implies that the single- vortex state is therniodynamically unstable 



|176j . In fact, in the absence of rotation, the centered vortex would decay by spiraling outward 
due to dissipation by the thermal atoms. 

As demonstrated in the experiments, imposing rotation on the system is a direct way to achieve 
vortex stabilization. If the system is under rotation, it is convenient to consider the corresponding 
rotating frame; for a rotation frequency fl = ilz, the integrand of the GP energy functional p.4p 
acquires an additional term: 



E' = I dr^* [~^ + V,. + ^\^\'-nL,]^, (3. 



where Lz = —ihixdy — ydx) is the angular momentum operator along the z-axis. The corresponding 
GP equation becomes 



-^^ + 14x + <?l*(r,i)|'-r!i. 
Zm 



^{v,t). (3.9) 



The energy associated with a single vortex line is predominantly contributed by the kinetic 
energy of the superfluid flow in the vortex, which can be estimated as 

El = / -mnvldr ~ —2Rz / vhnrdr = q^R^ In -F^ , 3.10 

J2 2 m \ S, J 

where n is the mean uniform density. Since Ei (x q^, the energy cost to create one q — 2 vortex is 
higher than that to create two q = 1 vortices, thus vortices with q > 1 are energetically unfavorable. 
Therefore, a stable quantized vortex usually has q = 1, and we will mainly concentrate on the q = 1 
vortex in the following discussion. 

If there is a quantized vortex along the trap axis, (L^) = Nh, the corresponding energy of the 
system in the rotating frame is E[ = Ei — Nhft. The difference between E[ and the vortex- free 
energy E'q determines the energetic favorability of a vortex entering the condensate. Since E'q 
is equal to the energy Eq in the laboratory frame, the difference is given by AE' = E[ — E'q — 
El — Eq — Nhfl. Thus, the critical rotation frequency 57c for the energetic stability of a vortex line 
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is given by flc = {Ei — EQ)/Nh. When fl exceeds flc, the single vortex state is thermodynamically 
stable. 

To calculate Ei more precisely, it is necessary to take into account the inhomogeneous effects 
of condensate density [177] . In the TF limit, for a condensate in a cylindrical trap = (an 
effective 2D condensate), the critical frequency is given by = {2h/mR']_)l-n.{0.888R±/£_). For 
an axisymmetric trap Vcx{r,z), the critical frequency is flc = (5fi./2TOi?^) ln(0.671i?j^/^). For a 
nonaxisymmetric trap, flc is slightly modified by a small numerical factor [1781 1179] . 

3.2.5 Vortex nucleation 

The critical rotation frequency Vic only indicates the energetic stability of the central vortex state. 
Vortex nucleation in a non-rotating condensate occurs when the trap is rotated at a higher fre- 
quency than flc to overcome the energy barrier that stops the transition from the nonvortex state 
to the vortex state [180) . The threshold of the rotation frequency for instability, leading to vortex 
nucleation, is related to the excitation of surface modes of the condensate |181[ I182[ I183[ 11841 1185) . 
According to the Landau criterion for rotationally invariant systems, the critical rotation frequency 
is given by fl^ = min(a;i//), where uji is the frequency of a surface mode with the profile ~ e*'^ 
and the quantum number I of the azimuthal {6-) direction. Above Vly, some surface modes with 
negative energy appear in the spectrum of a non-rotating condensate |180l 1182) , and their growth 
may trigger vortex nucleation. The negative-energy modes can grow exponentially in the presence 
of dissipation, caused, for example, by interactions with thermal atoms )184) . 

The frequency can explain vortex nucleation by a rotating thermal gas )128II183] . but cannot 
yet explain the results of external stirring potentials |119l I120[ 11211 1122) . For example, in the case 
of the ENS group )123) . the number of nucleated vortices had a peak near fl = 0.7uj±. Their 
experiments demonstrated that instability occurs when a particular surface mode is resonantly 
excited by a deformed rotating potential. The rotating potential of the ENS group mainly excited 
the surface mode with I = 2 (quadrupole mode). In a rotating frame with frequency fl, the 
frequency of the surface mode is increased by —Ifl. The resonance thereby occurs close to the 
rotation frequency VL = uJi/l. In the TF limit, the dispersion relation for the surface mode is given by 
LOi — \/lu}j_ with the trapping frequency u}j_ [186) . Hence, it is expected that the quadrupole mode 
with Z = 2 is resonantly excited at f2 = uj±/^/2 ~ 0.707cjj^. A theoretical study revealed that when 
the quadrupole mode is resonantly excited, an imaginary component in the excitation frequency 
appears in surface modes with high multipolarity [187) . This indicates that dynamic instability 
can trigger vortex nucleation even at zero temperature. This picture was further supported by an 
experiment of the MIT group [165) . in which surface modes with higher multipolarity (/ — 3, 4) were 
resonantly excited using multiple laser-beam spots. The largest number of vortices were generated 
at frequencies close to the expected values fl = uj±/Vl- However, a recent theory proposed an 
alternative mechanism, noting that the single vortex state can be regarded as the first Zeeman- 
like excited state, and its resonance frequency is quantitatively consistent with all of the above 
experiments [188] . 

The observation of vortex nucleation and lattice formation [121) has been well reproduced by 
numerical simulations of the time-dependent GP equation ([311) pMl [MHIWl HMlirggl [TMIIT^ . 
The simulation results clarified the following dynamics. When a sufficiently high rotation is applied 
to an initially non-rotating condensate, it becomes elliptical and undergoes a quadrupole oscillation 
[Fig. I3.2f a)]. Then, the boundary surface of the condensate becomes unstable, and generates 
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Figure 3.2: Dynamics of vortex lattice formation in a rotating BEC. The graphs show both the 
condensate density (upper panels) and phase (lower panels), being at (a) t = 300 ms, (b) 370 ms, 
(c) 385 ms, (d) 410 ms, and (e) 550 ms after the start of rotation. The phase varies continuously 
from (white) to 2 tt (black). The simulations were performed using the 2D GP equation with 
the experimentally accessible parameters; see Ref. [189) for details. 

ripples that propagate along the surface [Fig. I3.2f b)]. It is also possible to identify quantized 
vortices in the phase profile. As soon as the rotation starts, many vortices appear in the low- 
density region outside of the condensate [Fig. I3.2f a)]. Since quantized vortices are excitations, 
their nucleation increases the energy of the system. Because of the low density in the outskirts of 
the condensate, however, their nucleation contributes little to the energy or angular momentum. 
Since these vortices outside of the condensate are not observed in the density profile, they are 
called "ghost vortices". Their movement toward the Thomas- Fermi surface excites ripples [Fig. 
I3.2r c)]. It is not easy for these ghost vortices to enter the condensate, because that would increase 
both the energy and angular momentum. Only some vortices enter the condensate cloud to become 
"real vortices" with the usual density profile of quantized vortices [Fig. I3.2f d)]. eventually forming 
a vortex lattice [Fig. I3.2f e) and (f)]. 

3.2.6 Dynamics of a single vortex 
Precession of an ofF-centered vortex 

Precession of a vortex core upon the off-center of the condensate is a simple example of vortex 
motion. Core precession can be described in terms of a Magnus force effect. A net force on a 
quantized vortex core creates a pressure imbalance, resulting in core motion perpendicular to both 
the force and the vortex quantization axis. In the case of trapped BECs, these net forces can be 
caused by either condensate density gradients |1781 11961 1197) or drag due to thermal atoms |198| . 
The former may effectively act as a buoyancy of the vortex. Typically, the total buoyant force is 
towards the condensate surface, and the net effect is a precession of the core around the condensate 
center via the Magnus effect. The latter causes radial drag and spiraling of the core towards the 
condensate surface due to energy dissipation. 

Core precession was experimentally investigated by Anderson et al. [199] . Initially, a vortex was 
prepared to be off-centered, and its precession frequency was determined from the vortex position 
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subtracted directly from the snapshots of the density profile. The vortex core precessed in the same 
direction as the vortex fluid flow around the core. These features can be understood as a result 
of the anomalous mode (with negative energy and positive norm); when the energy is negative, 
the precession has the same direction as the vortex flow. The observed excitation frequency also 
agrees with the theory [178] and the numerical simulations [1961 11971 1200) . 

Vortex wave 

As in the studies of superfluid helium, vortex waves are also an interesting subject in trapped 
BECs. The dispersion relation of the vortex wave can be obtained by BdG analysis of a single 
vortex state, as shown in Fig. l3.3f aV In this calculation, the perturbation can be taken as 
S-^ = e*^[ufe^,ie*(''"^+'^^"*^ - vl^ ^g-j(fc^z+/e-(^t)] the cylindrical system, where and I refer to 
the wave number and the angular quantum number along the z-axis. The lowest modes along the 
radial direction can be classified into three groups: the Kelvin wave {I — —1) [Fig. I3.3f b)]. the 
varicose wave {I = 0) [Fig. I3.3f c)]. and the surface waves {I ^ 0,-1). Because of the finite-size of 
the system, the dispersion of the Kelvin wave is well-described by |201j 



where Tc is on the order of the core size, fco ^ 2n/ the smallest wave number, and ujq < Q the 
frequency of the anomalous mode described above. The varicose wave with / = is known in 
classical fluids as the axisymmetrically propagating mode, in which the core diameter of a vortex 
oscillates along the vortex line. 



Figure 3.3: (a) Dispersion relation of the axisymmetric single-vortex state of Kelvin waves with 
I ~ —\ (dashed line), varicose waves with I — (dotted line), and several surface waves with 
I > (solid lines), calculated for a cylindrically symmetric condensate with a vortex (see Ref. 
|204| V (b) and (c) show schematic illustrations of the Kelvin and varicose waves of the vortex line, 
respectively. 

The collective modes of the condensate have a strong influence on the excitation of vortex waves 
due to the finite size effect. Bretin et al. observed that when superimposed I = ±2 quadrupole 
modes were excited with equal amplitudes in the single vortex state, the oscillation of the I — —2 
mode decayed faster than that of the I = +2 mode |125j . This indicates that the I — —2 mode 
decays to Kelvin modes through a non-linear Beliaev process, which is supported by theoretical 
analysis based on the BdG equation |202] and by numerical simulations [201) . 




(3.11) 
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When an additional velocity V is applied along the z-axis, the dispersion of the vortex waves 
behaves as — Vk^, and the frequency can become negative above the critical velocity. This 
is known as the Donnely-Glaberson instability in superfluid helium |203] . and it results in the 
amplification of Kelvin waves |204j . This could induce reconnections of adjoining vortex lines, and 
eventually a turbulent state. 

Splitting of a multiply-charged vortex 

As seen in Eq. p.lOp . the energy cost to create a g > 1 multiply-charged vortex is unfavorable 
compared to q singly-charged vortices. This raises an interesting question: What happens when 
such an unstable vortex is created? A multiply-charged vortex can be created experimentally in 
atomic BECs by using topological phase imprinting (205) . and causes interesting disintegration 
dynamics [206) (The physics of topological phase imprinting will be described in Sec. 13. 4|) . 

The stability of a multiply-charged vortex in a trapped BEG shows an interesting interaction 
dependence, which originates from the finite size effect |207l . BdG analysis for the cylindrical 
system revealed the appearance of complex eigenvalue modes, which implies that the multiply- 
charged vortex is dynamically unstable. For a vortex with a winding number q, the conservation 
of angular momentum leads to constraint of the normal mode functions Uj{r) = Uj{r)e^^'^^'^^^ and 
Wj(r) = Uj(r)e*^^~^-'^. For q = 2 there are alternating stable and unstable regions with respect to 
the interaction parameter auz — a j \^\'^dxdy\ the first (second) region appears for < ariz < 3 
(11.4 < an^ < 16). Numerical simulations demonstrated that, when the system is in the unstable 
region, a doubly-charged vortex decays into two singly-charged vortices |208| . 

The experiment by Shin et al. investigated the splitting process of a doubly-charged vortex 
and its characteristic time scale as a function of anz=o |206) . The results showed that the doubly- 
charged vortex decayed, but that its lifetime increased monotonically with ariz^o, showing no 
periodic behavior. Numerical simulation of the 3D GP equations also displayed this mysterious 
phenomenon, demonstrating that the detailed dynamical behavior of a vortex along the entire 
z-axis is relevant in the characterization of the splitting process in an axially-elongated condensate 
[2091 1210| . The vortex begins to split near both condensate edges in the z-direction, and then the 
splitting propagates to the center. 

Isoshima et al. studied the splitting dynamics of quadruply-charged {q = 4) vortices, both 
theoretically and experimentally [211] . In this case, a rich variety of splitting patterns with l- 
fold symmetry was expected [212| . They initially prepared an off-centered q — A vortex, and 
observed that the vortex split into an array of four linearly-aligned q — 1 vortices. This is because 
the perturbation with two-fold rotational symmetry (/ = 2) is the most unstable mode under 
the experimentally relevant parameters. These observations were well-reproduced by numerical 
simulation of the GP equation, where the initial displacement of the vortex position was crucial to 
causing splitting instability due to the conservation of angular momentum. 

3.2.7 Vortex lattices 

At very high frequencies VL, a rotating superfiuid mimics rigid body rotation with V x Vg = 2f2, by 
forming a triangular vortex lattice which arranges itself parallel to the rotation axis |213j . These 
properties closely resemble the magnetic fiux-line lattice in type-II superconductors predicted by 
Abrikosov (2141. Here, we summarize the physics of vortex lattices in trapped BECs. Since various 
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features of the rapidly rotating BECs were reviewed thoroughly by Cooper [130) and Fetter {163| , 
in this section we describe only the theoretical basics that will be useful in the following sections. 

Basics 

For vortex lines parallel to the z-axis and located at = [xi,yi), the vorticity is given by the form 
V X -Vs — (ri)z;. Therefore, the average vorticity per unit area is given by V x Vg = n-n.^z, 

where n„ is the number of vortices per unit area. Hence, the density of the vortices is related to 
the rotational frequency U a,s — 2n/K ("Feynman's rule") [TU]. This relation can be used to 
estimate the maximum possible number of vortices in a given area as a function of fl. The vortex 
lattice can be characterized by the nearest-neighbor lattice spacing ^ b = {h/mny/'^ defined by 
the area per vortex = 7r6^, and by the radius of each vortex core Tc ^ 

Thomas-Fermi regime 

The GP energy functional of Eq. (j3.8l) in a rotating frame can be rewritten as 

+ Kff|^P + ||*|') , (3.12) 

where VcS = rn{ijLi\ — r2^)r^/2 + muj^z^ 12 is the effective trapping potential combined with the 
centrifugal potential; the rotation effectively releases the radial potential, and causes it to vanish 
at 17 = w_L. Because the first term in Eq. (|3T2)) reads jlm + to(vs - ^ x r)2|*p/2, it 

can be neglected when both the TF limit and the rigid-body rotation limit = Jl x r are satisfied. 
This TF approximation is valid when the lattice spacing h is much larger than the vortex core size 
Tc, giving an effective coarse-grained description of the rotating BECs. The TF radius is given by 
i?_L(r2) = i?j^/[l — (ri/ojj^)^]^/^" with _Rj^ for a non-rotating condensate, providing an aspect ratio 
Arb = -Rj_(fi)/i?z — A/[l — (r^/wj^)^]^/^. Thus, measurement of Arb can give the rotation rate of 
the condensate |128[ 1165] . Also, in the high rotation limit wj^, the condensate flattens out 

and reaches an interesting quasi-2D regime; current experiments have reached il/wj^ w 0.995 |166) . 

Mean-field quantum Hall regime 

Note that the first term in Eq. p.l2p can be identified as the Hamiltonian = {—ifi^ — 
eA/c)^/2m of a charge — e particle moving in the xy plane under a magnetic field Bz with a vector 
potential A — {mc/e)Q, x r. If interaction is neglected (3 = 0), the eigenvalues of the Hamiltonian 
of Eq. (I3.12P forms Landau levels as n^/^ = + n{uj_\_ + + m{oj± — ft) + (n^ -I- l/2)wz, 
where n is the Landau level index, m indexes the substates within the n-th Landau level, and 
Uz is the index of the states along the z-axis. The lowest energy states of two adjacent Landau 
levels are separated by h(ij± + Q), whereas the distance between two adjacent substates in a given 
Landau level is h{uj± — fl); when $7 = Ld±, all states in a given Landau level are degenerate. This 
corresponds to a situation where the centrifugal force exactly balances the trapping force in the 
x-y plane, and only the Coriolis force remains. The system is then invariant under translation, 
and therefore exhibits macroscopic degeneracy. This formal analogy has led to the prediction that 
quantum Hall-like properties would emerge in rapidly rotating BECs; see [13011131] and references 
therein for the details of this strongly correlated phase in rapidly rotating bosons. 

The quantum-Hall formalism provides a useful mean-field description of rotating BECs with 
vortex lattices [215) . Interaction will lead to the mixing of different (n, m, Uz) states. However, 
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because the averaged density n of the system drops as 17 — > uj±_ , the interaction energy ~ gn can 
become low compared to 2huj± and hujz- In this hmit, particles should condense into the lowest 
Landau levels (LLL) with n = 0. The system then enters the "mean field" quantum Hall regime, 
where the wave function can be described by only the LLL orbitals with the form 



LLL = ^ a,„Vm(l") = ^n(^ - ^j)^^' 
m>0 



(3.13) 



where z — x + iy, Zj are the positions of vortices (zeros), and A is a normalization constant. By 
minimizing Eq. p.l2p using the ansatz ^'lll with respect to Zj, one can analytically describe 
the vortex structure [^11 Hill EUl IUHl ElOl 1221] , although this method is effective only near the 
limit il ~ LOj_. However, this would be impractical for numerical calculations because of the time 
and accuracy required. A similar approach was applied to investigate the ground state of (not 
rapidly) rotating BECs with very weak interaction |222[ 12231 1224) , where the coefficients a„i of the 
harmonic oscillator basis ipm were minimized. 

Schwaikhard et al. created rapidly rotating BECs by spinning up the condensates to Q/uj± > 
0.99 [170j . and found some evidence that the condensate entered the LLL regime: (1) The 2D 
signature of a rapidly rotating BEC was confirmed by the excitation of an axial breathing mode 
{niz = 0). For a BEC in the axial TF regime, an axial breathing frequency wb = V^uJz has been 
predicted in the limit fl/u^ — > 1 [225 , whereas lob — Sa;^ is expected for a noninteracting gas 
with /i < hujz ■ Schweikhard et al. observed a crossover of ojb from VSujz to 2ujz with increasing 
f2 (/z ~ Shujz)- (2) A signature of the LLL regime is that the vortex core size Tc is similar to 
the vortex separation distance b = {h/rnQy/^. The theory P^I^^I^^HI^ predicted that the 
vortex cores begin to shrink as the intervortex spacing becomes comparable to the healing length 
^, and eventually the core radius becomes proportional to the intervortex spacing. The saturation 
of the fractional area A = r'^/b^ with increasing was actually observed, indicating that the 
system was in the LLL regime. (3) The observed global density profile remained a parabolic TF 
profile even for ft — > uj±, although it was expected to be a Gaussian j215j (see Eq. (I3.13p ) in the 
LLL regime. However, very small distortions of the vortex lattice from a perfect triangle could 
result in large changes in the global density distribution, such as from a Gaussian to a TF form 



Tkachenko oscillation of a vortex lattice 

It should be possible to propagate collective waves transverse to the vortex lattice in the superfiuid, 
in so-called Tkachenko (TK) modes [229) . For incompressible superfluids, the dispersion law is 



given by WTK(fc) = ^/hDjAm k. The TK mode of a vortex lattice in a trapped BEC was observed 
experimentally [168] . and was identified by the sinusoidal displacement of vortex cores with an 
origin at the center of the condensate. 

To explain the observed frequency of the TK mode a;(„ where the quantum number (n, m) 
refers to the radial and angular nodes, the effects of compressibility should be taken into account. 



According to the elastohydrodynamic approach developed by Baym [230) . the TK frequency is 
described by the compressional modulus Ci and the shear modulus C2 of the vortex lattice included 
in the elastic energy 
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where e(r, t) is the continuum displacement field of the vortices from their home positions. In 
the incompressible TF regime, C2 = —Ci = nhfl/8. Compressibility yields two branches in 
the energy spectrum: The upper branch follows the dispersion law = 417^ + c^A;^, where 
c = y^gn/m is the velocity of sound, which is the standard inertial mode of a rotating fluid with 
a gap at k — 0. The lower frequency branch corresponds to the TK mode with the dispersion 
= {hil / 4m)[c^ k'^ / (Afl^ + c^fc^)]. For large k, this reproduces the original TK frequency wtk, 
but for small k it exhibits the quadratic behavior aj_ ~ y^h/lGmflck^. The transition between k"^ 
and k dependence takes place at A; ~ fi/c > R^^. This suggests that the effects of compressibility, 
which characterize the fc^ dependence, play a crucial role in the TK mode. Thus, this regime 
is distinguished from the usual incompressible TF regime, and is referred to as the "soft" TF 
regime. Including finite compressibility, this can explain the observed values of W(i,o) |230ll2311l232) . 
First-principles simulations based on the GP and BdG formalism also agreed excellently with the 
experimental data [2331 1234] . 

3.2.8 Various topics on vortices in single-component BECs 
Vortices in an anharmonic potential 

For a rotating condensate with a frequency i7 in a harmonic potential (l/2)ma;^r^, the centrifugal 
potential effectively reduces the confinement, preventing a BEG from rotating at beyond oj±. 
This restriction can be avoided by introducing an additional quartic potential [2351 12361 1237) , so 
that the combined trapping potential in the xy plane becomes Vc^{r) = (l/2)mu;^(r^ + Ar^/a^^), 
where the dimensionless parameter A represents the relative strength of the quartic potential. 
The properties of a rotating condensate in an anharmonic potential have attracted considerable 
theoretical attention |239l 12401 12411 1242] . The vortex phases in an anharmonic trap are quite 
different from those in a harmonic trap, as shown in Fig. 13.41 because it is possible to rotate the 
system arbitrarily fast. For small fi, the equilibrium state is the usual vortex lattice state. As 
increases, the vortices begin to merge in the central region, and the centrifugal force pushes the 
particles towards the edge of the trap. This results in a new vortex state consisting of a uniform 
lattice (multiple circular arrays of vortices) with a central density hole. The central hole becomes 
larger with increasing f2, and the condensate forms an annular structure with a single circular array 
of vortices. A further increase of 57 stabilizes a giant vortex, where all vortices are concentrated in 
the single hole. 

A combined harmonic-quartic potential was constructed by Bretin et al., who superimposed 
a blue detuned laser with the Gaussian profile [129) . Since the waist w of the beam propagat- 
ing along the z-axis was larger than the condensate radius, the potential created by the laser 
C/q exp(— 2r^/w^) can be written as U{r) ~ Uq{1 — 2r^/u>^ + 2r^/'w^). The second term leads 
to a reduction of the transverse trapping frequency wj^, and the third term provides the desired 
quartic confinement, giving uj^/2t: = 65 Hz and A ~ 10^'^ for Vc^{r). For f2 < a vortex lattice 
was clearly observed. When $7 > wj^, however, the vortices became gradually more difficult to 
observe and the images became less clear for Q. = 1.05wj^. The most plausible explanation of this 
observation is that the vortices were not in equilibrium, and were strongly bent |240) . 

Vortex pinning in an optical lattice 

Rotating BEGs combined with a co-rotating optical lattice are an interesting system, which has 
two competing length scales: vortex separation and the periodicity of the optical lattice. The 
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Figure 3.4: Vortex structure of a rapidly rotating BEC in a combined harmonic-quadratic potential. 
This calculation was done by 2D simulations of the GP equation with the parameter values m = 
msTfj^j,, uj± = 10 X 27r Hz, a — 5.61 nm and A = 0.5. The upper and lower panels correspond to 
iV = 10^ and iV = 10^, respectively. 

structure of the vortex lattice is strongly dependent on the externally applied optical lattice. 
Various vortex phases appear, depending on the number of vortices per pinning center, i.e., the 
filling factor. Tung et al. created a rotating optical lattice using a rotating mask |243j . which 
provided a periodic pinning potential that was stationary in the corresponding rotating frame 
of vortices. The authors observed a structural crossover from a triangular to a square lattice 
with increasing potential amplitude of the optical lattice, which was consistent with theoretical 
predictions |244[ 1245] . In a deep optical lattice, the condensates are well localized at each potential 
minimum, so that the system can be regarded as a Josephson-junction array (JJA). A rotating 
Josephson-junction array can have characteristic vortex patterns with the unit-cell structure [246j . 
and was realized experimentally by Williams et al. |247j . 

Thermally-activated vortex generation 

In 2D systems with continuous symmetry, true long-range order is destroyed by thermal fluctu- 
ations at any finite temperature. For 2D Bose systems, a quasi-condensate can be formed with 
a correlation decaying algebraically in space, where superfluidity is still expected below a cer- 
tain critical temperature. This 2D phase transition is closely connected with the emergence of 
thermally activated vortex-antivortex pairs, known as the Berezinskii-Kosterlitz-Thouless (BKT) 
phase transition occurring at T = Jbkt [2481 1249| . For T < Tbkt, there are no isolated free 
vortices; vortices only exist as bound pairs, consisting of two vortices with opposite circulations. 
The contribution of these vortex pairs to the decay of the correlation is negligible, and the alge- 
braic decay is dominated by phonons. For T > Tbkt, the vortex pairs are separated, and the 
free vortices form a disordered gas of phase defects and give rise to an exponential decay of the 
correlation. 

Recently, a BKT transition was experimentally observed in ultracold atomic gases [2501 1251| . 
In the experiment, a ID optical lattice was applied to an elongated condensate, splitting the 3D 
condensate into an array of independent quasi-2D BECs. The interference technique revealed the 
temperature dependence of an exponent of the first-order correlation function of the fluctuating 2D 
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bosonic field |252j . A universal jump in the superfluid density characteristic of the BKT transition 
was confirmed by observing a sudden change in the exponent, in which a finite size effect caused 
a finite-width crossover, rather than a sharp transition. The microscopic origin of this transition 
was clarified from the image of the interference of the two 2D condensates. If isolated free vortices 
were present in either of the two condensates, the interference fringes would exhibit dislocations. 
Such a dislocation was observed in the high-T region of the crossover, as predicted theoretically 
by classical field simulations [253) . 

Vortex formation during a rapid quench 

One universal element of many continuous phase transitions is the spontaneous formation of topo- 
logical defects during a quench through the critical point, known as the Kibble-Zurek mechanism 
|254l 1255) . The rapid symmetry-breaking phase transition in the early universe is expected to have 
produced causally disconnected domains of a true vacuum state, and their phase mismatch would 
leave behind topological defects like cosmic strings or monopoles. Since the cosmological phase 
transition may be second-order, it is especially significant to perform experiments in a superfluid 
system whose transition is also second-order |256[ 1257) . 

Ultracold atomic BECs offer an ideal opportunity for studying the Kibble-Zurek mechanism 
because the temperature, strength of interaction, and external parameters such as the magnetic 
field and trapping potentials can be changed in a time shorter than the characteristic time scale 
of topological defect formation. A quench through the EEC transition temperature was realized 
by using a rapid ramping of the rf-field for evaporative cooling. Weiler et al. studied the sponta- 
neous formation of vortices in EEC both experimentally and theoretically )258) . Vortices appeared 
spontaneously during condensate formation, without the application of any external rotation. This 
experiment was compared to a simulation of EEC growth dynamics, calculated using the stochastic 
projected GP equation [259) , 

dV'(r, t)=V |--^LgpV'(i-, t)dt + - iGp)V'(r, t)dt + dVK^(r, t)| (3.15) 

which was derived from first principles using the Wigner phase-space representation |260j . Here, 
the states of the trapped system are divided into the "coherent" region and the "incoherent" region 
by the high-energy cut-off of the harmonic oscillator mode. The first term on the right describes 
unitary evolution of the classical field in the coherent region according to the GP equation. The 
second term represents growth processes, that is, collisions that transfer atoms from the thermal 
bath to the classical field and vice versa, and the form of 7(r) may be determined from kinetic 
theory [261) . The third term is the complex- valued noise associated with condensate growth, 
which is consistent with the fluctuation-dissipation theorem. The projection operator V restricts 
the dynamics to the coherent region. The simulation results show that, after the temperature 
quench, the condensate domains begin to grow in various regions of the fluctuating classical fleld, 
and several vortices and anti- vortices (like a quantum turbulence) are formed. Eventually, a single 
vortex remains in the well-developed condensate. 

Vortex formation during the interference of BECs 

In a related work, vortices were formed by merging three uncorrelated EECs that were initially 
separated by a triple- well potential [262) . This may verify one side of the KZ mechanism: Can 
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the causally disconnected domains generate a net circular flow, i.e., vortices, through a merging 
process? When the potential barrier created by a repulsive laser beam was removed, the three BECs 
merged. Depending on the relative phases between the condensates and the rate of barrier removal, 
vortices formed stochastically even without an applied rotation. If the condensates have random 
initial phases, a vortex should appear with a probability of 0.25 |2621I263] . which was observed with 
a slow removal rate of the barrier. With quick barrier removal, the resulting interference fringes 
could generate pairs of vortices and anti- vortices via "snake instability" [2641 I265j . but these 
soon disappear due to pair annihilation. These properties were studied by extensive numerical 
simulations of the GP equation [2661 1267) 

3.3 Vortices in two-component condensates 

We now address the properties of vortices in multicomponent BECs. Although there have been 
many theoretical predictions of exotic vortices in multicomponent BECs, much experimental work 
remains to be performed. Here, we primarily describe theoretical studies of the topic, along with a 
few related experimental results [113 [IMl [2M1 [IZQl [271] ■ In this section, we focus on the properties 
of vortices in two-component BECs. 

3.3.1 Multicomponent Bose-Einstein condensates 

Before discussing vortices, we summarize the types of multicomponent BECs and describe how 
they are obtained. In this system, one can load and cool atoms in more than one hyperfine 
spin state or more than one atomic element in the same trap, thereby experimentally realizing 
multicomponent condensates. The simplest and most intuitive example is a mixture of different 
kinds of atoms. However, due to the effect of strong inelastic collisions between different kinds 
of atoms, this example is generally difficult to realize. One successful example is the *^Rb-'*^K 
mixture [144] . Recently, a BEC mixture of isotopes of Rb atoms (^^Rb-^^Rb) was realized by two 
groups [1481 1149) . These experiments demonstrated that the interacting properties of two atomic 
species were tuned via a Feshbach resonance, opening a new avenue to study superfluid mixtures 
in a well-controlled manner. 

A common method used to realize a condensate mixture of ultra-cold atoms is to utilize the 
hyperfine spin states of atoms. The resulting BECs can have internal degrees of freedom, due 
to these hyperfine states. A hyperfine- Zeeman sublevel of an atom with total electronic angular 
momentum J and nuclear spin I may be labeled by the projection mp of total atomic spin F = 
I + J on the axis of the field B and by the total F, which can take a value from \I — J\ to 
|/ -|- J|. This is because the hyperfine coupling, which is proportional to I • J, is much larger 
than the typical temperature of an ultra-cold atomic system. The hyperfine state is denoted by 
\F, mp) with mp = —F, —F -I- 1, • • • , F — 1, F. The simultaneous trapping of atoms with different 
hyperfine sublevels makes it possible to create multicomponent (often called "spinor" ) BECs with 
internal degrees of freedom, characterized by multiple order parameters [1411 11421 11431 11451 11461 
1147) . The order parameter of a BEC with hyperfine spin F has 2F + 1 component with 
respect to the basis vectors \mp) defined by Fz\mp) — mp\mp), which can be expanded as 

I*) = E™,=o,±i,...±f*-^I™f). 

Two-component mixtures of BECs with different hyperfine spin states have been created in 
the laboratory for the systems of ^^Rb [14ll [27l [273] and ^SNa [274]. In Ref. [Ml], the two- 
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component BECs of ®^Rb atoms consisted of the hyperfine spin states \F = l,mp = —1) = |1) 
and \F — 2,tof = 1) = |2). The MIT group used an optical trap, ah spin manifold F — 1 
spin states of ^'^Na (tof = 0, ±1) were simultaneously trapped in a single trap. In this case, the 
spin exchange interactions (see Sec. 13. 4[) gave rise to spin mixing dynamics among the different 
components. They made use of a quadratic Zeeman effect for atoms by applying a magnetic field 
to prevent the \mp = —1) component from appearing. As a result, the system can be regarded as 
a two-component BEC. 

When the trapping potential V^^^{r) {i = 1,2) is produced by a magnetic field, the atoms of 
each hyperfine state have their own potential energy because the potential energy depends on the 
magnetic moment of the atom. For example, the trapping potential for atoms in the \Fi,mF.) 
state is expressed as a function of the magnetic field |-B(r)| as |275) 



^^(r) = gF.fj.BmF,\B{r)\ ~ gF^HBiriF^ 
1 



Bo + ^{K,x^ +Kyy^ +K,z^) 



= V^ + -m,{^ix'+^fyy'+ujtz^), (3.16) 

where g^. is the f;-factor of the i-th atom, fis is the Bohr magneton, Vq — gFifJ-B'm-FiBo, and 
^ik {k — x,y,z) is the trapping frequency satisfying the relation niiLuff, — gFilJ^BfnPiKk- The 
potential minima Vq' for BECs with different values of gF^mFi do not coincide. This positional 
shift is also influenced by a gravitational effect, the difference in nuclear magnetic moments, and 
nonlinearity of the Zeeman shift |140[ 1141] . If the potential is created by a polarized optical laser 
|1421 11431 11451 1146| , atoms in all hyperfine states share the same confining potential. 

3.3.2 Coupled Gross-Pitaevskii equations 

First, we describe the mean-field theory for two-component BECs, which is a natural extension 
from that of single-component BECs. The condensates are confined by trapping potentials V^J^r) 
(i = 1,2), which are assumed to rotate at a rotation frequency about the z axis as H = ilz. The 
CP energy functional of the two-component BEC in the rotating frame is 
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Here, is the mass of the i-th atom. The coefficients gi and gi2 represent the atom-atom inter- 
actions between atoms in the same component, and atoms in different components, respectively. 
These are expressed in terms of the corresponding s-wave scattering lengths oi, a2, and ai2 as 
gi = infi^ai/mi, and 512 — 27r/i^ai2/mi2, where = m^^ + "^2^ is the reduced mass. The 
time-dependent coupled CP equations for two-component BECs can be obtained using a variational 
procedure ihd'^i — 5E/5^* as 

(^-^ + Kx(r)+5d*.P+.9i2|*3-.p-f^^.j*., (3.18) 
Substituting 4'i(r, t) = $i(r)e~*'^'*/'' yields the time-independent coupled CP equations 

4^ + V:^ir) + + gi2|$3-z|' - = (3.19) 

2m,: / 
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where we have introduced the Lagrange muUipher /x^, which represents the chemical potential 
and was determined so as to satisfy the conservation of particle number Ni = J dr\^i\^ for each 
component. 

It is instructive to examine the ground state without rotation (f2 = 0). Compared to single- 
component BECs, two-component BECs exhibit a rich variety of ground state structures, depending 
on the various parameters of the system [275] . The intercomponent interaction gi2 plays an 
especially important role; phase separation of homogeneous two-component BECs occurs when 
9i2 > v^5iff2 |276l 1277) . where strong intercomponent repulsion overcomes the intracomponent 
repulsion. In *^Rb, the scattering lengths are nearly equal (oi = 5.67 nm, 02 = 5.34 nm, and 
ai2 — 5.50 nm), thus .g^2 — \/9i92, but the small relative displacement of the trapping minima leads 
to a spatial separation between the two components [141) . Even without the trap displacement, 
the two components phase separate into oscillating ring structures |273| . In contrast, for the 
^■^Na-mixture, oi = 2.75 nm, 02 = 2.65 nm, and 012 = 2.75 nm, thus 51^2 > 1/3152, causing a 
well-separated spin domain to form |142[ I274| . A recent experiment demonstrated the control of 
the miscibility of a *^Rb-*^Rb mixture via Feshbach resonance |149j . 

3.3.3 A single vortex in a two-component BEC 
Experimental creation 

The first observation of a quantized vortex in an atomic BEC was achieved in a two-component 
BEC consisting of ^^Rb atoms with hyperfine spin states \F = 1, mp = —1) = |1) and \F = 2, itif = 
1) = |2) |118| . which were confined simultaneously in almost identical magnetic potentials. Initially, 
condensed atoms were trapped in one state, such as the |1) state. Then, a two-photon microwave 
field was applied, inducing coherent Rabi transitions of atomic populations between the |1) state 
and the |2) state. For a homogenous system in which both components have uniform phases, the 
interconversion takes place at the same rate everywhere. However, time variation of the spatially 
inhomogeneous potential changes the nature of the interconversion. Let us consider a co-rotating 
frame of an off-centered perturbation potential at the rotation frequency Q' . In this frame, the 
energy of a vortex with one unit of angular momentum is shifted by Ml,' relative to its value in the 
laboratory frame. When this energy shift is compensated for by the sum of the detuning energy 
of an applied microwave field and the small chemical potential difference between the vortex and 
non- vortex states, a resonant transfer of population can occur between the non- vortex |1) state and 
the vortex |2) state. The rotating perturbation created by a spatially inhomogeneous laser beam 
was rotated around the initial non-rotating component By adjusting il' and the detuning, 
the 1 2) component was resonantly transferred to a state with unit angular momentum by precisely 
controlling the time at which the coupling drive was turned off |278j . This procedure resulted in a 
"composite" vortex, where the |2) component had a vortex at the center, while the non- rotating 
|1) component occupied the center and acted as a pinning potential to stabilize the vortex core. 

It is possible to prepare the initial condensate into either the |1) or |2) state, and to then form 
a vortex in the |2) or |1) state. It was observed that the stability of the vortex states was strongly 
dependent on the atom-atom interactions [118) . The vortex was stable only when the vortex was 
in the |1) state and the core was in the |2) state. The other state, with the vortex in the |2) state, 
exhibited an instability, with the |2) vortex sinking inward toward the trap center and breaking 
up. The details are discussed below. 
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Equilibrium structure 

Here, we present theoretical studies on the structures of a single vortex state in two-component 
BECs. The analysis is based on the coupled GP equations p. 191) . As observed in Rcf. [118| . 
the vortex structure consists of one circulating component that surrounds the other non-rotating 
component, where both wave functions have axisymmetric profiles. We assume that the condensate 
wave functions have one quantized vortex at the center as ^^(r) = (f>.i{r, z)e^'''^ , where Ui — <j)f is 
the condensate density and the trapping potential is given by V^^{r) — mi{ujfj_r'^ +ujf^z'^) /2. Thus, 
the axisymmetric (singly-charged) vortex states are characterized by {qi, 52) — (1, 0), (0, 1), (1, 1). 




Figure 3.5: Stationary solutions of two-component BECs with winding number {qi,q2) — (0,0), 
(1,0), and (0,1). The parameters of the ^^Rb atoms in |1,— 1) and |2, 1) states were used. The 
trapping potential was oblate as uj± —5x2?: Hz and ujz — 10 x 27r and the particle number 
Ni — N2 = 10^. In (a)-(c), the radial density profile is shown. The corresponding (pseudo)spin 
textures in the z — and y — planes are shown for (d) (1, 0) and (e) (0, 1) states. 

For ®^Rb atoms in non-rotating traps, in order to decrease the intracomponent mean field 
energy, the component with the larger intracomponent scattering length ai forms a shell outside 
the ^'2 component occupying the central part, as shown in Fig. I3.5f a). For (91,(72) = (IjO) [Fig. 
I3.5r b)]. a density depletion associated with the vortex core appears in the density of ^'i along the 
z-axis. The ^'2 component then fills the vortex core, because of intercomponent repulsion. As a 
result, the characteristic size of the core becomes larger than the healing length ^ = h/ ^2migini 
of the single component. In this case, the centrifugal force associated with the vortex causes the ^'i 
component to expand radially, which allows a decrease in the intracomponent mean-field energy 
of the '^i component, rather than that of the nonvortex state. 

On the other hand, the structure of (gi, 92) = (0, 1) is different from the (gi, 52) — (1, 0) solution 
[see Fig. I3.5f c)]. The core size becomes smaller than that of the (1,0) solution. While a fraction 
of the component fills the vortex core, the excessive Vfi component extends outward, resulting 
in an increase in size of the region of coexistence with the rotating \E'2 component. As discussed 
below, this configuration is dynamically unstable |279[ 1280] . Also, the (91,92) = (1,1) state is 
unstable because the two components fully overlap, which is energetically unfavorable. In the 
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presence of a slow rotation, the thermodynamically stable configuration of (1,1) is that in which 
the two vortices in each component are displaced from the center, which decreases the size of the 
overlapping region. 

It is instructive to represent the spinor order parameter of two-component BECs using the 
"pseudospin" [2811 12821 1283) . which allows us to analyze this system as a spin- 1/2 BEC. We 
introduce a normalized spinor Q to represent the two-component wave functions as [5'i,^'2]"^ — 
where 



Ci 

C2 



cos(6l/2)e-*'^/2 
sin(6'/2)e^'^/2 



(3.20) 



and ICiP + IC2P = 1- Here, we assume that ^1 and 4*2 represent the up and down components of 
the spin-1/2 spinor, respectively. Hence, the four degrees of freedom of the original wave functions 
= (their amplitudes nj and phases dj) are expressed in terms of the total density 

nx — ni + n2, the total phase O = 9i + 62, and the polar and azimuthal angles (0, ip) of the local 
pseudospin S ~ {Sx, Sy, Sz), defined as 



Cl*C2 + C2C1 

-*(Cl*C2-C2*Cl) 

IC1P-IC2P 





sin 9 cos if 




sin sin ip 




cos 9 



(3.21) 



where cr is the Pauli matrix, cos 9 = (ni — ri2)/riT, (p = 92 — 9i, and |Sp = 1. 

In the pseudospin picture, the axisymmetric vortex is analogous to a spin texture, extensively 
studied in superfluid '^He. When the wave function for the (1, 0) case is parameterized as (^"1, ^'2) = 
y^riT(e*^ cos(/3(r)/2), sin(^(r)/2)), the configuration satisfying the boundary condition /3(0) — n 
and /3(oo) = is referred to as an Anderson- Toulouse vortex |284| . and /3(cxd) = tt/2 is a Mermin- 
Ho vortex [285) . The spin texture for the (1,0) solution is shown in Fig. 13.51 (d). Here, the ^I*! 
component vanishes at the center, so that the pseudospin points down at r = (/3(0) — tt) according 
to the definition of Eq. (I3.2ip . The spin aligns with a hyperbolic distribution as {Sx,Sy) oc {—x,y) 
around the center in the x-y planes. At the edge of the atomic cloud, the ^2 component vanishes, 
and the pseudospin points up. In between, the bending angle /3(r) — cos^^ Sz decreases smoothly 
from TT at r = as r increases, as seen in the profile in the y ^ plane. For the (0, 1) state, the 
texture exhibits a radial-disgyration as shown in Fig. 13.51 (e), where the spin at x = y = points 
up and aligns as {Sx, Sy) oc (x, y) in the x-y plane. Note that the spin rotates 7r/2 as going radially 
outward from the center. This configuration corresponds to a Mermin-Ho vortex. In the case of 
superfiuid ^He, an MH vortex is stabilized by the boundary condition of the 1-vector imposed by 
a cylindrical vessel. However, in an atomic-BEC system there is no constraint at the boundary; 
the value f3{r) at the boundary r = R should be determined self-consistently. The value /3{R) 
can change arbitrarily with D, by varying the ratio Ni /N2 or gi / g2 1283) , which implies that an 
intermediate configuration between a Mermin-Ho vortex (/3(i?) = 7r/2) and an Anderson- Toulouse 
vortex {f3{R) — 0) can be thermodynamically stable. 

In terms of the spin-1/2 BEC, since the density ut does not vanish at the vortex core, these 
vortices can be called careless vortices. We can define an effective velocity field for the spinor BEC 
as 



Vcff 



2im 



i=i,2 



2m 



(Ve-cos6iV(p), 



(3.22) 



which depends on the gradient of the total phase 9 = 6*1+^^2 and that of the angle of the pseudospin. 
By neglecting the singular contribution of the velocity due to the total phase 8, we can obtain the 
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vorticity 

u;eff = V X Veff = —(V9) X (sineV ip) = — e„«^S'aVS'« x VSy (3.23) 
2m 4to 

with the Levi-Civita symbol eajSj- The last equality is known as the Mermin-Ho relation |285) . 
The vortex was characterized by the topological charge (or the Pontryagian index) 1286) 

Q = J d?rq{v) J d^ru.fi, = ^J (fr {ea^p^Sa^V x VS^)^ , (3.24) 

where q{r) is the topological charge density. The topological charge density q{r) characterizes the 
spatial distribution of the vorticity of the coreless vortex, which is distributed around the center, 
and |vcft| vanishes at the center, contrary to the case of a conventional vortex in a single-component 
condensate. 




Figure 3.6: A simulated (1,1) state, referred to as a vortex molecule |288| . The parameters of the 
^^Rb atoms in |1, —1) and |2, 1) states were used. The trapping potential is oblate uj± = 5 x 27r Hz 
and uJz — 20 X 2tt and the total particle number = 2 x 10^. To stabilize the vortex molecules, 
the rotation il = 0.2wj^ and the Rabi coupling with ujr = 0.1aj_L is applied. Then, the balance of 
the population (A^i ~ N2) is controlled by the detuning A = 0.08, which denotes the difference 
between the frequency of the rf-field and the energy difference between the two internal states, (a) 
Cross section of the condensate density along the x-axis. (b) Grayscale plot of the relative phase 
around the vortex molecule in the z = plane, (c) The spin textures in the z = Q and y = 
planes. 

It is interesting to consider the nonaxisymmetric (1,1) state in terms of the pseudospin, where 
the pair of coreless vortices can be identified as a pair of Mermin-Ho vortices (meron-pair) [285! , 
because the spins are oriented along the x-y plane, far from the two vortices. An example is shown 
in Fig. 13.61 Note that the boundary condition for the pseudospin is not well-defined in our system. 
If we introduce the internal Josephson coupling between the two components, the spin should 
be oriented along the x-axis, far from the defect. This can be achieved by applying an external 
driving field that couples the two hyperfine levels of the atoms. When the atoms are condensed, 
coherent Rabi-oscillation between the two components occurs [287] . If the field strength is increased 
gradually from zero and its frequency is gradually ramped to resonance, a stationary nearly- weight 
superposition of the two components can be obtained. This Josephson coupling energy, given by 
Ej = — ^p{(5'^\I'2 -t- c.c) — -'2hujB.y/nin2 cos (p — —itku^riT^Sx with the Rabi frequency wr, acts 
as an external field that aligns the spin along the a;-axis. Therefore, the meron-pair is stabilized 
|288j . In the presence of the Josephson coupling, the two vortices are bounded by the domain wall 
of the relative phase [289], shown in Fig. I3.6f b). whose binding energy can be controlled by the 
Josephson coupling strength. 



53 



Vortex stability and dynamics 

Unstable vortex dynamics were experimentally observed to depend on which component had a 
vortex |118j . As seen in the ground state configuration of Fig. I3.5f a). if a component with smaller 
gi is displaced from the trap center, it will tend to return to the trap center. This effect is sometimes 
called the "buoyancy" effect in trapped two-component BECs |197] . and plays an important role 
in the vortex dynamics. 

Quantitatively, the instability can be understood from the BdG analysis of the vortex state 
|279l 1280] . Numerical simulations showed that for equal populations Ni ~ N2 — N, and with the 
coupling constants gi > gi2 > (?2, the axisymmetric vortex (1,0) state is stable, while the (0,1) 
state is dynamically unstable. Among the normal modes of the (1,0) case, there is a negative 
eigenvalue. This mode, which is a dipole-type mode, belongs to a perturbation that causes an 
initial displacement of the vortex from the trap center. However, the drift instability is affected 
by dissipation; the configuration is dynamically stable without dissipation, and furthermore, it 
can be energetically stabilized by rotation |283j . In the (0, 1) case, on the other hand, there are 
normal modes with complex frequencies. The unstable modes have a similar shape as the negative- 
eigenvalue modes of the (1,0) case, which means that the perturbations also displace the vortex 
from the center. The imaginary part of the eigenvalues implies that vortices with unit charge 
in 1 2) are dynamically unstable under a generic perturbation of the initial configuration and this 
instability grows without dissipation. This result is consistent with the JILA experiments |118j . 

To study the dynamics beyond the linear stability, one begins with the time-dependent coupled 
GP equations p.l8p . Numerical simulations of the nonlinear vortex dynamics caused by large 
perturbations were performed by Garcia- Ripoll and Perez-Garcia |279| . The linearly stable state 
(1,0) is robust and survives even under a wide range of perturbations; the vortex only shows a 
precession around the center. In contrast, the unstable configuration (0, 1) gives rise to recurrent 
dynamics: the density of the |1) component and the vortex in |2) oscillate synchronously. These 
oscillations grow in amplitude until the vortex spirals outward. Since Eq. (j3.18l) conserves the 
total angular momentum (L^) = / (ir(5'^L^4'i -I- 5'2L2^'2), the vortex is transferred from |2) to 
|1). Although the dynamics are not completely periodic, this mechanism exhibits some recurrence, 
and the vortex eventually returns to |2). The unstable dynamics of |2) with a multiply-quantized 
vortex such as (91,92) = (0^2), (0,3) were investigated by Skryabin [280] . 

The effect of the intercomponent interaction was also demonstrated by the precession of an 
off-centered vortex |199) . The experiment observed that the precession frequency of the filled-core 
(coreless) vortex was slower than that of the empty core. The slower precession of filled cores can 
be understood in terms of the buoyancy effect. Because of its slightly smaller scattering length, the 
1 2) component has a negative buoyancy with respect to the |1) component, and consequently tends 
to sink inward toward the center of the condensate. With increasing amounts of the |2) component 
in the core, the inward force on the core begins to counteract the effective outward buoyancy of 
the vortex (described in Sec. I3.2.6p . resulting in a reduced precession velocity. It was predicted 
that with a filling component of sufficiently negative buoyancy in the core, the core precession may 
stop, or even precess in a direction opposite to that of the fluid fiow |197] . 

3.3.4 Vortex lattices in two-component BECs 

Vortex lattices in rotating two-component BECs are more complicated than those in single- 
component condensates because the system has two-component order parameters and are coupled 
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by intercomponent interactions. To simplify the problem, each component has an equal number 
of bosons, and the trapping potentials of the two components are identical. Then, the two com- 
ponents will have the same size and the same density of vortices. In this case, one would expect 
that each component will contain identical but staggered vortex lattices, with one lattice displaced 
relative to the other because of the intercomponent repulsive interaction. We will focus on the case 
of gi = g2 7^ <7i2, and vary the parameter 5 — gi2/ y/gigi- As described above, the condition of 
phase separation for non- rotating condensates is given as 5 > 1. The presence of a vortex lattice 
naturally modulates the density of each component, with the high density regions of one compo- 
nent overlapping with the low density regions of the other (the vortex core is filled by the other 
component). Thus, the system is effectively phase-separated whenever staggered vortex lattices 
are present, even for 5 <1. 

The first theoretical study of the lattice structure was done by Mueller and Ho within the 
mean-field quantum Hall regime [290] . extending the analysis in Sec. 13.2.71 to a two-component 
system. In this analysis, (i) the rotation frequency was assumed to be very close to the radial 
trapping frequency (J7 ~ wj^), (ii) the vortex lattice was assumed to be perfectly regular, (iii) the 
range 5 > 1 could not be described, because a non-periodic structure appears as shown below. For 
5 < 0, the two components prefer to overlap and form a single triangular lattice. As 5 increases, 
the lattices undergo several structural transitions. There are five distinct phases of the lattice: 
triangular, displaced triangular, rhombic, square, and rectangular with increasing 5. Kegeli and 
Oktal calculated the dispersion relation of the TK mode for each lattice type, and discussed its 
relation to the structural phase transition [291] . 
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Figure 3.7: (a) fl-S phase diagram for vortex states in rotating two-component BECs, obtained 
by 2D numerical simulations of Eq. (j3.19p : A: triangular lattice, o: square lattice, o: vortex 
sheet; see Ref. '292' for details. Because of the continuous change between triangular and square 
lattices, their boundary is shown by a combination of A and o. The plots at = 1 show the 
results obtained by Ref. |290j based on the LLL approximation, (b) Typical vortex states in each 
phase for fl = 0.8ujj_ and 6 = 0.5, 0.8, 1.1. 

Beyond the above restrictions (i)-(iii), one needs full numerical calculations of the coupled GP 
equations p.l9p [252^ . Figure 13.71 shows the numerically obtained phase diagram of vortex states 
in the {6,fl) plane. In the overlapping region S < 1, two types of regular vortex lattices were 
obtained. For 6 = 0, the two components did not interact. Therefore, triangular vortex lattices 
formed, as in a single component BEG. As S increased, the positions of the vortex cores in one 
component gradually shifted away from those of the other component and the triangular lattices 
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became distorted. Eventually, the vortices in each component formed a square lattice rather than a 
triangular one. The two vortex lattices were interlaced in such a manner that a peak in the density 
of one component was located at the density minimum of the other. The stable region of the square 
lattice depended not only on 6, but also on f2. While an increase in S did cause the deformation of 
the lattices from triangular to square, the transition occurred at a significantly higher value of 6 
than that of the LLL result {S — 0.373 [290) ): for example, the transition occurred at (5 ~ 0.65 for 
fl — 0.7. This implies that an increase in rotation frequency causes a transition from triangular to 
square lattices. The excitation spectrum of these equilibrium states were calculated numerically 
by Woo et al. [293j . showing complicated behavior consisting of the TK modes, the hydrodynamic 
modes, and the surface modes. 

When S exceeds unity, the condensates undergo phase separation to spontaneously form do- 
mains having the same spin component. Concurrently, vortices begin to overlap. In the strongly 
phase-separated region 5 > 1 (c2 < 0), the domains of the same spin component, at which the 
other-component vortices are located, merge further, resulting in the formation of "serpentine" 
vortex sheets |292j . A typical example is shown in Fig. 13.71 Singly-quantized vortices line up in 
sheets, and the sheets of components 1 and 2 are interwoven. The sheet distance can be estimated 
from the competition of the kinetic energy of the superflow and the surface tension of the domain 
wall. The detailed discussion of the vortex sheet was reported in Ref . [2941 . 

Schweikhard et al. observed an interlaced vortex lattice in a two-component BEC consisting of 
two hyperfine levels of the *^Rb atom: |1) and |2) [269) . They initially created a regular triangular 
vortex lattice in a BEC in the |1) state. Then, a short pulse of the coupling drive transferred a 
80 — 85% fraction of the population into the |2) state. After a variable wait time, they took images 
of the condensates, revealing interesting dynamics of the structural transition of the lattices. For 
the first period 0.1 — 0.25 s, there was little dynamical behavior in either component, and certainly 
no structural transition in the vortex lattice. From 0.25 — 2 s, turbulent behavior appeared in 
both components in which vortex visibility degraded significantly, exhibiting a transition from 
overlapping hexagonal vortex lattices to interlaced square lattices. From 2 — 3 s, square lattices 
emerged from the turbulent state. From 3 — 5.5 s stable square lattices were observed in both 
components. At this stage, despite the large (80 — 85 %) initial population transfer to the state 
1 2), the number of |2) atoms decreased rapidly because of a larger inelastic scattering loss of the 
1 2) state than of the |1) state. As the |2) state population continued to decay, the vortex lattice 
reverted to a hexagonal lattice in the |1) state. 

3.4 Vortices in spinor condensates 

An optical trap removes the restriction of the confinable hyperfine spin states of atoms, thus 
realizing new multicomponent BECs with internal degrees of freedom, referred to as "spinor BECs" 
|1421ll43lll451 11^1147) . In contrast to the two-component BECs, interatomic interactions allow for 
a coherent transfer of population between different hyperfine spin states (spin-exchange collisions), 
which yields a fascinating physics of both ground-state properties and spin dynamics. Recently, 
an excellent pedagogical review appeared, addressing various problems on spinor BECs |295j . In 
this section, we focus on the studies of quantized vortices in such spinor BECs with F ^ I and 
F = 2 hyperfine spin states. 
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3.4.1 General formulation of spin-F Bose systems 

To describe the zero-temperature properties of spinor BECs, a generalized GP mean-field model 
was introduced by Ohmi and Machida |296j and Ho [297] for F = 1 spinor BECs. Here, we consider 
the most general case: spin-i^ bosons under a magnetic field B. The second quantized Hamiltonian 
of the system is H — Hq + Hint ■ The single-particle part Hq can be written as 

Ho^ dr J2 ^l(r)(/i + ffF^BB-F„^)v;^(r), (3.25) 

•' a,(3=-F 

where V'alr) is the field annihilation operator for an atom in the hyperfine state \F^a) {a ~ 
~F, ~F + 1, ■ ■ ■ , F — l,F), satisfying the commutation relation [^/'^(r), ■(/'^(r')] — Sai3S{r — r') 
and [-00. (r), ^^(r')] = 0, and h = — fi,^V^/2m -I- V^xlr) with an optical potential Kix(r) that 
can confine the atoms independently of their hyperfine spin states. The spin matrix Fap = 
[{Fx)ai3, {Fy)ai5, (-^2)0/?] is the spin-F representation of the SU{2) Lie algebra. 

The interaction Hamiltonian Hi^t — J drirfr2VF(ri — r2) is governed by a two-body interaction 
Vpiici — r2) that is invariant under spin rotation and preserves the hyperfine spin of the individual 
atoms. In a low-energy limit, this is of the form 

2F 

Vpivi - r2) ^S{n~r2)J2 9 f'Pf, (3-26) 

where <?/ = Ai:h^af/m for atoms with an s-wave scattering length ay corresponding to a chan- 
nel with a total hyperfine spin / of two colliding atoms. The operator Vf projects the pair of 
atoms to a state with total hyperfine spin /, written in the second quantized form as 7^/ = 
Y.L=~fO\^a^f,a, where Of. a = Lai.aa (/"l^' "15 "2)'i/'ii V'ca with the Clebsch-Gordon coeffi- 
cient (/a|F, ai; F, 0:2), forming a total spin / state from two spin-F particles. The sum is taken 
for even numbers of /, because only symmetric spin channels are allowed for bosons. 

3.4.2 Mean-field theory for spin-1 BECs 
Generalized GP equation 

For the spin-1 case, the field operator has three components 
be expressed as 

1 1 

^^=l = Y E V'lV'^V'/sV'c + I E i^ii^l^c.a'-Vpp4p4a'. (3.27) 

where we have omitted (5(ri— r2). The collision coefficients are written as g„ — {go+2g2)/5 and — 
{92 — 5o)/3. The term proportional to Qn is symmetric in the spin indices and represents the spin- 
independent contact interaction. The term proportional to gs, on the other hand, is spin-dependent 
and represents the short-range spin-exchange interaction. The sign of gs determines the nature of 
the spin-exchange coupling: a negative (positive) gs represents ferromagnetic (antiferromagnetic) 
coupling. The spinor BEG of F = 1 ■^■^Na atoms was demonstrated to have an antiferromagnetic 
interaction |142j . while the F = 1 ®^Rb condensate has a ferromagnetic interaction [143) . Although 
9n S> gs in most systems, the t^j-term plays a crucial role in determining the properties of spinor 
BEGs. 



(i/'i , V'o I V"- 1 ) ■ The interaction can 
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In the zero-temperature limit, the total Hamiltonian reduces to an energy functional by taking 
the expectation value E ~ (H) and introducing the condensate wave function = (^a); which 
amounts to neglecting the effects of the non-condensed gas. For the condensates trapped by an 
optical potential and in a uniform magnetic field B = Bz, we can obtain 



1 



.Q = -l 



Y: nh^.-p{F.)+q{F!) + 3-^n' + '-^\S 



(3.28) 



where nT(r) — J2a=-i l^ct('')P ^^'^ total particle density, and S = {{Fx), {Fy), {F^)) is the spin 
density vector defined by {Fi{v)) = Yl\t i3=-i^*a{^){Pi)ap^ p{^) (* = 2;,?/, z). In addition to the 
linear Zeeman term with p = —gp^sB, we have introduced an additional quadratic Zeeman term 
characterized hy q = {gp^sB)^ / E\i{ with the hyperfine splitting Ei^i; these coefficients can be 



varied arbitrarily in experiments |142) . and play an important role in determining the properties 
of spinor BECs. 

The time evolution of the wave function is governed by ihO"^ a{^) / dt = 5E/5'^*^; substituting 
Eq. p.28p into the right-hand side gives 

1 

ih—^ = h^o. + 5„nT*a + ^ [-p{F,)o.p + q{F^)o,p + 5«S • F„^] ^p, (3.29) 

/3=-l 

which are the multicomponent GP equations that describe the mean- field properties of spin-1 BECs 
[HUMS- In a stationary state, we substitute *a(r,t) $a(r)e-'''*/'* into Eq. (P?^ to obtain 

1 

ap + 5sS • Fa^] ^fj — fl^a- (3.30) 

The chemical potential /i is determined so as to conserve the total particle number N = J drnT{r). 
In addition, the spin dynamics governed by Eq. (j3.29p conserves the total magnetization M — 
J dr{Fz{r)) , because the interactions conserve the total spin of two colliding atoms. Therefore, 
when we search for the ground state at a given magnetization A/, we must replace E with E — AM, 
where A is a Lagrange multiplier, p is then replaced with p = p + X, which is determined as a 
function of M. 



Ground state properties 

It is easy to identify the ground state of a uniform system when the Zeeman terms are negligi- 
ble {p = q = 0), because we can set the kinetic energy and t4x(r) to zero. From the energy 
functional (I3.28p . the ground state is ferromagnetic (|S| = tt-t) for 5^ < 0, and polar or antifcr- 
romagnetic (|S| = 0) for gs > 0. Introducing the normalized spinor ^ as ^{r) — y/ nT(r)C(r) = 
•\/nT(r)(C+(r), Co(r), C- (i"))"^: where \C\ = 1, one can write the order parameter of the ferromagnetic 
phase as y^nT(l, 0, 0)"^, while that of the polar (antiferromagnetic) phase is given by y/Ti{0, 1, 0)"^ 
(\/'^t/2(1, 0, 1)-^) [2961 1297] . Here, the polar and antiferromagnetic state are related to the n/2- 
rotation of the spin vector about the j/-axis, which is degenerate in the absence of a magnetic field. 
Because of the symmetry of the Hamiltonian, all states obtained by the gauge transformation e*-^ 
and the spin rotation U {a, /3, 7) = e~'^^"e~*^!'^e~*^'='*' of these ground states are degenerate, where 
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(a, P, 7) are the Euler angles. Thus, the general forms of the order parameter can be written as 



for the ferromagnetic and polar states, respectively. The symmetry group of the ferromagnetic 
state is iS'0(3), since the distinct configurations of C, are given by the full range of Euler angles. 
Note that the phase change by x S'Hd the spin rotation by —7 play an equivalent role, which 
represents a "spin-gauge" symmetry [298]. That of the polar state is U(\) x S'^ jZ^-, where UiX) 
denotes the phase angle x ^-nd S"^ is a surface of a unit sphere denoting all orientations (a, /3) of 
the spin quantization axis. Moreover, a discrete symmetry Z2 indicates that 'J'P°''*'' is invariant 
under tt gauge transformation combined with tt spin rotation around an axis perpendicular to the 
quantization axis n = (cos a sin /3, sin a sin /3, cos i.e., X X + ^-nd n — > — n [299) . 



Figure 3.8: -p-q phase diagram of a uniform spin-1 BEC with (a) 5s > and (b) < 0. For (a) 
Qs > 0, the ferromagnetic phase and the polar phase are separated by a clear boundary \p\ = q + gs, 
which indicates that the \l/o and ^±1 components are immiscible. For small bias fields, with 
q < gsTiT/^ and \p\ < ggTiT, there is a region (IV), where the ^±1 components are mixed and the 
ratio of the populations is given by |\E'+p/|\I'^| = {2gs+p)/{2gs—p). The boundary between (III) 
and (IV) lies at \p\ = 2y^. For (b) gs < 0, The polar state occurs for \p\ < \/q{q — i\gs\), and the 
ferromagnetic state occurs for \p\ > q. Here, all three components are generally miscible, because 
the polar state is skirted by a region of a broken- axisymmetry phase (V) where all components 
have non-zero population. In (c) , the in-situ observation of spontaneously generated spin vortices 
in a ferromagnetic spinor BEC |270] is shown, where a magnetic field was quenched along the 
arrow in (b). The picture represents the profile of the magnitude of the transverse magnetization 
A, its orientation (f), and the magnitude of the longitudinal magnetization Ai [Sadler et al: Nature 
(London) 443 (2006) 312, reproduced with permission. Copyright 2006 by the Nature Pubhshing 
Group] . 

In the presence of an external magnetic field, the ground state becomes more complicated due 
linear and quadratic Zeeman effects. The ground-state phase diagram in a parameter space of 




(3.31) 




A 
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{p,q) was studied in Ref. 1142) . Since gnriT can be absorbed into fi by defining fl = ii — gnnT, 
the ground state corresponds to the spinor, which minimizes the spin-dependent portion of the 
energy density Eg/V = gs\S\'^ /2 — p{Fz) + q{F^). A phase diagram of the ground state is shown 
in Fig. 13.81 (a) and (b); the details of which are described in Ref. f295[ 1300] . Briefly, there are five 
distinct states (I) - (V) depending on the values of {gs,p,q): (I) and (II) are ferromagnetic phases, 
whereas (III) is polar. In (IV), the longitudinal magnetization depends on p as Sz = p/gs- In 
(V), the magnetization is tilted from the z-axis and its azimuthal angle is spontaneously chosen. 
This is referred to as the broken-axisymmetry phase [300) . The rotational symmetry of the order 
parameter about the magnetic-field axis is broken in states (III) and (V). 

3.4.3 Topological vortex formation in spinor BECs 

The internal degrees of freedom in BECs provide a peculiar way to generate vortices, which is 
called "topological phase imprinting" [3011 1301 ^M- The MIT group 1103 IMH] and the 
Kyoto university group |211[ 1305) used this method to create a vortex in a trapped BEC. In a 
magnetic field, atoms can be trapped by the Zeeman interaction of the electron spin with an inho- 
mogeneous magnetic field. Atoms with electron spins parallel to the magnetic field are attracted 
to the minimum of the magnetic field (weak- field seeking state), while those with an electron spin 
antiparallel are repelled (strong- field seeking state). In the experiment at MIT, ^'^Na BECs were 
prepared in weak-field seeking states, a to^) = |1, —1) or |2, -1-2) hyperfine state, and confined 
in a loffe-Pritchard magnetic trap, described by B = i?^(a;x — yy) + Bz{t)i. An essential idea in 
topological phase imprinting is azimuthally dependent adiabatic inversion of the bias field (t) . 

To interpret the mechanism, let us consider a BEC with F = 1. The basis vectors in the 
representation (\I>i, ^o, ^-i) are {|±), |0)}. We introduce another set of basis vectors |x), \y) and 
\z), which are defined by Fj;\x) = Fy\y) — Fz\z) — 0; these vectors are related to the previous 
vectors as |± 1) = =f(1/-\/2) [\x) ± i\y)) and |0) = \z). When the z-axis is taken to be parallel to the 
uniform magnetic field, the order parameter of the weak-field seeking state takes the form -i = V 
and = ^1 = 0, or in vectorial form as '4' = (ip/V^) (x — iy). When the magnetic field points in 
the direction B = i?(cos a sin/3, sina sin/3, cos a rotational transformation with respect to the 
Euler angle (a,/3,7) gives 'S' = (V'/-\/2)e~*'''(m -I- in), where m = (cos a cos /3, sin a cos /3, — sin /3) 
and n = (sin a, — cos a, 0). The unit vector 1 = rh x n = (— cos a sin/3, — sina sin /3, — cos /3) 
corresponds to the direction of the spin S, which is antiparallel to B. The three real vectors 
{1, m, n} form a triad, analogous to the order parameter of the orbital part of superfiuid '^He. The 
same amplitudes in the basis {|0), |±)} are ^-i = (?/'/2)(l-cos/3)e"*"-*'', *o = -(-0/^2) sin /3e-'^, 
and = (VV2)(1 + cos/3)e*"-*'^. 

When the field is strong compared to the transverse quadrupole field B'j^r, the trapped 
condensate initially has an order parameter ^ — {i/j/\/2)(i — iy) without vorticity, i.e., — ip 
and ^1 = ^0 = 0. This configuration corresponds to /3 = and j ~ a ~ 9, where 9 is the 
azimuthal angle. The adiabatic condition with respect to the change of Bz is required for atoms 
to remain in the weak-field seeking state, so that 1 is always antiparallel to B. Thus, when Bz 
is gradually changed in the opposite — z direction, the 1-vector points +z so that /S — tt, but 
7 = a = 6* is maintained by the effect of B^. Substituting these angles into and 5*0, we 
obtain = = and '^i — ipe~'^^^ , which corresponds to a vortex with a winding number 
q — 2. This result can be reinterpreted in terms of Berry's phase j205l 1304] . In general, when the 
hyperfine spin is F, we obtain a vortex with a winding number 2F since "^f^p and acquire 
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phases F{a — 7) and F{—a — 7), respectively [2061 1305] . 

When the bias field is turned off at i?^ = 0, a spin texture known as cross disgyration appears 
in the loffe-Pritchard trap. Here, the angle f3 increases from to tt/2, whereas -f = a = 9. The 
1- vector aligns with a hyperbolic distribution around the singularity at the center by following the 
direction of a local magnetic field of the loffe-Pritchard trap. This spin texture has been observed 
as a coreless vortex composed of three-component order parameters ^'±1,0 |268| . and its theoretical 
account will be given below. 

If one inverts the bias field back to the initial state after the vortex creation, the induced vortex 
unwinds itself, and the resulting state is non-rotating. On the other hand, if higher order magnetic 
fields are employed, it is possible to gain more than 2Fh of angular momentum per particle through 
a single inversion of the bias field. Hence, by using, e.g., a hexapole field in the radial direction 
and inverting the bias field B^, a 4_F-quantum vortex is produced. Switching to a quadrupole field 
and inverting the bias field back to the original value, one loses only 2Fh of angular momentum 
and the final state has a total angular momentum of 2Fh. This process can be repeated, and each 
cycle increases the angular momentum of the system by 2Fh per particle. The idea of cyclicly 
pumping vortices into a condensate was theoretically proposed in Ref . [3061 1307) . 

3.4.4 A single vortex (a few vortices) in spin-1 BECs 

Next, we consider the vortex states in slowly rotating spin-1 BECs. The characteristics of the vortex 
states are strongly dependent on the sign of gs, i.e., ferromagnetic {gs > 0) or antiferromagnetic 
{gs < 0). This can be seen by their superfluid velocity v = {h/m)(^W(; from Eq. p.3ip . they are 

,^fcrro ^ ry / _ i ^ ^polar ^ i^y^^ .332) 

m m 

where the angles are assumed to be spatially varying functions. For the ferromagnetic case, the 
supercurrent can be induced by not only a spatial variation of the phase Xi but also by that of the 
spin, similar to the pseudospin interpretation of two-component BECs described in Sec. 13.3.31 On 
the other hand, the polar case resembles that of scalar BECs. However, it allows half of the unit 
of quantization due to order parameter symmetry, as shown below. 

The structure and stability of the axisymmetric vortex states for various fl and total magneti- 
zations M were thoroughly investigated in Ref. |308[ 1309] . The equilibrium state can be obtained 
by minimizing the energy in the rotating frame: E — \M — Vl{Lz). The axisymmetric vortex 
configuration is classified by the combination of the winding number qj of the condensate wave 
function j = •y^nj(r)e*(^3+*J^) with the cylindrical coordinate {r,9,z), where the homogeneity of 
the wave functions along the z-axis is assumed, and Xj is the overall phase of the j-th component. 
The phases are constrained to minimize the spin-exchange interaction energy in Eq. (I3.28P : 

Es = ^|S|2 ^il{{n,~ + 2no [ni + n_i + 2V^ITHrTcos(x + qO)] } , (3.33) 

where x = Xi + X-i ^ 2xo and q = qi + q-i — 2qo. To minimize Es, Xj and qj should satisfy 
2xo = Xi + X-i + and 2qa = qi + q-i, where n' is an odd (even) integer for antiferromagnetic 
(ferromagnetic) interaction. The global phase x has no effect on the vortex structure, and was 
therefore set to zero. The possible combination of qj satisfying the above relation gives this 
system a characteristic vortex structure. If the values of qj are restricted to \qj\ < 1, we have 
(gi,goi9-i) = ±(1: 1,1)1 i(liOi~l) and ±(1,1/2,0), where cases with negative sign are omitted 
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in the following. For (1,1/2,0), the value 1/2 is not allowed in an axisymmetric system, so that 
the ^'o-component vanishes; we denote this state as (1, x,0). 

Unfortunately, the assumption of axisymmetry restricts the possible ground states, limiting the 
stable vortex configurations. The vortex states without assuming axisymmetry were considered by 
some authors [3101 13111 13121 13131 13141 1315) , who constructed phase diagrams of vortex states for 
various parameter spaces in various situations. The details for each vortex state are summarized 
below. 

(1, 1, 1) vortex 

In this case, all components have singly-quantized vortices at the center, and the densities of 
the three components are fully overlapped. As a result, the total density is equivalent that of a 
conventional BEC with an empty vortex core. 

With a fixed A'l , this axisymmetric vortex state cannot be stabilized in any parameter region for 
both (?s < and 5s > ,.309] . However, as in the (1,1) vortex state in two-component BECs shown 
in Fig. 13.51 by displacing the vortex cores of each component from the center of the trap, the (1,1,1) 
vortex can be stabilized in a nonaxisymmetric coreless configuration [312] . For > Q, the vortex 
cores are displaced such that the condensates have two singularities through an overlap of the 
and ^'-1 components, or they have three singularities that form a triangular configuration. For 
gs < 0, two singularities of '^i and are displaced from the center, while ^'o with a singularity 
at the center prevents phase separation. For all of these cases, the breaking of axisymmetry leads 
to a smooth variation of the total density by decreasing the overlap area of'^±i and 'So for > 0, 
and and for gg < 0. 

(1, 0, —1) vortex 

In this configuration, the non-rotating '^q component occupies the central region of the vortex core, 
which is made up of the ^'i and components with opposite circulations. This vortex state is 
referred to as a "polar-core" vortex. Since the condensate at the center of the trap consists only 
of the polar state, the spin texture has a singularity, although it can vary continuously around the 
vortex core. This vortex state is thermodynamically stable in the high magnetization region, which 
is independent of the strength and sign of the spin exchange interaction [312] . Also, it becomes 
stable in a non- rotating harmonic trap with an applied loffe-Pritchard magnetic field [314] . 

Interestingly, this vortex state can be generated by the intrinsic dynamical instability of a fer- 
romagnetic spinor BEC. Saito et al. demonstrated that an initially prepared without vorticity 
can generate transverse magnetization due to ferromagnetic interaction [316j . Since spin conser- 
vation prohibits the appearance of longitudinal magnetization, the dynamical instability results in 
a polar-core vortex by spontaneously breaking the chiral symmetry of the initial state, where the 
^±1 component can take either q — ±1 circulation. Since oppositely rotating VP^i components 
have the same amplitude due to M = 0, the vortex carries no net mass current but a spin current 
with one quantum of circulation. This spontaneous formation of such polar-core "spin" vortices 
was demonstrated experimentally by the Berkeley group, using an instantaneous quench from the 
polar phase to the ferromagnetic phase r270j [See Fig. 13. 8| . 
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(1, x,0) vortex 

This vortex state consists of a non- rotating component at the center fiUing the vortex core 
of the component. This is similar to the (1,0) vortex of two-component BECs, but its in- 
terpretation is quite difFcrcnt. In the (l,x,0) vortex, the spin component is suddenly reversed 
near the vortex core of the VPi component because the absence of the vPq component means that 
Sx = Sy — 0, which cannot yield a continuous variation of the spin texture. This vortex state is 
thermodynamically stable only for the polar case (gs > 0) |3091 1312] . 

This vortex state is called the half-quantum vortex state ("Alice state") |3091I281] . To see this, 
let us consider 'S'P°''*'' and a loop that encircles a vortex with a fixed radius. Each point on the 
loop can be specified by its azimuthal angle 9. The single- valuedness of the order parameter vI'P°''^'' 
in Eq. (I3.3ip is met if we take, for example, —a — x = 9^/2 and (3 — Tr/2. The order parameter at 
r — cx) is given by 'S'P°'^'' = y^iT/2(— e*'^, 0, 1). Thus, from Eq. (j3.32L the circulation becomes 
^.^poiar . ^£ _ (^fi/2m)q, quantized in units of h/{2m) rather than the usual h/m. The underlying 
physics for this half quantum number are based on the discrete Z2 symmetry of 'S'P°''''', which 
is invariant under tt gauge transformation (x ~^ X + combined with a tt spin rotation around 
an axis perpendicular to n = (cos a sin/?, sin a sin /3, cos /3): n — !■ — n. Thus, the polar phase of a 
spin-1 BEC can host a half-quantum vortex |299l 1281] . The dynamical creation of half-quantum 
vortices under an external rotating potential was demonstrated in numerical simulations j317l 1318] . 

The half-quantum vortex state is unstable for ferromagnetic spinor BECs. The axisymmetry 
of the transverse magnetization is spontaneously broken to form three- fold domains |319] . This 
originates from the topological structure of the half-quantum vortex and spin conservation. 

(0, 1, 2) vortex 

In ferromagnetic interactions, a more stable vortex configuration occurs when the winding number 
qj exceeds unity. Mizushima et al. found that ferromagnetic interaction supports the thermody- 
namic stability of the (0, 1, 2) axisymmetric vortex state over wide parameter regions |311[ 1312] . 
The central region of the harmonic trap is occupied by a non-rotating ^E*! component. The 
component with go = 1 is pushed outward, while the component with g_i = 2 occupies the 
outermost region. This configuration is favorable for the ferromagnetic case, because it is more 
favorable than spatial phase separation of 5'i and '5-1 ; the presence of vortices with different 
qj effectively causes phase separation in the radial direction. Only a (0,1,2) vortex can have a 
nonsingular continuous spin texture under slow rotation. Mizushima et al, performed extensive 
numerical studies to determine the fi-M phase diagram of the vortex states, including the (2, 1, 0) 
state referred to as the mixed- twist texture |315j . 

When the wavefunction is parameterized by the spatially-dependent bending angle /3(r) 

and X = and a = -7 = 61 in Eq. (15311) . we have = ^/m(r) (cos^ e*^V2 sin ^ cos 

g2ie gjj^2 where /?(r) is an increasing function of r from /3(0) = and < (3{r) < n. The 

spin direction is related to the 1- vector, and is given as l(r) — zcos/3(r) -|-sin/3(r)(cos0x-|-sin6'y), 
where f3 varies from /3{0) = to /3(i?) = tt/2 {— tt) for a Mermin-Ho (Anderson- Toulouse) vortex. 
As in two-component BECs [see Sec. 13.3.3] . the (0, 1, 2) vortex can be interpreted as a vortex with 
an intermediate boundary condition (tt/2 < P{R) < tt); we can control the value of /?(i?) from 
the Mermin-Ho condition to the Anderson- Toulouse condition by merely changing the total mag- 
netization M [311] . In another numerical study, Martikainen et al. analyzed the coreless vortex 
state as a function of rotational frequency without fixing the total magnetization [313] , finding that 
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P{R) increases with increasing ft and that the upper value of f3{R) is 37r/4, above which additional 
vortices nucleate. This implies that an Anderson- Toulouse vortex can never be the ground state 
of the system. 

Because this coreless vortex contains a doubly-charged vortex in one of the components, it is 
interesting to study whether the coreless vortex state inherits the dynamical instability of doubly- 
charged vortices in scalar condensates, as shown in Sec. 13.2.61 Pietila et al. studied the low-lying 
excitations of the (0, 1, 2) state for < in the presence of an loffe-Pritchard trap based on the 
BdG approach. The phase diagram of the dynamical stability was explored for the bias field Bz 
and the transverse quadrupole field i?^ |320j . Takahashi et al. considered the same problem for 
both f^s < and > without the loffe-Pritchard trap, finding that the dynamical instability is 
suppressed (enhanced) for < (g^ > 0) [321) . The unstable dynamics are associated with the 
vortex-splitting of the q = 2 vortex, and phase-separation in the gs > case. 

3.4.5 Vortex lattices in F=l spinor BECs 

Vortex lattices in rapidly rotating spin-1 BECs can possess richer vortex phases than those in 
two-component BECs, but their classification is very complicated. Kita et al. studied vortex- 
lattice structures of antiferromagnetic spin-1 BECs using the phenomenological Ginzburg-Landau 
equation |322j . which is formally similar to the mean-field quantum Hall approach. This study 
revealed that the conventional Abrikosov lattice with hard vortex cores is unstable, and the vortex 
cores shift their locations. The system has many metastable configurations of vortices, depending 
sensitively on the ratio gs/gn- The vortices are characterized by the distribution of magnetization 
and the difference in the number of circulation of quanta per unit cell, all of which makes the 
characterization of the vortex lattice structure complicated. 

Further discussion on the phase diagram of the ground state in rotating spin-1 bosons was 
reported in Reijnders et al. [323] . This paper also describes a study of vortex-lattice structures with 
the mean-field quantum Hall approach over wide parameter regions and the exact diagonalization 
study of the quantum liquid phase. A similar analysis was done by Mueller for slowly rotating 
spin-1 BECs [282j . In addition, Mizushima et al. performed numerical simulations of the CP 
equations p.30p for a ferromagnetic spinor BEC in a fast-rotating regime and also in an external 
magnetic field |315j . For M = 0, the equilibrium state is a square lattice constructed from two 
sublattices of the Mermin-Ho vortex and the mixed-twisted vortex, and the local spin on the two 
sublattice sites are locked in alternate directions z and — z. For M 7^ 0, the composite lattice of a 
coreless vortex and a polar-core vortex was found. 

3.4.6 Quenched spinor BEC: Kibble-Zurek mechanism 

Topological defect formation via the Kibble-Zurek mechanism can occur in a ferromagnetic phase of 
a spinor BEC. In an experiment by the Berkeley group |270j . ^^Rb atoms (with gs < 0) in the m — 
state were prepared by applying a magnetic field along the z-axis. This initial state corresponds 
to the point in region (III) of the phase diagram in Fig I3.8f b). The strength of the magnetic 
field was then suddenly decreased toward region (I) along a constant -p line of the phase diagram. 
The m = state was no longer the ground state, and a spontaneous transverse magnetization 
emerged due to the conservation of longitudinal magnetization. Since this phase transition was 
triggered by a sudden change in a magnetic field, this can be regarded as a quantum quench at 
r = 0. The magnetization dynamics were observed by a spin-sensitive in situ measurement to 
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form complicated ferromagnetic domains. Remarkably, polar-core spin vortices were identified in 
some snapshots of the spin distribution, as shown in Fig. I3.8f c). consistent with the prediction 

m- 

Theoretical analysis of this experiment was done by several authors [3241 13251 13261 13271 1328) , 
in relation to the Kibble-Zurek mechanism. The numerical simulations revealed that the creation 
of the spin vortices was caused by instability of the initially generated solitons, whose pattern 
and characteristic length scales were strongly dependent on not only the quench time but also the 
properties of the initial noise. 

3.4.7 Vortices in F=2 spinor condensates 

For the F — 2 case, the field operator has five components (%l)2,'^i,ipQi'^-i,'^-2)- The interaction 
becomes 

a,P=-2 a,Q',/3,/3' = -2 ^ 

+bc2i^i^{2, a; 2, a'|0, 0)(0, 0|2, /3; 2, • (3.34) 

The collision coefficients are given by co = (452 + 334)/?, ci = —(32 ~54)/7, and C2 = ((?o — 54)/5 — 
2(52 ~ 54)77. The energy functional E = (H) for spin-2 BECs can be written as 

Em = jdv^Y. + Y"T + ylSp + |l^2oP^ , (3.35) 

where nT(r) = Y^a=-2 I^q(i")P the total particle density and S — {{Fx), {Fy), {F^)) is the spin 
density vector defined by (^^,(r)) = El^=-2 *a(i-)(^^)a/^*/3(i-) {i = x,y,z). The new feature 
of spin-2 BECs, not found in spin-1 BECs, is the presence of the spin-singlet pair amplitude 
^20 = (25'24'-2 — 2^'i5'_i -t- \I'q)/'\/5. Figure [379r a) shows a mean-field phase diagram in the 
absence of a magnetic field. There are four distinct phases: ferromagnetic, uniaxial nematic, 
biaxial nematic, and cyclic, where the order parameters of the last three phases are characterized 
by a discrete symmetry. According to mean-field theory, the uniaxial nematic phase and biaxial 
phase are degenerate in the absence of a magnetic field, however, quantum fluctuations induce a 
quantum phase transition between the two phases, lifting the degeneracy [33011331] . For ci > and 
C2 > 0, both S and A20 vanish, which yields the cyclic phase concentration characterized by the 



spin-singlet "trios" peculiar to spin-2 BECs [332) ; one representation for the cyclic phase is given by 
^^=5"= = y/nT{i/2, 0, 1/V2, 0,i/2)^ . Although several experiments demonstrated that spinor BECs 
with an F = 2 manifold of ®^Rb atoms were characterized by biaxial-nematic spinor dynamics, 
[1451 1147j , the possibility of the cyclic phase has not yet been excluded due to complications arising 
from quadratic Zeeman effects and hyperfine-spin-exchange relaxations [333] . 

The vortex states in spin-2 BECs are strongly affected by the symmetry of the order parameters 
in each phase. For ferromagnetic phases, the properties of the vortex states are not so different 
from those of spin-1 BECs. However, the symmetry of the biaxial nematic phase and the cyclic 
phase can yield the exotic vortices governed by the non- Abelian group [3341 13351 13361 13371 1338) . 
Here, we focus on the vortex states in the cyclic phase, and briefly summarize the salient results. 



65 



(a) 



Uniaxial Nematic: 

W' = ^iij (0. 0. 1. 0, 0)' 



Biaxial Nematic: 

W"=-^^(1.0. 0. 0. l)7V2 



Cyclic: 

C2 



Ferromagnetic: 

If""" = 7^(1, 0, 0, 0, 0)' 




Figure 3.9: (a) Mean-field phase diagram of spin-2 BECs. (b) Scfiematic illustration showing 
that the symmetry axes of the cyclic order parameter constitute a tetrahedron, (c) Numerical 
simulation of the collision dynamics of non-Abelian vortices. After the collision, two vortices get 
connected and a rung appears between them. The topological charges for each vortex line are 
denoted [Kobayashi et al: Phys. Rev. Lett. 103 (2009) 115301, reproduced with permission. 
Copyright 2009 by the American Physical Society]. 



Vortices in cyclic condensates 



As in the spin-1 BEC, through the U{1) gauge transformation and the 5*0(3) spin rotation, it 



is possible to transform from one to another cyclic state, 'Jf'^y^' — 



ne^fcy<=, where n and 



are the unit vector of the rotational axis and angle of the spin rotation, respectively. The order 
parameter vP'^y'^ is invariant under the following 12 transformations: 1, Ix — e^^" 



ly — C 



F^+F,+F^)/3V3^ C\ i^c, lyC, hC, hC\ IyC\ and hC^ [SM]. The 
overbar was added to emphasize that the operation includes not only a rotation in spin space but 
also a gauge transformation. These 12 transformations form the non-Abelian tetrahedral group T, 
shown in Fig. 13.91 (b). 



Topological charges of vortices can be classified by 12 generators 1, Ix, ly 



of the non- 



Abelian group T. Vortices are also classified into four conjugacy classes: (I): integer vortex; 1, 
(II): 1/2 - spin vortex; Ix, ly, and 4, (III): 1/3 vortex; C, IxC, lyC, and hC, and (IV): 2/3 
vortex; C^, IxC^, lyC^, and IzG^ ■ Topological charges in the same conjugacy class transform into 
one another under the arbitrary global gauge transformation and the spin rotation, denoted by S. 
The order parameters for straight vortices along the z-axis in each conjugacy class can be written 
in cylindrical coordinates (r, 9, z) as 



(i/(r)e2*("^+i)^ 0, ^/2h{r), 0, i/(r)e-2*("^+i)^ 

(z/(r)e^(2„. + l)e^ 0^ ^^^^y-^i^2n, + l)e 

1_ (/(,.)e^(2n.+i)e, 0, 0, y2g(r)e— ^^ 0' 
A (/(r.)e*(2»^-i)^ 0, 0, V25(r)e-"^^ 



(3.36) 



for (I), (II), (III), and (IV), respectively, where the vortex is located at r = 0. Here, ni and n2 
are the integer winding numbers, and /, g, and h are real functions that satisfy [/^ -I- /i^]/2 = 
[P + 2(7^]/3 — 1 and /(r — > 00) — g[r — > 00) — h{r — > 00) — 1. At the vortex core, the cyclic 
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order parameter changes to that of a different phase. With the minimum windings (rii — 712 = 0), 
we obtain the core structure of each conjugacy class by taking /(r = 0) = as: (I) (II) ^' oc 
S'(0, 0, 1, 0, 0)^ and (III)(IV) * cx S{0, 0, 0, 1, 0)'^, i.e., the core of (I) and (II) vortices have a finite 
spin-singlet pair amplitude { A20 7^ 0), and that of (III) and (IV) vortices have a finite magnetization 
{Sz 7^ 0). The 1/3 and 2/3 vortices are a hallmark of the cyclic state with tetrahedral symmetry. 
The energetic stability of these fractional vortices in rotating potentials was examined in Ref. [338] 

A salient feature of non-Abelian vortices emerges in their collision dynamics |337) . When two 
Abelian vortices collide, there are three possibilities: (i) they reconnect themselves, (ii) they pass 
through each other, or (iii) they form a rung that bridges the two vortices, depending on the 
kinematic parameters and initial conditions. Usually, reconnection occurs due to the energetic 
constraint. However, when two non-Abelian vortices collide, only a rung can be formed; recon- 
nection and passing through are topologically forbidden because the corresponding generators do 
not commute with each other. In fact, the nonzero commutator of the two generators gives the 
generator of the rung vortex. Figure I3.9f c) illustrates a typical rung formation. 

In other studies, vortex molecules consisting of a bound pair of fractional vortices was discussed 
in Refs. [3361 1338j . The structure of vortex lattices in cyclic condensates was discussed in Ref. 
|335j . where the transitions of the various lattice geometries were demonstrated as a function of 
applied magnetic field and temperature. 



3.5 Vortices in dipolar condensates 

A dipolar condensate refers to a condensate in which atoms interact via dipole-dipole interactions 
(DDIs), in addition to the usual s-wave contact interactions. For two particles 1 and 2 with dipole 
moments along the unit vector ei and 62, and relative position r, the DDI has the form 

Cdd ei -62-3(61 ■f)(e2 -r) 
^dd(r) - — . (3.37) 

The coupling constant Cdd is given as /ioMd particles with a permanent magnetic dipole moment 
/Xd (/xo is the magnetic permeability of vacuum) and (P/eo for a permanent electric dipole moment 
d (eo is the permittivity of vacuum). Recently, BECs of chromium atoms have been created [339], 
in which the atoms exhibited a larger magnetic-dipole moment (/id — S/ie) than typical alkali 
atoms {fi(i ~ /ib). This enables the study of the effects of anisotropic long-range interactions in 
BECs. The rapid progress in studies of dipolar BECs can be seen in the review article of Ref. 

[313]. 

In particular, when the dipole moments are polarized under an external field along the z- 
axis, the interaction potential between two magnetic dipoles with separated by r is given by 
Vdd(r) = (/ioMd/47r)(l — 3 cos^ 0)/r^, where Bz-r = cos 6*. Thus, the system is described by the scalar 
condensate wave function, and the DDIs contribute to the CP equation as a nonlocal mean-field 
potential 

K=x + 5l*l' + / dr'ydd(r'-r)|vI/(r')P 



^9* 



2m 



(3.38) 



Since the contact interaction g — AirfPa/m can be tuned to zero by the Feshbach resonance 
technique, we can obtain novel quantum ferrofluids dominated by the dipole-dipole interaction 
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Before studying the vortex states, it is useful to examine the properties of non-rotating dipolar 
BECs. The principal effect of Vdd on the equilibrium shape of a trapped condensate is to cause 
distortion of its aspect ratio |340[|341] . Since the dipoles are aligned along the z-axis, the attractive 
interaction of dipoles becomes significant in a cigar-shaped trap, and the atomic cloud is elongated 
along the z-axis. Although this may appear counterintuitive, it can be understood by considering 
the energy density of the DDI. It has the anisotropic form of a saddle, with negative curvature 
along the direction of the dipoles |342j . Because of attractive interactions, the condensate will 
eventually collapse, once the particle number exceeds some critical value Nc, similar to conden- 
sates with attractive contact interactions |155j . For a pancake trap, the condensate is elongated 
similarly along the z-axis, but stable against the collapse because the DDI is mostly repulsive. 
The instability condition also depends on the strength of the contact interaction. These collapse 
dynamics were demonstrated by tuning the s-wave scattering length into the unstable regime, 
showing an anisotropic explosion characterized by d-wave symmetry |344j. 

Another interesting feature of a dipolar BEC is that the Bogoliubov excitation spectrum exhibits 
a "roton-niaxon" behavior, caused by the attractive nature of the DDI. Santos et al, found that, 
for a quasi-2D system, harmonically confined in the z-direction and free in the x and y directions, 
the excitation energy has a local maximum and minimum dependent upon the in-plane momentum 
q [3451 . Although this behavior is reminiscent of that found in the excitation spectrum of liquid 
helium II, its physical origin is somewhat different. This can be intuitively understood as follows. 
When the in-plane momenta q are much smaller than the inverse size L of the condensate in the 
z direction, excitations have a 2D character. Because the dipoles are oriented along the z-axis, 
the interaction is effectively repulsive and the in-plane excitations become phonons. For q <^ 
excitations begin to have a 3D character and the interparticle repulsion is reduced due to the 
attractive force in the z direction. This decreases the excitation energy with an increase in q. The 
excitation energy reaches a minimum (roton) and then begins to increase with an further increase 
in q, eventually leading to a quadratic behavior oc q. 

As the dipole interaction becomes stronger, the energy at the roton minimum decreases, and 
then reaches zero, representing an onset of instability. This "roton instability" is probably towards 
the formation of a density wave. Ronen et al. numerically studied the equilibrium solutions 
for a pancake geometry, finding that the condensate always has unstable criterion of the dipolar 
interaction strength, even for the pancake-trap limit [346] . They observed that, near the instability 
threshold in a certain range of the trap aspect ratio, the equilibrium density possesses "biconcave" 
shape, where the density modulates radially and its local minimum appears at the center. This 
density modulation is associated with the roton instability. Note that this physical origin of the 
instability is different from that discussed in the dipolar collapse experiment; the roton instability 
exists even for a 2D pancake geometry. A more detailed discussion and references can be found in 
the review article |343) . 

3.5.1 Vortices in spin-poralized dipolar BECs 

Here, we summarize some of the theoretical results on vortices in a spin-poralized dipolar BEC, in 
which the above features of the excitation spectrum have a strongly influence on vortex states. 
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Structure of a single-vortex 

The effect of dipolar interactions on the single vortex state was first studied by Yi and Pu by 
numerically solving Eq. p.38p |347j . While there was no significant difference of the vortex 
structure for the repulsive contact interaction, the density profile exhibited a "biconcave" structure 
around the vortex core when the condensate had attractive contact interactions (a < 0) and was 
trapped in a highly prolate trap. More detailed numerical study was performed by Wilson et al. 
|348| , who found that these biconcave density ripples around the vortex core emerge at milder trap 
aspect ratios and with purely dipolar interactions (a = 0). By using perturbation theory, they 
related these density oscillations to the roton mode near the instability. 

Vortex stability and dynamics 

The thermodynamic critical rotation frequency flc for dipolar BECs in the TF limit was investi- 
gated by O'Dell and Eberlein [349]. The value of ilc decreased for a condensate in a pancake-shaped 
trap, while it increased in a cigar-shaped trap, compared to that of a conventional BEC. This is 
because, for the pancake trap, the DDI is repulsive on average, and results in the slight increase 
of both the axial size Rz and the radial size Rj_. Thus, fie decreases since oc i?]]^- For the 
cigar-shaped trap, the DDI caused the slight increase in the axial size Rz but decreased the radial 
size R±, which increased Vic- The TF approximation is invalid for pure dipolar BECs. The full 
numerical calculations for the single-vortex state were done by Abad et al. 350 , and their results 
were good agreement with Ref. |349| . 

Another important issue is the critical rotation frequency for vortex nucleation in rotating 
dipolar BECs. Following the analysis for the conventional BEC [187) . Bijnen et al. studied the 
dynamical instability of a rotating dipolar BEC in the TF limit |351l 1352] . The dynamically 
unstable region of the rotation frequency O versus trap aspect ratio w^/wj, nontrivially depends 
on the dipole interaction strength. Interestingly, the critical frequency fic for a cigar-shaped trap 
can become larger than the onset of the dynamical instability of a rotating condensate. This is an 
intriguing regime, in which a rotating dipolar BEC is dynamically unstable but vortices will not 
enter. 



The BdG analysis for a single vortex state in dipolar BECs was done by Wilson et al. |353j 



As discussed in Sec. 13.2.41 the single vortex state in a conventional BEC (without rotation) is 
thermodynamically unstable, but dynamically stable. However, the vortex states in dipolar BECs 
can be dynamically unstable, depending on the dipolar strength and the aspect ratio of the trap 
potential. For a pancake trap, the instability is associated with the radial and angular roton 
modes, which induce a density wave and local collapse. For a cigar-shaped trap, the 3D character 
of the DDI becomes important. Thus, the vortex wave (Kelvin mode) is strongly affected by the 
DDI. Klawunn et al. showed that, when applying a periodic potential along the vortex line, which 
changes the effective mass along this direction, the dispersion of the Kelvin mode has a roton-like 
minimum |354j . As the dipolar interaction increases, this minimum can reach zero, forming an 
instability similar to the Donnelly-Glaberson instability (see Sec. 13.2.61 for the discussion of vortex 
waves). 

Vortex lattices 

Rapidly rotating dipolar BECs exhibit a rich variety of vortex phases characterized by different 
symmetries of the lattice structure |355l 13561 1357] . Cooper et al. treated this problem within the 
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LLL approximation |355) . In particular, with increasing dipole interaction strength or with a high 
fiUing factor, the vortex lattice may undergo transitions between different symmetries: triangular, 
square, striped vortex crystal, and bubble state. In addition, for vortex lattices in double-well 
potentials, the competition between tunneling and interlayer DDI should lead to a quantum phase 
transition from a coincident phase to a staggered phase [356) . 

3.5.2 Vortices in spinor dipolar condensates 

When the spin degree of freedom is taken into account for dipolar BECs, extremely rich physics 
appear. Let us consider the F — 1 case. Typically, the energy associated with spin-exchange 
interactions is much smaller than that of contact (spin-preserving) interactions, i.e., gn ^ gs- 
Hence, the spin-exchange interaction may become comparable to the DDI. As a consequence, even 
for alkali spinor BECs (in particular ^^Rb) the DDI plays a significant role [358) . 

Under a uniform magnetic field B, the second quantized Hamiltonian for spin-1 dipolar BECs 
is H = Hq + Hint + ^^dd, where Ho and Hint = f dridr2VF=i(ri — is given by Eqs. (|3.25p and 
(j3.27p . respectively. The dipolar contribution becomes 



where ri2 — (r — r')/|r — r'| is a unit vector. The spin-exchange term and the dipolar term describe 
two types of spin-dependent interactions. Contrary to the contact interaction, the DDI does not 
necessarily conserve the spin projection along the quantization axis due to the anisotropic character 
of the interaction. Note that the direction of the spin of the order parameter is not polarized by 
an external field. In such a situation, the system can develop a nontrivial spin texture due to the 
interplay and competition between these two terms. 

Einstein-de Haas effect 

If the atoms are initially prepared into a maximally polarized state, e.g., mp = —F, contact 
interactions cannot induce any spinor dynamics due to the conservation of total magnetization 
M . DDIs, on the contrary, may induce a transfer into an uip + 1 state. Here, the DDI conserves 
the total (spin -1- orbital) angular momentum of the system. If the system preserves cylindrical 
symmetry around the quantization axis, this violation of the spin conservation is accompanied by 
a transfer of angular momentum from the spin to the orbital part. Due to this transfer, an initially 
spin-polarized dipolar condensate can dynamically generate vortices in a process resembling the 
Einstein-de Haas effect [55^1550] . 

Unfortunately, the Einstein-de Haas effect is prevented in the presence of even a weak magnetic 
field, which favors the spin polarization. However, for a ferromagnetic spinor BEC such as _F = 1 
^^Rb BECs, a significant Einstein-de Haas effect may be expected under specific resonant conditions 
|361j . The population transfer away from the initially prepared mp = 1 state is typically very 
small, but can be significantly enhanced by applying a resonant magnetic field, such that the 
linear Zeeman energy of an atom in the mp — 1 state is totally transferred into the kinetic energy 
of the rotating atom in mp = 0. 




1 



37/;l(r),AUr')(F„„. •fi2)(F/3^. •ri2)^/3'(r')^a'(r) 



(3.39) 
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Spontaneous circulation in the ground state 

The ground state structure of the dipolar spinor condensate was first studied by Yi et al. 1362] 
under the single mode approximation (SMA), which assumes that all components share the same 
spatial wave function. For spinor BECs, the SMA is valid when the spin-dependent interaction 
is sufficiently weaker than the spin-independent interaction ((;„ 3> Qs for spin-1 case). However, 
in the presence of even weak DDIs, the dipole moments are spatially modulated due to the long- 
range anisotropic nature of the interactions, thus the application of the SMA is questionable. It 
is interesting that, in this system, the spatial spin variation, caused by the DDI, induces the 
supercurrent due to the spin-gauge symmetry, as seen in Eq. (I3.32[) . 

Systematic numerical studies of the ground state are described in Refs. |363l [364], revealing 
the spontaneous emergence of spin textures in the ground state. Yi and Pu constructed a phase 
diagram of the ground state of a spin-1 dipolar BEC with respect to the trap aspect ratio A and the 
dipolar interaction strength Cdd, for both the ferromagnetic case and the antiferromagnetic case 
|363] . For Cdd smaller than gn, the SMA is valid. Otherwise, three distinct spin textures appear, 
which can be classified as: (a) for a pancake geometry (A > 1), the polar-core spin vortex with 
(—1,0,1). (b) for a cigar-shaped geometry (A < 1), the coreless vortex with (0,1,2), where the 
spins twist around the z-axis. (c) for a spherical geometry (A ^ 1), the non-axisymmctric state. 
Kawaguchi et al. studied ferromagnetic spin-1 BECs by varying gg and Cdd, finding three distinct 
phases: a polar-core vortex, a flower phase, and a chiral spin- vortex phase [364] . Hence, dipolar 
spinor BECs can have spontaneous circulation in the ground state, even without rotation. 

3.6 Quantum turbulence 

As discussed in the previous chapter, the study of the turbulent state in quantum fluids and its 
relation to classical turbulence (CT) is an intriguing physical problem. Although the study of 
quantum turbulence (QT) has a long history, only superfluid "^He and '^He systems have been used 
to realize QT until recently. Recently, atomic Bose-Einstein condensates have become another 
candidate for QT research, since a turbulent state was realized in this system. 

In this section, we briefly review some of theoretical studies and an important experimental 
work for QT in atomic BECs. 

3.6.1 Theoretical work on quantum turbulence in atomic BECs 

A major problem in the study of QT in atomic BEC is the difficulty of applying a velocity field to 
generate the turbulent state. Thus, alternative methods are required to realize QT. 

BerlofF and Svistunov suggested the realization of QT in the dynamics of the formation of a BEC 
from a strongly degenerate nonequilibrium gas [3661 1367) . Such dynamics can also be described by 
the CP equation in the framework of a classical field description. 

Parker and Adams suggest the emergence and decay of turbulence in an atomic BEC under a 
simple rotation, starting from a vortex- free equilibrium BEC 1194] . Starting from the vortex- free 
steady state with a constant potential, they performed a numerical simulation of the CP model 
with realistic experimental parameters for ^^Rb BECs, and showed that the total dynamics can be 
divided into four distinct stages as follows: (i) Fragmentation - the quadrupolar oscillation of the 
condensate breaks down, ejecting energetic atoms to form an outer cloud, (ii) symmetry breaking 
- the twofold rotational symmetry of the system is broken, allowing the rotation to couple to 
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additional modes, thereby injecting energy rapidly into the system, (iii) turbulence - a turbulent 
cloud containing vortices and high energy density fluctuations (sound field) is formed, and (iv) 
crystallization - the loss of energy at short length scales coupled with vortex-sound interactions 
allows the system to relax into an ordered lattice. In the third region of turbulence, the system 
produces a Kolmogorov energy spectrum (|2.10l) demonstrating a classical analog of QT. 

To realize steady QT, Kobayashi and Tsubota suggested combining rotations around two axes 
[368]. Starting from a vortex- free initial state, they performed numerical simulations of the GP 
model with realistic parameters for ^''Rb BEG, and obtained a steady turbulent state with no 
crystallization, but with highly-tangled quantized vortices. The incompressible kinetic energy 
spectrum satisfies the Kolmogorov law in the inertial range. 

White et al. numerically investigated QT in atomic BEGs, and found nonclassical velocity 
statistics |369j . The probability density function (PDF) of velocity is another important statistical 
measure to study turbulence, and for the case of GT, the PDF obeys a Gaussian distribution 
[3701 13711 I372| . They performed numerical simulations of the GP model (|2.5|) using realistic 
experimental parameters for a ^"^Na condensate. The method to create turbulence was "phase- 
imprinting" in the condensate. They calculated the PDF of the superfluid velocity, and found 
that the velocity statistics were non-Gaussian and had a power-law dependence. Their result 
is consistent with the high-velocity tails found experimentally by Paoletti et al. in turbulent 
superfluid ''He [UTT) . 

3.6.2 Experimental work on quantum turbulence in atomic BECs 

Recently, a turbulent state was realized in atomic BEGs by two methods. Weiler et al. performed 
a rapid quench of an *^Rb gas through the BEG transition temperature |258j . which is a similar 
situation to that examined by Berloff and Svistunov |367| . Through the high density fluctuation 
regime (weak turbulence) in a short period, several vortices and anti-vortices were formed to make 
the system turbulent. This experiment was interpreted in the previous section fsubsection 13.2.8)) . 

The turbulent state created in the above method strongly depended on the initial uncontrollable 
thermal state. Furthermore, the turbulent state is merely intermediate, and the final state that 
contains only a few vortices is not QT. As a method with better control of the turbulence, Henn 
et al. introduced an external oscillatory perturbation to an ^^Rb BEG [3741 1375] . This oscillating 
magnetic field was produced by a pair of anti-Helmholtz coils which were not perfectly aligned to 
the vertical axis of the cigar-shaped condensate. Additionally, the components along the two equal 
directions that result in the radial symmetry of the trap were slightly different. This oscillatory 
field induced a coherent mode excitation in a BEG. For small amplitudes of the oscillating field 
and short excitation periods, dipolar modes, quadrupolar modes, and scissors modes of the BEG 
were observed, but no vortices appeared. Increasing both parameters, vortices grew in number, 
eventually leading to the turbulent state. In the turbulent regime, they observed a rapid increase 
in the number of vortices followed by proliferation of vortex lines in all directions, where many 
vortices with no preferred orientation formed a vortex tangle (Fig. 13.101) . Another remarkable 
feature is that a completely different hydrodynamic regime was followed: the suppression of aspect 
ratio inversion during free expansion, despite the asymmetric expansion (from a cigar-shaped to 
pancake-shaped) of the usual quantum gas of bosons, or the isotropic expansion of a thermal 
cloud. Although the theoretical understanding of this effect remains incomplete, it represents a 
remarkable new effect in atomic superfluids. The emergence of QT with vortex reconnections. 
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excitation of Kelvin waves, and the existence of cascade phenomena are expected to be directly 
observed in this well-controllable QT system in the near future. 




200/im 



Figure 3.10: Atomic optical density after 15 ms of free expansion (left) and a corresponding 
schematic diagram showing the inferred distribution of vortices obtained from the image (right) 
[575] [Henn, Seman, Roati, Magalhaes, and Bagnato: Phys. Rev. Lett. 103 (2009) 045301, 
reproduced with permission. Copyright 2009 by the American Physical Society]. 
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Chapter 4 

Summary and conclusions 



In this chapter, we have reviewed recent topics on quantized vortices in superfluid heUum and 
atomic Bose-Einstein condensates. Quantized vortices were discovered in superfluid ^Hc in the 
1990s. However, they have recently grown in importance, for two reasons. The first reason is that 
the research of QT entered a new era since the mid 1950s. A primary interest is to understand the 
relationship between CT and QT. The second reason is the realization of atomic BECs. Modern 
optical techniques have enabled the direct visualization of quantized vortices, and multi-component 
BECs have further enriched the world of quantized vortices. 

The research of vortices in classical fluids has a very long history. When Leonardo da Vinci 
drew his famous sketches of water turbulence in the 16th century, he certainly recognized that 
turbulence consists of many eddies of different sizes. It was in 1858 that Helmholtz proposed the 
idea of a vortex filament. However, it is not straightforward to understand the role of vortices 
in the classical hydrodynamics of turbulence. One reason is that vortices in classical fluids are 
unstable and not very well defined. Quantum condensed systems create more stable and well- 
defined vortices, namely quantized vortices. 

This field is currently very active, attracting the interest of scientists from many different fields. 
We hope that this chapter will contribute to further development and progress in this fascinating 
field of physics. 
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